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1. — Introduction. 


The question of measurability of the electromagnetic field strengths was 
first, raised by LANDAU and PEIERLS (1) who concluded that the measurability 
of electromagnetic field strengths is subject to limitations which go beyond 
the assumptions of quantum electrodynamics, and that quantum electro- 
dynamics gets away from its physical foundations.  BoHR and ROSENFELD (?) 
in their eminent work pointed out that Landau and Peierl’s result does not 
necessarily affect the consistency of the theory because in studying the meas- 
urability of electromagnetic field strengths one has to use test bodies of finite 


(*) A thesis submitted to the Faculty of the University of North Carolina in 
May, 1960 in partial fulfillment of the requirements for the degree of Doctor of Phi- 
losophy in the Department of Physics. 

(1) L. Lanpavu and R. Priers: Zeits. f. Phys., 69, 56 (1931). 

(2) N. Bour and L. RosenreELD: Det Kgl. Dans. Videnskab. Selskab, 12, 8 (1933). 
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distributions of charge, and not point charges; 1.e., not the field components 
taken at definite space-time points as treated by LANDAU and PEIERLS. Ac- 
tually quantum electrodynamics is itself independent of the atomic structure 
of matter, as is apparent from the fact that Planck’s constant and the velo- 
city of light are the only constants which are insufficient to fix space-time 
dimensions appearing in the commutation relations of electromagnetie field 
strengths. Only by introducing the rest mass of the elementary particle can 
one fix the space-time dimension. By using extended test bodies whose linear 
dimensions are chosen large enough compared to atomic dimensions, one can 
regard their charge density to be approximately constant. Then the radiation 
reactions of the test body on itself, which would be uncontrollably large if 
we treated the problem as LANDAU and PEIERLS do, will now be almost neg- 
ligible. As emphasized by BoHR and ROSENFELD, the passage from classical 
field theory, which is characterized by the measurability of field strengths at 
definite space-time points by means of a point charge, to quantum electro- 
dynamics, in which the field strengths are represented by the average values 
of the idealized field strengths over space-time regions, corresponds to the 
transition from a point charge to an extended one. 

Some progress has been made during the past few years in the theory of 
quantization of the gravitational field, such as linear approximation by 
LAURENT (?) and GUPTA (*), Feynman quantization by MISNER and ROSEN (5) 
and quantization by the canonical formalism by Dirac (5). We shall study 
the question of measurability of gravitational field strengths in the case of 
the linearized theory, 7.e., consider a Lorentz space and treat the actual metric 
of the Riemannian manifold deviating from this flat case as small. In this 
theory the commutation relations of the quantized gravitational field strengths 
are c-numbers containing three universal constants—namely, the universal gra- 
vitational constant G, Planck’s constant h, and the velocity of light e. Similar 
to the case in quantum electrodynamics, one characterizes quantum gravito- 
dynamics by the quantized gravitational field strengths, t.e., the quantized 
Riemannian tensors, which are no longer represented by means of definite 
Space-time point functions but by means of functions averaged over space-time 
regions. Therefore, we are not concerned with the atomic structure of matter. 

It is significant that the usual description of a gravitational field in terms 
of the Riemann curvature tensors at definite Space-time points, which are 
characteristic of classical general relativity and which are measurable by means 
of an idealized point particle, cannot directly apply to quantum theory. By 


(*) B. E. LAURENT: Nuovo Cimento, 6, 1445 (1946). 

(4) S. N. Gupra: Rev. Mod. Phys., 29, 334 (1957) and papers cited there. 

() C. W. MISNER: Rev. Mod. Phys., 29, 497 (1957); G. Rosen: Ph. D. Thesis 
(Princeton, 1959). 


(9) P. A. M. Dirac: Proc. Roy. Soc. London, A 246, 333 (1958). 
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passing to quantum gravitodynamies in which the quantized gravitational field 
strengths, or the quantized Riemann curvature tensors, are no longer regarded 
as finite space-time point functions, but as functions of space-time regions, 
one finds the way to avoid the difficulties associated with a point mass, such 
as the infinite self energies of elementary particles inherited from classical 
field theory. 

Hence, the purpose of the present work is to study whether the limitations 
on the measurability of gravitational field strengths agree with those deter- 
mined by physically possible measurements so as to give a basis for the validity 
of quantizing the gravitational field. 

In a manner analogous to the electromagnetic formalism, the Riemann 
tensors can be divided into two classes: The gravitational field strengths of 
the electric-type components and those of the magnetic-type components. In 
the measurement of field averages of the magnetic-type components, one has 
to measure the momentum imparted to a rigid test body. Since the meas- 
urement is very similar to that in quantum electrodynamics, the magnetic-type 
components are to be averaged over a space-time region in the usual sense. 
However, in the measurement of field averages of the electric-type components, 
for which one requires a quite different type of physical object—the meas- 
urement of the strain of a test body, the field averages over a space-time region 
are no longer averages in the ordinary sense, but are weighted averages over a 
space-time region. 

If one ignores all limitations inherent in the atomic structure of the test 
body, then despite two kinds of field averages required in the mesurement 
of gravitational field strengths, a complete agreement between the limitations 
on the measurability of field averages and the quantum gravitational formalism 
can actually be shown. 

Similar to quantum electrodynamics, in which the formalism is understood 
in perturbation expansions in power series of the fine structure constant e*/he, 
one has in the quantum gravitational formalism expansions in power series 
of the constant Gm2/he, where m denotes the mass of an elementary particle. 
In the first approximation, therefore, the formalism is understood by quan- 
tizing a given pure gravitational field, and not by quantizing the interacting 
fields between a test body and the given gravitational field. The second ap- 
proximation, which is not treated in this paper, is the quantization of the 
field associated with mass particles. The radiative corrections (*), which do 
not essentially affect the conclusions attained in the first approximation, can 
be neglected. Their neglect can be justified in this case much better than in 
the case of quantum electrodynamics because the constant @m?/he is very 


(*) The radiative correction in the case of the quantum electromagnetic formalism 


has been treated by E. CORINALDESI (7). 
(7) E. CorinaLpEsI: Suppl. Nuovo Cimento, 10, 83 (1953). 
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small: of the order of magnitude 10-*° taking the value of the mass of a 
nucleon, d.e., m~10-** g. 

In Section 2 Einstein’s theory of general relativity in the linear approxi- 
mation is reviewed and commutators for the quantized field components as 
definite space-time point functions are given. In Section 3 some treatment 
of the quantum fluctuations is carried out and commutators between field 
averages are explicitly written down. Assumptions required for a test body 
in order to measure the electric-type components and a detailed analysis of 
the measuring process are given in Section 4, while a similar treatment for the 
magnetic-type components is given in Section 5. The uncertainties of field 
components in one region induced by the measurement of an electric-type com- 
ponent in another region are treated in Section 6 and those induced by the 
measurement of a magnetic-type component are discussed in Section 7. In 
Section 8 the measurement of a single field average is given. This is extended 
to the measurement of two field averages of the same type in Section 9 and of 
different types in Section 10. Finally, Section 11 is devoted to the study of 
the measurability in the more general case where three respective principal 
axes of two cubical test bodies are not parallel to each other. In every case 
the agreement between the quantum gravitational formalism and the limi- 
tations imposed on the measurability of the test body is attained successfully. 


2. — Commutators for the quantized gravitational field components in the linear- 
ized theory. 


Even though several theories of gravitation have been proposed, the most 
widely accepted is the theory of general relativity by Einstein. One denotes 
the quadratic form of the line element of the 4-dimensional Riemannian mani- 
fold by dS°=g,,dx“dx", where the Greek indices take the values 1, 2, 3, 4 
and the space-time co-ordinates are denoted by L = (X15 La, #3, ict), Einstein’s 
classical field equations in the presence of matter are given by 


ada K 
(2.1) RS eB 


MV 9 


7 a ara Dh a OS d 

where Ro are the Ricci tensors, R= lee the curvature sealer, Li the energy- 
momentum stress tensors of matter, G the universal gravitational constant, 
and ec the velocity of light. Eq. (2.1) can also be written in the form 


(2.2) Rpg STES T = Te 


where the usual summation convention is applied to the indices which appear 
twice. 
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As one knows that the field equations are non-linear with respect to the 
field variables Ius it 18 possible to linearize the theory with a considerable degree 
of accuracy. Hereafter we will limit our discussion to the linearized theory 
unless otherwise stated. In the linearized theory we shall assume that the 
deviations of the actual metric from that of flat space are so small that the 
actual metric can be approximated by Inv = Ou Pur Where ò,, is the Min- 
kowskian metric tensor which has the components 


Pelee} \ 


(2.3) ee 


and g,, are small quantities of the first order. The Riemann tensors are given by 


| en 


(2.4) RI 


pvot 2 (ges VE P ur vo uri Pro,ur in Oo) ’ 


where the commas followed by indices mean ordinary derivatives. It is worth- 
while to note that eq. (2.4) is invariant under a gauge transformation 
Po = Pun + SN which is induced by co-ordinate transformations of the 
type x, =, €,(,) which change the co-ordinate values by small quantities 
of the first order. 

By repeatedly using the Bianchi’s identities 


(2.5) li ey xg LEE: SF Res =:0 
and referring to Einstein’s eq. (2.1), one obtains 


== K 
(2.6) Roue ABS he re + Ove He ar (Tr 5 One 4 WS = 


ri = 5 Our SERRE A VS > Ope DO ’ 


the solution of which is given by 


= 


È Ja Di va (Ta F0, T), no A IE, ee 5 Ove T),ur Fa 


(2,7) RvoNT) nt o. 
Era (Ti 5 Our 14 ln LI FRE 4 do DIA ’ 


where D(x — x’), the Green’s function of the two space-time points # and 2’, 


satisfies the equation 


(2.8) C2 Dia — a’) = O(a — x'). 


In the absence of matter, the classical linearized field equations can be 


written in the following form (*) 


(2.9) Ry =0, D es 


(*) Latin indices take the values 1, 2, 3. 
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As a preliminary for quantization, it is more convenient to divide the grav- 
itational field strengths into two parts; namely, the electric-type compo- 
nents £,, and the magnetic-type component H,;, which are defined by 

v 


Hi => 


Li 


È Eirt 
(2.10) JD) rao ence Er Rei = — mae Rrvio ; 
pe] 


where £,,, is the 3-dimensional permutation symbol. It is worth-while to men- 
tion that the Riemann tensors are symmetric, i.e., H,,=0, E,;=E;;, and 
that by eq. (2.9), E.,=0 and H,,=H;. 

So far we have limited our discussion to classical general relativity theory. 
By applying quantum theory to general relativity, the quantized Riemann 
tensors, or the quantized gravitational field strengths, are no longer regarded 
as c-numbers, but are regarded as operators which operate on a state vector. 
The quantized Riemann tensors at two space-time points (x, y, 2,¢) and 
(2°, y’, 2, t') (*) no longer commute. The commutators have been given by 
Dr. B. S. DEWITT (8) as 


(2.11) nat), Lega (©) ue Ke| (dux doy + Òuy dar Ons Oxy) Dare s(@ — ®') — 
RO OR 0,0 e ay Get ae 
MOLOLEL do on _ 8,,4,,) D apo — -æ') + 

Pe (Oy! O <I Opt Oph O ND PAIE 
(0 ga a) epee 3 bs.) Dara ea a’) + 
+ (6.4 ae aE de da 7 Opy) D is) 
+ (0,5 os 0; dati Or Op, JD st til 
— (0,8 Ot dde — Oye Opy) D naga dj 
— (On Ons Gn adi 3) D sep (@ — x) + 
HOG set Ohh Ongar OO pe) Di ei 
(02040 de 3, DATE ieee 
Fa RUE EE O45) D popyl® Sh 
+ (842905 + Ono bag Ons Ops) Dea 
om (CAT dus On: Ò 1085) DE ia 
ns (0,8 do + 9,5 Oo D983) D peo (€ a!) + 
Te NO pg O 0. OD 


where /% denotes Planck’s constant divided by 2x. 


* 


() We may also denote two space-time points by (15 Yr. 21 4) and (x, Yo, 22 bo) 
(*) B. S. DeWrrr: unpublished manuscript. 
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Consequently the commutators between the gravitational field of electric- 
type components and those between the gravitational field of magnetic-type 
components are given by 


(2,12) [P(@), En(æ)] + [H;;(), Hyx(e')| = 


WO = È 
= A Ke(dad;; + ddr — di;du) VD (x 


a") , 


and the commutators between electric-type components and magnetic-type 
components are 


BS 5 ih |, ONE 
(2.13) [Z:(&), Hy(x')] LT Ke(dinE;im + Ain€iim — dij Exim) oot an Da—a'), 
where 
0 el o? 
(2.14) dix = da = Va, Vi =. 
Om" oak Ox* da' 


In the case of equal times of two field components, these become 


(2.15) [EF is(y), Bely)| = (AY), Hay] = 0 
and 


> Taal NG th 2 ove Ji 
(2.16) [E.;(Y), Hx(Y )| = me Ke(— der FE din Eitm e di; Exim) Ogm d(Y Tr Yi ) a 


Here H,,(7) and H,;(x) denote the gravitational field components of the electric 
type and the magnetic type respectively at the space-time point (+, y, 2, t). 

Our approach to field quantization, in which the first approximation is 
the quantization of pure source-free gravitational fields, enables us to study 
the measurability of the quantum gravitational formalism from the commu- 
tation relations as given above. 

The lack of commutability of the field components which can be regarded 
as the typical character of quantum gravitodynamics is a pure effect of the 
unity of classical general relativity with quantum field theory. It has no con- 
nection with the problems of elementary particles such as singularities due to 
their self-energies because we can neglect the atomic structure of matter. 
The use of the test body as a measuring apparatus is always restricted to 
the extent to which one can handle its behaviour in the given field, as well 
as its reaction to the field on the basis of classical general relativity. 

As we have emphasized earlier, the quantum theoretical field strength has 
to be considered as a space-time average instead of a point function of the 
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classical field strength; in the case of gravitational fields it is of decisive im- 
: a È BARS TA TE ae r 
portance to consider two kinds of field averages over a space-time region. The 
first kind is an average of a magnetic-type components H(x) over a space-time 
. . . ì IT(R) 
region À in the ordinary meaning, denoted by H 


Tie L/2 1022 
1 
(2.17) fiw — |), Jarfavfa 2 Jarre). 
—Ll2 —Ll2 —L|2 T 


and the second kind is a weighted average of an electric-type component Æ(x) 
over a Space-time region À, denoted by 


"i DIL 
(2.18) Em — LT [asaya far (a Gear all Sco iF r.) E(x), 


—L12 —L/2 —L/2 | 


where L is the linear dimensions of the cubital space region V, x, the x co- 
ordinate of the center of mass of a test body covering this region, and 7 the 
time interval of the space-time region A. 


3. — Commutators for field averages. 


In the preceding section we have established the commutators of field 
components, which are regarded as definite space-time functions. As em- 
phasized earlier, we will ignore the atomic structure of the test body and hence 
we shall need the commutator between field averages. However, there are 
two kinds of field averages which play an essential role in the measurement 
of gravitational field strengths; the electric-type components are always taken 
as weighted averages over a space-time region, while the magnetic components 
are taken as ordinary averages over a space-time region. 

By taking weighted averages for the electric- type components and ordi- 
hary averages for the magnetic-type components of the commutators given 
by (2.12) and (2.13) properly over two given space-time regions I and IT ac- 
cording to (2.17) and (2.18), we get the commutation relations between two 
averages of electric-type components 


VihK È ‘ 
3.1 DO FO e ee ian 
| | ; Hire) il RO) 
il 


(dd, + dada — dda) V4 D(x'— x), 
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and the commutation relations between a weighted average of an electric-type 
component and an ordinary average of a magnetic-type component 


| = SK n a ! È 
= anne [22] de'ote) 


I II 


Cer DI Ign 
(3.2) [Bi . elba 
É (= dix E jun =a dir Eium 5 Exim) A DEA È 


and finally the commutation relations between two averages of magnetic-type 


components 
(3.3) [ES Hob] FE Le bre Data d;,dj,— d;;du) Vt D(z'— x) 
ALT: Liat nc i dA ai eda Fo 
I If 
Here 


(3.4) 


and 


(3.5) fate fata e [ar] arfar Jar. 


: is se in Wang Che 


Since we shall always choose center-of-mass co-ordinates for each test body, 
we shall drop hereafter x, and 3 the co-ordinates of the center of mass of 
each test body with respect to the space fixed co-ordinates, in expressions (3.4). 

Since we are interested only in the commutators between the weighted 
average of the Æ,, component and all field averages over the space-time re- 
gion IT, we write down the expression explicitly : 


(3.6) E AA AE) 
and 

[Ba HI = 0, 
(3.7) ; À 

| ca, ay = an Ba Rom, 


where i and j in (3.7) can not take the value 1 at the same time. 
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Here for compactness we set 


ii 9K ae 
x ap — "_ | dx | d4a' w(x, x')(d, VID" — x) 
P Ai LT LA di | n (du ; 
4 I I 
ve ALD = _ IK PA fata E (a, ')(24%, — dd») V4 D™*(w'— x), 
Si D MT Tue, 
Ne : I II 
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te: Ata = IT, Te foe o(#, ©) (dy, dy.) V4 D™*(x'— x), 
it I 
ZII dix | dix’ w(a dis) V4 D™**(a’— x) 
1113 ane | ? æ')(d 11 13) ? 
RO ee SS [are Gun #')(2dy2d43 — dur das) V4 D™*(x'— a), 
Eire LT Ly Tne | 
| if 
and 
Da _ 3K Sw 11 o(æ)(d;;V?) se DE aa). 
ne ere TiLiTue, cot oz 
IT 
= 3K i e re ! DI AD 
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II 
= SK | ; t / 
(LD) ee Re cea 4, 1 ET 
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I II 
= 3K f 0? 
RENEE Bei = 4%! «(x soe DIES 
dato SATTA] * far o(x)(dy,V? rai (Ala), 
| I If 
Th ‘ 
HET SK È OV: 
LU Bam _ —— , 19 49) o(: Val ; — dd ne si (I, A 
bs Bucs DIRT TAT ae | d fa L'0(L)(2d,9 E13 11 £937) Cot Ox, (x 2) 
Ù I Il 
12 
tr from which the remaining expressions A™ and 8! can be obtained by the 
> appropriate interchange of the suffixes in the retarded Green’s function 
a DT — a). 
“ On the other hand, the products of the complementary uncertainties of two 
o weighted field averages follow immediately from eq. (3.6) and (3.7): 
(6.0) AGR ARE 1 |A Amp 
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(: FO AD __ 4 TI) À FO) | DAI BAL 
(3.11) AB? AH® — 0, NECA A (Sai 


where 7 and j in (3.11) can not take the value 1 at the same time. 

Quite similar to the case of the quantum electromagnetic formalism (except 
one has instead of electromagnetic field averages the weighted gravitational 
field averages of the electric type, which in fact do not affect the characteristics 
obtained in the former case), several important results follow from the above 
expressions: The quantities 4°" and B®' change in a continuous way under 
continuous displacements of the boundaries of the regions I and IT as long 
as the dimensions of the region are finite. The differences A°™— 4%” and 
BB vanish without discontinuity when the two regions gradually 
coincide. 

We turn now to the commutators between the average of the H,, com- 
ponent over the space-time region I and all field averages over the space-time 
region II. From (5.2) and (3.3) we can write down the commutators explicitly 
as follows: 


(3.12) (Hay let Cie One) 
and 
(3.13) [HE Et == 0, ARR A pep Dar), 


where i and j in (3.13) can not take the value 1 at the same time. 
Here for compactness we set 


| (ue dix | dix'(d,, V2}? Dt(æ'— x 
Ci AT L date] 4 e (du V?)? ( is 
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iba = Lys pito, Dis (200,03) VED (a — x), 
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from which the ©" and D""% can be obtained by appropriately inter- 
changing the primed and unprimed quantities in the retarded Green’s function 
D'"(x'— x) in the above expressions. 

The corresponding products of uncertainties are 


(Se16) AB} ABS = h| O87 — O82 | 
and 
(3.17) AHO ADS (e GAH ARO Sh pt pe 


where © and j in (3.17) can not take the value 1 at the same time. Just as 


È r 5 È : “are à Fa ul PL II 
in the discussion in the preceding part of this section on the As and Bg, 


the OS and Ds have the same property of continuity under continuous 
displacements of the boundaries; also the differences CTP G™ and De — 
— D™ vanish continuously if the boundaries of the two regions are brought 
gradually to coincide with each other. 

From (3.6) and (3.7) (or (3.12) and (3.13)) it follows that the weighted 
averages of any two electric-type components (or magnetic-type components) 
over the same space-time region commute with each other and therefore can 
always be measured simultaneously with an arbitrary accuracy independent 
of one another. This can be seen from the Symmetry properties of the integrands 
in relations (3.6) and (3.7) (or (3.12) and (3.13)), say A‘ — 4@» (*) and 


Bu» — fav as follows. If the two time regions are identical, 7,= 77, from 
(*) We have denoted any one of the integrands of (3.8) by 40. and similarly for 
Bo», 64.2 and DA», 
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the antisymmetric properties of the differences À4:2— 4,» (or Cr (12,0) 
with respect to t, and ft, we find that the time integrals in the differences 
Ae — AT (or O°” OU) vanish. Hence the weighted field averages (or 
the ordinary field averages) of the same type over different spatial regions com- 
mute provided their associated time intervals coincide. Therefore 


(3.18) DECANE DT) 0s AHO TVA VII 


If the two spatial regions are identical, V.=V,, then owing to the sym- 
_ metry of the Green’s function D(a,—«,) under interchange of the spatial 
points (7, Y,,%) and (%, Y:, 2) together with the antisymmetry of the deri- 
vatives in the differences BU? Bar (or pa» — Pr) with respect to the 
spatial points, it may be inferred that the two field averages of different types 
over the same spatial volume but arbitrary time regions commute. Thus 


(3.19) NESTOR GAZA 


From (3.18) and (3.19) we get the required result that any two field components 
averaged over the same space-time region can be accurately measured simul- 
taneously. 

However, if one follows the corresponding process treated by HEISENBERG (*), 
in which a test body is considered to be a point particle in the limiting case, 
by the same reasoning as in the quantum electromagnetie formalism the dif- 
ferences Bo BUD (or DT DD) suffer discontinuous changes under con- 
tinuous relative displacement of the two regions I and II. Due to this dis- 
continuity the averages of two field components of different types over the 
same space region do not commute. Since we have ignored the atomic structure 
of matter we are not interested in this case, but are interested only in the case 
where the spatial dimensions, the order of magnitude of which we shall denote 
hereafter by L, is much larger than e multiplied by the time dimension 7’. 

The case L<cT is actually of little interest in the domain in which clas- 
sical theory is valid because upon taking a spatial average of the field over 
a region V the peculiarities of the field cancel out. An analogous cancellation 
of the fluctuation phenomena occurs in principle in the quantum domain. 
The statistical fluctuation resulting from the character of the formalism there- 
fore plays a minor role in the case L> ct, but enters in essentially for the 
case L<et, which we shall discuss below. 

The distribution density /(p) of a given graviton with isotopic polarization 
and given momentum hp and energy Æ = hep represented in the form 


DE | 1 
(3.20) as Pm) bip — PI) = Papa), 


(*) W. HEISENBERG: Physical Foundation of Quantum Theory (1930). 
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in terms of the annihilation and creation operators of a graviton a;(p) and 
a,(p) belonging to a definite polarization A of momentum state p and |P) 
belonging to the state with definite graviton composition. 

Although the expectation value of a weighted field average, for example 
E, vanishes, its expectation value of the mean square fluctuation as well 
as its vacuum expectation value of the mean square fluctuation does not 
remain finite but diverges. 

Since we are not interested in the classical domain, we will not discuss 
the problem of fluctuations further, but will turn our attention to the truly 
quantum domain of which the critical field strength w, which is detined as the 
square root of the order of magnitude of expression (3.10), can be found by 
a simple order of magnitude estimation in either case; when L< eT 


DE RAD RE LOT AA 
(3.21) VAKGAE — | SIAT: U, 
and when L> e7, 
TES = hG 
9: / (ID) Pa 
(3.22) VABE® ABD ~ | io È 


In the case where field strengths are larger than « we reach the domain of 
the validity of classical general relativity in which quantum features lose their 
significance. 


4. — Measurement of electric-type field averages. 


The measurement of electric-type gravitational field strengths appears quite 


| different from that of electromagnetic field strengths. It is better to start our 


discussion in a manifestly covariant form and later on go to the linearized 
theory. 

It is more convenient to begin by considering two adjacent particles of 
mass g which describe in space-time world lines 0 and C’, which are not geo- 
desics, in general, due to their mutual interaction. If we signify C and C’ 
respectively by a(t) and «'(t’), t and 7’ being proper times along C and C’ 
respectively, and set up a one-to-one correlation between points of the two 
curves and introduce the vector 7"= a" — x", we get 


(4 1 ) 57" Lu v oO 2 v (e) 
: Sa (U#.,U).cN° + RU UT, 
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where “= dx"/dr, 8/57 denotes absolute differentiation and the dots followed 
by indices, covariant differentiation. Thus far, we have considered the devia- 
tion of co-ordinates between two idealized mass points, the masses of which 
vanish everywhere except in the interior of the two small world tubes C and 0°. 
These tubes are supposed to be concentrated near the world lines so that their 
dimensions are much smaller than the dimensions of 7". 

On the other hand, measurements of gravitational fields require an appa- 
ratus consisting of an elastic medium made up of a continuous distribution 
of mass over an extended volume. In order to treat these macroscopic pro- 
perties of matter, rather than its microscopic behaviour, one may regard 
eq. (4.1) as an equation which describes the deviation of co-ordinates between 
the two centers of mass of two volume elements with mass density 9 which is 
uniform throughout the whole volume of a test body. World lines are thus 
supposed to be the world lines of the centers of mass of these volume elements. 
Thus we regard eq. (4.1) as the dynamical equation of motion of a system of 
continuous matter. When no volume forces act on the test body the dyna- 
mical behaviour of the elastic medium is, in fact, limited by the conditions 


(4.2) (octul’u' 442), = 0, and DEU, =) 


where ¢“” (*) is the internal stress tensor of the continuum which is due to the 
mutual forces exerted between the particles contained in the two adjacent 


volume elements. Eq. (4.1) together with (4.2) completely determines the | 


dynamical behaviour of the system. In the case t“"=0, that is, when inte- 
ractions between the particles contained in any two volume elements are 
neglected, C and ©’ are geodesics. The first term on the right-hand side of 
eq. (4.1) therefore vanishes and the equation reduces to the equation of 
geodesic deviation (1%). By eliminating the first order covariant derivatives 
of uw“ from eq. (4.1) and (4.2) we have (8) 


(4.3) (Of + c-2 yu! Uy) ——=— £ (Si + c-2 UH Uy)[ o- "(Op + wua)t9.:].02 + 


Up till now we have kept the calculation manifestly covariant. We now 
linearize the above equation so that hereafter we shall not need to distinguish 
subscripts and superscripts. Furthermore, if we neglect terms containing 7 
higher than quadratic and if we suppose the velocities of the particles are 


(9) P. G. Beramann: An Introduction to the Theory of Relativity (London 1947), 


chap. VIII. 
Ù A N " / 
(20) J. L. SYNGE and A. SCHILD: Tensor Caleulus (Toronto, 1949). 
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much smaller than the velocity of light, we have 


d?7; 


(4) dt? 


9 j =, il j 
=—¢ E 050 9° — 08 binennl] 


Let us now denote the components of Euler co-ordinates (the space-fixed 
co-ordinates) of the centers of mass of the two volume elements of the test 
body by % and 2’ and their differences by dv’. Then the above equation 
becomes 


7 3 / 


(4.5) 0,— (Chios 6 * ben a) (Ly — %;) ; 
where the dots denote ordinary differentiations with respect to local time. 
In Lagrangian co-ordinates (the co-ordinates moving with particles) where the 
co-ordinate components of elongation are denoted by dx,;, we get 


(4.6) CAC + dex) — &;(0,) = do, + Ô:,; da, 


where dx, is independent of time. By substituting (4.6) in (4.5) and neglecting 
higher order infinitesimal quantities on the right-hand side of (4.5), we get 
the fundamental dynamical equations of motion 


(4.7) ÔË, , — — celo =F On bine ps ; 


from which we shall start our discussion of the measurability of a given gravi- 
tational field by considering limitations imposed by the Uncertainty Principle. 
In the linearized theory we measure the components of the Riemann tensors 
because they are invariant in the gauge transformation as discussed after 
eq. (2.4). Equation (4.7) apparently fits our purpose since it contains the 
electric-type components Y,,. 

In measuring gravitational field strengths one finds it indispensable, in 
order to attribute a definite meaning to the field strength, to assume that the 
mass distribution of a test body is uniform, Considering that a test body has 
atomic structure and therefore consists of many elementary particles, it is 
apparent that the assumption of a uniform mass distribution in the test body 
is already an idealization which consequently imposes limitations on its ap- 
plicability. 

Let us now consider a cubical test body with linear dimensions L of an 
elastic medium which not only has a uniform mass distribution, but also is 
constrained, for the sake of simplicity, rigidly in the y and 2 directions and 
constrained to stretch uniformly in the # direction. #,, is then Just the internal 
stress needed to maintain this constraint. In order to measure an electric-type 
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gravitational field component, for example H,,, averaged over a volume V 
and a time interval 7, one has to locate the test body to cover the spatial 
region Y. For the measurement of the electric-type component it is also of 
fundamental significance to consider a hyperbolic weighted field average over 
the space-time region. Since the field equations are invariant under the gauge 
transformation @,,= @,yt Su En) We May suppose that the center of mass 
of the test body is always at the same position insofar as the measurement 
of electric-type components are concerned. 

By putting dr =e., dy=0 and de=0 into eq. (4.7), one has the equation 
for measuring the non-uniform electric-type component in a single direction x, 
with its origin at the center of mass of the test body, 


(4.8) En CEn- ty a. 


From the boundary condition that the internal stress vanishes on the surface 
of the test body, t,;= 0, and from the condition that the internal forces can 
do no net work, 7.¢., 

LI2 LI2 LI2 
(4.9) far] dy | det one = (}. 

119 — 7/2 —L/2 

we find 

LI2 LI2 ZI2 


(4.10) far] do | donne OI. 


=f/2 —D[2) —L/2 


Taking the weighted average of equation (4.8) over the volume J’, we 
obtain the dynamical equation of motion 


(4.11) Eo = CH, , 
where we have used (4.10) together with the expression 


(4.12) VO Ti far ay fau Pre 


— 112 —LI2 — 112 


and the relation 


i f Ls 
(4.13) Jarfay| accu =. 


—DI2 —LI2 —L/2 


Here e, is by definition the strain, or elongation per unit length in the x di- 


Vi 
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rection, which is evidently of decisive importance in the measurement of elec- 
tric-type components. 

The Lagrangian of the test body system is then defined by the total work 
done on the system, %.e., if | [oda dar dy dz dt: 


M 
(4.14) $= 5 lies, 


where M is the total mass of the test body and é, the time rate of change of 
the strain at some instant ¢. Following the usual procedure for partially dif- 
ferentiating the Lagrangian with respect to the time rate of change of the 
given canonical variable ¢,, we find z,, the momentum canonically conjugate 


tO £,, as 
Wi be 
(4.15) Ris = D JE 


By locating the test body to fill the volume V and by measuring the cano- 
nically conjugated momenta x, and x, of the test body in the x-direction at 
the beginning # and the end #” of the time interval 7, we obtain the equation 
for the balance of the strain momentum of the test body with the weighted 
field average Æ,, over the space-time region 


(4.16) a, 7, = GME, T . 


Here it is assumed, of course, that the time intervals required for the strain 
momentum measurements, whose order of magnitude we shall denote by aa 
can be regarded vanisingly small compared with 7; and that the magnitude 
of e, is much smaller than unity. In other words, compared to the linear di- 
mensions Z of the volume we can neglect the total elongation ÔL of the test 
body system. This elongation results not only from the measurement of the 
strain momentum, but also from the corresponding acceleration of e, imparted 
to it by the gravitational field, which is to be measured during the time 
interval 7. 

Because the strain e, and its canonically conjugate momentum x, of the 
test body system satisfied the quantum commutation relation 


(4.17) [ele 


each measurement of the strain momentum z, carried out with an accuracy 
Ax, is accompanied by a sacrifice in knowledge of its strain to a certain extent 
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Ae, which, according to the Uncertainty Principle, is given by 
(4.18) Ae, An, of . 
By using (4.16) and (4.18) and by neglecting At and ÔL compared to T and 


L respectively, we get the corresponding inaccuracy of the field measurement 
AF,, in the order of magnitude 


(4.19) AB 12h 
raat i” MITE, 


For any value of Ae,, no matter how small, AE,, can be measured with 
arbitrary accuracy if the mass M of the test body is chosen sufficiently large. 
Aside from the mass of the test body, the accuracy of the field measurement 
also depends on the absolute magnitude of the value of #,, itself. This is 
because for given fringes of At and Ae, the value of £,, from (4.16) is affected 
by an uncertainty of the strain and an uncertainty due to the inaccuracy of 
the boundaries of the measuring region, both of which go beyond all limit as 
E,, becomes infinite. Another difficulty arises when W is large; the accurate- 
ly measured gravitational field so far described not only contains the exter- 
nal field that we wish to measure, but also the so-called self-reaction of the 
test body, #.e., the gravitational field produced by the test body itself during 
the measurement-time interval. The self-reaction is very large if M is very 
large. We shall discuss this problem later on. 

Now let us consider the validity of the linearized theory and hence vacuum 
fluctuations in a space-time region of dimensions ZL. Consider the Feynman sum, 


(4.20) > exp 


where the summation is carried out over all possible histories, the order of 
magnitude of the exponential only being given. We see that the contributions 
will add in phase for metric fluctuations smaller than Ag which is given by 
putting the quantity in the exponential equal to unity: 


(4.21) Agira 


Corresponding fluctuations in the Riemann tensors are 


/hG 1 


(4.22) AR = / c3 Ta Py 


Thus, the minimum linear dimensions allowed for the validity of the linearized 
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theory are given by 


hG 
4.23 —,-~1, 
| GLi 
Lin being of the order of 107? cm. 


On the other hand, this is also the limit for which the magnitude of €, is 
much smaller than unity throughout the measuring time interval 7°, because 


beets 
(4.24) AR L? ~ FO 


Di 


is much smaller than unity if Z is much larger than L In the domain 


min * 


where LZ is smaller than ZL there will be a non-linear effect which will go 


min ? 


beyond our discussion. As mentioned in Section 3, we are interested only in 
the truly quantum domain in which we must have 


(4.25) AH: <4 


for L> eT where u=L,,,/\V L’cT’. Moreover, defining a dimensionless factor 
À by AK,;-°3Z,u we have 


(4.26) ; MZ Go x 
220 Vai cae U, 
MG Ae, 


which is an important factor for the estimation of the accuracy of measure- 
ment. The requirement that 4, as well as Ae, be very small as compared to 
unity means that the test body must consist of a very large number of ele- 
mentary particles of mass m, this number being 


A M ID TE he 
4.2 Ne cata | pi 
(agi m el eT | ma’ 


which is very large if Z>cT and 7,«&1 are assumed. 

The factor Gm?/he whose smallness is of the order of magnitude 10-* 
taking the value of the mass of a nucleon, i.e., m = 10-24 g, appearing in the 
last expression is the important constant appearing in the quantum gravita- 
tional formalism corresponding to the fine structure constant in quantum 
electrodynamics. This smallness justifies our formalism very successfully. The 
assumed largeness of N is a necessary condition for the physical realization of 
a uniform mass distribution of the test body over the Space region, and the 
larger the cubical dimensions L of the test body compared with atomic dimen- 
sions, the better it fulfills the physical requirements in principle. 

In what follows, we shall discuss the reaction of the test body on itself 
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during the course of the measurement. Our main problem is to find the most 
physically attainable measuring arrangement so that the reaction of the test 
body can be minimized or compensated for to the largest extent. As we have 
mentioned in Section 1, we are able to treat this reaction classically and not 
as a quantized field owing to the smallness of Gm2/he and as far as it can be 
classically computed, the actual compensation of the reaction by properly 
chosen spring mechanisms and their force constants can really be shown. 

Since we use an extended test body whose dimensions are much larger than 
atomic dimensions, the corresponding difficulties on the measurement of the 
strain momentum, which arise in the electromagnetic formalism as considered 
by LANDAU and PEIERIS, do not arise in the measurement of either electric- 
type or magnetic-type components. 

Furthermore, let us consider a system of two test bodies. In this case not 
only a discussion of the self-reaction on each test body, but also a discussion 
of their mutual influence is necessary. The measurement of a space average 
of the field reaction, H,,,..4, at a point P, of the test body II in volume JV, 
as a result of the strain of the first test body contained in a small volume ele- 
ment dv, of the volume V, at a point P, is exactly equal and opposite to the 
measurement of #,,...4, the field strength produced by the reverse process 
at P,. This-fact implies an essential feature of the formalism and hence when 
the boundaries of the two volumes coincide, the mutual influence of the field 
reaction completely vanishes. Actually in the above discussion of the mutual 
influence of the disturbance, one may take into account this effect between 
two small volume elements of one test body. When the strain increases, there 
will be an increase in the metric which in turn increases an amount of the 
momentum-energy tensor. 

Owing to the fact that in the linearized theory the action is propagated 
with the finite velocity of light, we have to take into account this effect on 
the strain so that we are not permitted to regard the uniform stretching of 
the test body literally, but must regard the test body as a system of consti- 
tuent particles of small dimensions. We must then think that the measurement 
of the total strain momentum x, of the test body is being carried in such a 
way that, to a sufficient approximation, pach constituent particle suffers a 
displacement proportional to the magnitude of the x co-ordinate of its loca- 
tion from the center of mass during the strain momentum measurement. 
Therefore the strain, or elongation per unit length, remains the same on the 
average. The above requirement can actually be fulfilled so far as we can 
neglect the atomic structure of the test body. This, in fact, suggests that 
the strain momentum can be measured on a classical basis, for example, by 
a study of the Doppler effect in the reflection of photons, provided a bundle 
of photons sent to the test body is large enough to determine accurately the 
interaction between the test body and the radiation. 
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Since the space-time weighted average of the #,, component can be written 
as 12 SY 2,(p,— p.)/ML?eT, where the summation is extended over all com- 
ponent particles 7, we know that the strain momentum 7, = Nepi > as, 
where m, is the mass of a component particle. 

Making use of (4.14) the Hamiltonian of the test body system is given by 

6 È 


DIS H == J & : 
1520) TEE 


from which it follows that 

(4.29) NE == ENT. 

From this relation together with the uncertainty relation between the time 
interval At and the uncertainty AW in the knowledge of the energy, exchanged 


between the test body and the bundle of photons during the collision process 
in the measurement, %.e., 


(4.30) AH Ai~ hi; 
we get 
(4.31) Ax, |é, —& Ath. 


The accuracy of the strain may be known within an uncertainty 
29 7 . II at 
(4.32) Ag, =|, —é,|At, 


which is in complete agreement with the uncertainty relation (4.17), regard- 
Ù PAT Ahi x a 

less of whether |é,—é,| can be known accurately for a sufficiently heavy 

test body. Recalling e, to be the elongation per unit length, we have 


(4.33) A asa (7 | 


which together with the relation |e,|<1 gives rise to the relation 
4; c 
(4 104) INS << TA ; 


which limits the absolute accuracy to Az, for a given upper limit At in a 
measurement of the strain momentum. 

We can neglect the special-relativistic effect between any pair of adjacent 
constituent particles of the test body by using sufficiently heavy test bodies so 
that the relativ e velocities of any pair of adjacent constituent particles are 
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always much smaller than the velocity of light. Therefore we are also able to 
regard the uncertainty of the total elongation Le, of the test body, resulting 
from the strain momentum measurement, as being very small compared with 
¢ At, which aslo has to be chosen very small. Consequently, we can arrange 
the measurement so that 


(4.35) [Ae, | <& 1 


always holds. A test body must be thought to consist of a large number of 
component particles which interact with a bundle of photons. 

In order to measure the strain we use the Doppler effect by optical means 
quite similar to the considerations of BoHR and ROSENFELD. We imagine 
each component particle to be equipped with a small mirror which is per- 
pendicular to the x direction. In addition to the mirrors attached to each 
component particle, we also have properly situated mirrors to make the light 
path from the light source to each component particle the same so that by 
sending a bundle of photons to each component particle, we may know the 
momentum of each component particle, except for an unavoidable uncertainty, 
through the spectral analysis of the reflected photons. We can now imagine 
the whole apparatus to be so designed that the number of reflected photons 
from the mirror of a component particle is proportional to the magnitude of 
the x co-ordinate of the component particle referred to the center of mass, 
which is the origin of the co-ordinates of the test body. The other half of the 
test body is imagined to have equipment similar to the half described above 
together with the same light source. | 

We send a bundle of photons of duration Af so that each component par- 
ticle will receive a momentum which is proportional to the magnitude of the x co- 
ordinate from the center of mass of each test body. At the same time different 
component particles suffer different accelerations. Attention must be paid to the 
fact that when we send a bundle of photons of a definite frequency from the light 
source, the frequency will change when it collides with a component particle 
of the test body because of the gravitational red shift. This correction must 
be taken into account for a detailed analysis of the measurement, but we will 
not be concerned with this effect in the following analysis and will simply as- 
sume the condition that the stress uncertainty of the radiation bundle is small 


as compared to that of the test body, 7.¢., given by the condition VAG/e <1, 
which is also the condition for the validity of the linearized theory. 


> 


Let us use the same notation as reference (2) and follow the same method. 
We get 


hh In, oly 
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451) Vretôe T Ur.atta — mc ’ 
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where ihe and Ve» denote velocities before and after collision of a component 
particle at x, m, is its mass and the summation extends over the n, photons 
reflected by the component particle. We have denoted by »’ and »" the fre- 
quencies of the photons before and after collision with the component particle. 
The mean frequency v of the bundle is assumed very large compared to (At)-} 
and »’—y". 

Subtracting eq. (4.36) and (4.37) and dividing by the stationary length dx 
between the two adjacent component particles, we obtain 
(Vata — Urn) — (Vre+te—Vr0) 2 const hry 


(4.38) dai = m,¢ 


from which by summing over ail component particles and dividing by the total 
stationary length we have 


(4.39) eb è 


which means that the time rate of change of the strain is on the average the 
same between any two adjacent component particles. The impact of the 
photons therefore produces an uncertainty of the strain Ae,, which is uncon- 
trollable but the same between any two adjacent component particles on the 
average. 

By virtue of (4.30), (4.31) and (4.39), we have for the energy of the com- 
ponent particle 7 


(4.40) EM [ D w=") >» (y! —v)] 


and hence 


(4.41) Ar, Ne, o AD h [ > (v= »") — Do y") |} , 


1G Nr ,x+0 NT a 


Because the summation > extends over all component particles and the dif- 


tT 


ference of the summation, Y — Y, over all adjacent pairs of particles, the 
Nr eton "To 


part of relation (4.41) inside the curly brackets becomes (9) 


(4.42) Di'-v)-Yho —v) 


LATE "10 


() These cancellations are valid only if one can ignore the retardation of forces. 
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Or 


in which the last summation can be ignored, since we can always make an 
arrangement with 7,70, te, no photons being sent to a component par- 
ticle located at the center of mass of the test body. Then we can write 


(4.43) Anz Ae, => At A [> DI Rly’ — v')] 


E]. 


0) 


which is identical with the uncertainty relations (4.17), where Ae, is assumed 
as before to be vanishingly small and uniform over the whole test body. 

Despite the reflected radiation being measurable exactly, the incident beam 
can only be known to an accu racy complementary to At. Even though, in 
general, we cannot obtain information on the strain momentum of any pair 
of adjacent component particles, it is possible through the spectral analysis 
of the reflected photons to measure a statistical distribution of the strain mo- 
menta which will show a hyperbolic distribution. This fact is closely related 
to the hyperbolic weight we have to take in the averaging of electric-type com- 
ponents. We must further emphasize that the test body system suffers a uni- 
form stretching under the strain momentum measurement and that this fact 
is of importance not only in insuring the cancellation of the self reaction, but 
also in enabling us to regard all test bodies as not suffering strains outside 
the small time intervals for the strain momentum measurement. 

Because the change in strain velocities transferred to each pair of com- 
ponent particles by the initial impact is known, one can compensate for this 
effect by giving the test body system another oppositely directed impulse. 
This can be done after the initial strain momentum measurement, still during 
At, without a sacrifice of knowledge concerning the strain momentum of the 
test body; so that within an uncertainty LAe, the test body covers the space 
region, the field strength of which we want to measure. After the time inter- 
val T following the final measurement, the same process has to be carried out 
again. During these collision processes it is impossible to know the time inter- 
val and the strain more accurately than Aft and Ae, respectively. 

It is always necessary to know the exact position of every test body betore 
and after each measurement. In order to know exactly the position of each 
test body before and after its use as a measuring apparatus, outside the meas- 
uring time interval 7, one binds each component particle of the test body to 
each component particle of a rigid framework body which occupies nearly 
the same space region and on which the detecting apparatus of radiation 
bundles may be mounted. The mass of the framework body is supposed to 
be much larger than that of the test body so that the oscillation period of the 
framework body is much longer than the measuring time interval 7. The 
mass of the framework body, however, cannot be too large, otherwise it would 
violate the condition of the linearized theory, è.e., M,@/Le?*<1, where the 
mass of the framework body is M,. At the beginning of the time interval 
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T the binding has to be loosened so that the sum of the strain momentum 
of each component particle can be measured, while, by an immediate counter 
impulse, one must be able to bring the test body to rest with a degree of 
accuracy LAe, allowed by the uncertainty relation. At the end of 7, the rigid 
binding is imposed again on the test body so that it occupies exactly the 
same original position. 

In the measurement of electric-type field averages over two partially over- 
lapping space-time regions, one requires idealized test bodies which can be 
relatively displaced through each other without mechanical interaction. We 
are able at least in principle to imagine that the component particles of the 
test bodies are bound to the corresponding constituent particles of the frame- 
work bodies as in the preceding discussion. Only in so far as we can neglect 
the atomic structure of the test body can we justify these idealized assumptions, 
as well as compensating spring mechanisms to be introduced later on. 

The above discussion has shown that it is possible to distinguish the limi- 
tations imposed by the manipulation of the test body, which depend on Planck’s 
constant, from the limitations on the structure of the test body due to the 
atomic structure of matter. In spite of the fact that the constants appearing 
in our formalism are sufficient to fix the space-time dimensions, it is required 
for the validity of the linearized theory that one take the average values of 
the field strength over a space-time region. 

The weight one has to take in the case of electric-type field averages is 
closely related to the fact that in the radiation of gravitational waves one 
can always eliminate the dipole radiation by suitably choosing co-ordinates, 
leaving the quadrupole radiation as the first radiative term. 


5. — Measurement of magnetic-type field averages. 


The measurement of the magnetic-type gravitational field strength is quite 
analogous to that of the electric field strength. Consider a spinning test par- 
ticle of PAPAPETROU (11) whose equations of motion are derived under as- 
sumptions on the form of the energy-momentum tensor in a small world tube. 
In the limit of the tube becoming a line, these equations of motion are given by 


dò | 5 al 

~ (6 

( ay RESTA) Yup sal 

(5.1) (ru + Uy 5a SP ae Re UE ee 0 
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(7) A. PAPAPETROU: Proc. Roy. Soc. London, A 209, 248 (1951): A. PAPAPETROU 
and E. CORINALDESI: Proc. Roy. Soc. London, A 209, 259 (1951). 
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where S“” is a skew-symmetric tensor expressing the internal angular mo- 
mentum of the particle. If in addition to these equations we choose the con- 
dition (#7) S“’w,=0, which is a generalization of the conservation of momen- 
tum, we can write (5.1) and (5.2) in the form 


@ Sul Du 
(2.3) Le ÔT Di DEE E + Roo u” S27 = OF 
and 
è dm 
“«(5.4) ve == ) 7 
di 


which expresses the conservation of mass. From eq. (5.3), which is the funda- 
mental equation for a spinning test particle, we shall start our discussion on 
the measurement of magnetic-type components. If the velocity of the spinning 
test particle is much smaller than the velocity of light and if we consider the 
sase in which S?u,/dt? = 0, by linearizing the theory we have 


ane si Cc 
(ay) ML, = —> Fase Sane 


In order to measure the magnetic-type components we have to consider a 
test body consisting of many spinning particles. The cubical dimensions of 
the test body must be chosen much larger than the atomic dimensions so that 
we can assume that the test body has a uniform mass density m. Moreover, 
we assume that the test body has a spin density with the components o, = 
— §,,40, c,=0 and o,=0 distributed over the whole volume so that inside 
the test body mass currents cancel out one another leaving mass currents only 
on its surfaces. Then we may write eq. (5.5) as 


(5.6) Me, = — Cli joss Ox - 


For the measurement of a magnetic-type field strength, for example H,,, 


we need only the equation 

(0.7) Ni CH Ow 

which is exactly the same type of equation as that describing the motion of 
a charged particle of mass m and charge e in an electric field of strength L,, 


if we replace the magnetic-type component H,, multiplied by the spin density 0, 
by H,, and replace the velocity of light by the charge e. Taking the space-time 


(12) F. PIRANI: Acta Phys. Polon., 15,389 (1956). 
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average of eq. (5.7), one has 
"A 7 I oO 1 
(5.8) Pi Dy = CH aral « 


We can proceed with the whole discussion on the measurement of the magnetic- 
type gravitational field in a manner parallel to that of the electric field average 
according to reference (2), but we will drop the topic at this point. 


6. — Electric-type field reactions of the test body. 


In what follows we shall study the field reactions of the test body which 
accompany the measurement of the gravitational field strengths and show in 
fact that these reactions, which are proportional to the strain, can be auto- 
matically compensated for by adjusting the spring constants appropriately, 
just as Bonr and ROSENFELD did in matching the reaction of the electro- 
magnetic field by a mechanical spring. We intend to study the reactions of 
the test body produced only by a single field measurement in this section. 
To follow the foregoing procedure we shall assume that every test body has 
a uniform mass distribution and suffers a- uniform stretching during the meas- 
urement of the strain momentum. 

Let us now consider two space-time regions I and II having volumes V, 
and V,, and time intervals 7, and 7, so located that light signals emitted 
from some points of I can reach IT but signals from II cannot reach I. We 
imagine that the volume VW, is covered by the first test body and V,, by the 
second test body. The measurement of the weighted averages of H,, over the 
region I produces the gravitational field at a point (%, y., 2, #) of the re- 
gion II, by whose effect we can not measure the field strength in the region II 
accurately. To calculate thisi nduced gravitational fiel dat the point (72,Y2,22,t2), 
we may be able to follow a process corresponding to the calculation done by 
BoHR and ROSENFELD by finding the metric induced by the energy-momentum 
stress tensor and then the induced gravitational field, ete., but instead we 
shall follow a direct method to compute the induced gravitational field by 
using the solution (2.7). 

It is possible to assume that the center of mass of each test body is always 
at rest in so far as the measurement of electric-type components is concerned 
so that we can use the center of mass co-ordinates of each test body. At the 
beginning of 7, in the time interval from # to #+At,, the first test body 
suffers an increment of the strain Ae,; it remains at a constant strain Aéy 
during the time interval 7,, i.e., t;—t,. At the end of T,, in the time inter- 
ral from # to bo At, it recovers its original shape. We shall assume that At, 
is much smaller than T,=t—t and that Ae, is much smaller than unity, as 
well as smaller than c At/L. 


Consider now the energy-momentum tensor of the first test body system 


18 
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accompanied by the strain e,, instead of the strain Ae,, which is generated by 
a uniform stretching of the test body in the x direction and maintaining it 
unchanged in the y and ¢ directions. If we denote by L the linear dimensions 
of the cubical volume, then according to the mass conservation law it follows 
that the increase of the energy density Ty, is — oe?e, throughout the vo- 
lume V, and is oc?Le, Ô(æ—(L/2)}/2 and oe>Le,d(x+(L/2))/2 on the surfaces 
æ— L/2 and æ=— L/2 respectively during the time interval T,. Summing 
up these quantities, we have for the total energy density 


oe 5) 3 +) | 5 (x 2) 7 5 92 3) 


/ sé L 7 GE 
-O(t — th) 6(t! —1)6(y a 5)0(5—u)9(z+ ess | 


which can be written in the form 


(6.1) Ty = 00%, 


OPA 
(6.2) pO ae 2 
Ca? 
with 
oc L\ (L L\ ,[L VIE 
PSE (E) + Dot) 
TD dB 1 11 
ale + 5) jou mou —n, 


where e, is assumed to be a constant and 0(x) is a step function with the 
property that 0(x)=1 if «<0 and is zero if x>0. Since the momentum 
density 7’, is proportional to the magnitude of the + co-ordinate and is gen- 
erated only at t=f, in the x direction and at {= 1, in the opposite direction, 


we may write 


(6.3) Ty = — oexes[Ò(t — tr) — d( ole +: 3)! (5-2) 


L L 
UP —5)0(5—y\ o(: + i ; 


| GP 
(6.4) To RIGA 


a: 
OOO 


which also can be written as 


The momentum flux density 7,,, however, cannot be obtained from physical 
considerations directly, since we do not know the internal stress tensor ¢,; of 
the test body. Nevertheless, we can derive it from the conservation law 
OT ,,/ Ox, — (0), because all other components are zero. Hence 

OAR 
6.5 IE ; 
( ) Lil dae 


vO) 
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The energy-momentum tensor of the first test body induced by the uncer- 


tainty of its strain Ae, is thus given by 


A er 22 D 

oP, op, _  ®P, 
AG gp == INGE == NS VAN 2 
(6.6) AT vo Ante? 10 dx 0 x 


with the generalized polarization of the electric type 


(Bo dl Ae,(x) JE Le 5) f) E = 2) 0 ly = 5) 0 È _ ): 
TE +5) (3-00 moi ve 


Replacing the Green’s function D(x,— x,) by the retarded Green’s function 
D'*(x,— x;), the uncertainties of field components at (2, Y2, &, %) due to the 
increment Ae, or the first test body in the region I is 


K 
(6.7) ARwol@)= — 3 | die, D™ (07, 2 (AT, —4 Sig AD) ug — (AT yo — $y AT) pre — 


— (AT $ 0p: AT) yp + (AT o — Os NT) oe] - 


For the £,, component we have by means of (2.10) and (6.6) 
2 Q al K \ 4 ret 
(6.8) AE,(0,)= 8 d'a, D"*(%, — PENE + Pena + Peo000 + Pe1100 — AP, 10 | . 


Integrating by parts on the right-hand side of eq. (6.8) (*) and neglecting all 
terms except those which have fourth order differentiation of the function 
D'(æ, — x), we have 


Ko : ot di 24 | 
6.9 AF;,(@) = — Az, | dir, o(x., x, Lo uu —— ERA 
(6.9) (42) 8 1O(X, Wy) dai e? dn OF ae et (15 ia) 


By using eq. (2.8) and by the same reasoning as in the foregoing discussion, 
the terms involving the 6-function on the right-hand side vanish. Therefore 
we finally get 


KM 


(6.10) AF, (%) = SI? 


Ae, | dir 0(2,, 2) (dy, VID", — è) - 


I 


(*) See Appendix I. 
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Taking the weighted average over the space-time region II and signifying by 
indices (I, IT) the disturbance in the region II resulting from the measurement, 
in the region I, we obtain 


= SK M Ae, 
Valet / EURE Da 1p 4 la x x! 2\2 7)ret L'— % : 
(6.11) AE qe | rf a œ(æ, æ')(d, V2}? D'(x 1) 


il 


By procedures similar to those we have used in deriving AES. we find 
for the other uncertainties of electric-type components 


M 
ATE e = (Ant, Ate da (i a')(2di, — dy, dp) V4 D™*(a'— a) , 
3K M Ae, 2 
AES ID ~ 4312, T ac “fe dx! (x, x’) (2dis GE di, 433) Ve D(x a x) ; 


3K M Ae, 


ATEI av] di Bh ot È (dida) V4 De x) À 


(6.12) NES IT) 


ALS ID —_ = LI. de op! Ole , (dis) VE D (r— r) , 
4 IT 


3K M Ae, ; ho 
ABs; ti a = ammo Li d'o' w(a, æ')(24d35 dis — dy, dos) V4 D™**(x'— x). 


I 


For the uncertainties of magnetic-type components, we have similarly 


a (| 

AHS — 2 || æ' «(a#) (dy, V?) 7 ta D*(a'—a), 

AHL™ — nino fata fa oc w(x)(d;3V?) A DOG), 
(6.13) AH; Genta Slate fa: A æ'o(æ)(d;, V?) > D°'(a'— x), 

AH — sai] "| dix’ o(x)(d,, V?) - Fa D™*(a'— a) , 


II 


KM Ae, 0° re » 
AH ID — sara) dal 4x" 0) ( x) )(2dyE15i— di 1£231) V? 0a On. D ‘(a a). 


192 HSIN YANG YEH 


On the other hand, from the commutation relations (3.1) we get, for the 
product of the uncertainties of two weighted averages of electric-type com- 
ponents, 


Nago es fate auront, oldu tee Dm ea), 
ui LT: Tne | 
I 


il 


= = VIAL } y È i IN fe ret (pl 
ABV AHS BEDE. | dis) dix! w(x, x’) (24, — dd») V4 D™*(a'—®) 5 
dh K | 
ESAES & 14a | d4a’ w(a, a’)(2d$, — dy, dss) V4 D™*(x’— x), 
AEG AES LIT Li Tue C | 2" D ( ')(2di3 11 das) (0° ) 
(6.14) À 
AK | | 
(1) DE 1: 14x' 1 d d x Vi Dret Dr 
NED NES zan) d'a | dix (x, a) (dys dae) ( Nis 
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DEAN DS - [dix | dix’! a Ad) NEEDS 
AE AE LIT LR Tue, Cc | dir’ (a, £')(dy, d,s) (Lx), 
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= = dh K È 
(1) Gide SUA ; 74 Dretlg'—w). 
AESAES & Tre) dur fato (a, 0')(24,3d,5— d das) V4 D"*(e'— x) 
i ir 


Finally from (3.2) we get for the product of uncertainties of the Hy, average 
and the averages of magnetic-type components 


lABYARU — 0, 
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Because in this section we have assumed that no light signals from the 
region II can reach the region I, no disturbances resulting from the measure- 
ment in the region II, such as the one we have computed, can influence the 
measurement in the region I. Hence by comparing eqs. (6.11) and (6.12) with 
(6.14) and (6.13) with (6.15), we find that AZ® must take the value 
12h/M,c? TL; Ae (*). This result does agree completely with the uncertainty 
of the field strength £% given by (4.19), derived from the dynamical equation 
of motion (4.11) by taking into account the Uncertainty Principle (4.18). 

The expressions (6.11), (6.12) and (6.13) relate the disturbances in the 
region II resulting from the measurement in the region I of the #,, weighted 
average to the uncertainty in the strain Ae. Evidently from (6.11), (6.12) 
and (6.13) it follows that 


‘20 ARG” = = ML Ae? AGP , 
and 

AHI =O | 
(6.17) + 


2 c? D 
LID 2 (MEI 
AVES == 12 M LT Aes Pi; 
hod 


where i and j in (6.17) can not take the value 1 at the same time. Recalling 
the continuous natures of the A’s and B’s under the continuous changes of 
the regions I and II, we see that the uncertainties of field strengths in the 
region II are simply proportional to the uncertainty of the strain Ae” of 
the first test body. These properties are of decisive importance for the pos- 
sibility of cancelling the uncontrollable field reactions of the test body by 
means of suitably chosen force constants for springs joining each test body 
to its own framework body. 

In the foregoing discussion we have considered only the effect of the meas- 
urement of the Æ,, weighted field average. In the measurement of £,, and Éy, 
we can properly choose a test body which is uniformly stretched not only in 
the x direction, but also in the y and 2 directions. On the other hand, the 
weighted averages of non-diagonal components of the electric type, for example 
E,., cannot be measured purely by themselves, but are always accompanied 
by, and mixed with, the measurement of the diagonal components of the electric 
type, for example H,, and #,,. This fact can be understood by taking the 
polarization effect into account, in which H,, can be produced by a 45° 
rotation of either the H,, or 7», component. Later on we shall return to this 


(*) Hereafter we will use the indices I and II to signify explicitly the quantities 


in the regions I and II. 
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problem by rotating three principal axes of one test body to arbitrary direc- 
tions with respect to those of another test body so that their principal axes 
are no longer parallel to one another. 


7. — Magnetic-type field reactions of the test body. 


We have obtained the dynamical equation of motion of a spinning test 
particle in Section 5. We shall consider now a test body of cubical dimen- 
sions L which is made up of many uncharged spinning elementary particles 
and has a uniform spin density 0,70, o,—0 and o,=0. Therefore the mass 
current inside the volume is completely cancelled out, leaving surface currents 
at y= + L/2 and z=+L/2. The components of the energy-momentum ten- 
sor of the steady state of the system are thus given by 


(7.1) 


1 L L L L 
est (e+5)9(3-2)p(1+ 2)-2(5)]e (+35): 


If one shifts the test body in the + direction by an amount Ax at time # 
and shifts it back at time # , the increases of energy flux AT, and AT, during 
the whole time interval are 1 


A 7 Ti 
AT 29 = Ot — At ata) JE | 4 (5 DÌ 
L IE L 
bo 
(7.2) | 
Ax da 
Ne LE Ot i) > (a 5) (5 o) > 


and all other components vanish, where Dates 
Because during the whole time interval there is an increase of energy flux, 
there will be à momentum flux at the beginning and end of the time inter- 
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val 7. The components of this flux density are 


C70, L L L L 
Ae 7 A06 (2 + Ce) Aly am Cv) 
L L ! ) I 
fs 3) (5 :) [ d(t — ti) — 6(& —t)] , 
ATs =—"F*Ard(e+5)6(5—2); 


Lo | 5) JE ») a( + 5)9(5—-| Pi: 


From (7.2) and (7.3) one sees at once that the induced stress components can 
be written in the form 


CP. 

pato ce 

0? m 

RE oa 

(7.4) io 
pia 

OP, 
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with the generalized polarization of the magnetic type 


AIM Met DE 


CAS 


(7.5) P, 


It is easy to verify that the induced stress components as defined above satisfy 
the conservation law ÔT,/0æ,— 0. Hereafter we shall regard AT, as the un- 
certainty of the stress tensor produced by the uncertainty of position Ax. 

The uncertainty of the H,, components at a space-time point (#2, Yo, 2, ta) 
due to the uncertainty in position Ax of the first test body is, by a calcu- 
lation similar to the one used in computing the uncertainty of the electric-type 
component in the previous section, 


~ Kx OK SE 
7.6 AH (0%) = sii dr (A V2) D'(%: — Li), 
(7.6) LL 
i 
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where 2 — o,Lî is the total spin angular momentum of the test body covering 
the spatial region I. By taking the ordinary average over the space-time 
region II covered by the second test body, we see that 


HD Ax® 
TT AHL™ — = | 4 ee da V2?) D™* (a'— x 
(7.7) tL Ts | Me 


Similarly we find for the other uncertainties of magnetic-type components 
the expressions 
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For the uncertainty in the weighted averages of electric-type components 
over the space-time region II, we find similarly 


AB T= 0 
AES — ees : fa ‘| d'a' w(a')d,,V? - = > D'*(x'— x), 
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On the other hand, by taking into consideration only the retarded effect 


from the region I to the region IT in (3.16) and (3.17) and by the Uncertainty 
Principle, we have 
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AH®AHS® RI Ma 4 Appl 2\2 7)ret '_ x 
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Comparing (7.7) and (7.8) with (7.10) and (7.9) with (7.11), we find that 
AH must take the value 7/3 AxeT,. On the other hand, the equation 
of momentum conservation (5.8) gives the same result upon taking the 
Uncertainty Principle Ap Av~h into account. 

Expressions (7.7), (7.8) and (7.9) can also be written as 


ry (LIL Ty I) (1,11) 
(7.12) ND ) — Co AN DC 
and 
nn) NO TI (II) OSE D AO 7: Fam 
(7.13) Ari 0 Abi GL ADR) E 


where ? and j in (7.13) can not take the value 1 at the same time. This 
shows that the disturbances induced in the region II as a consequence of the 
measurement of À,, in the region I are proportional to the displacement of the 
position of the first test body. This result plays an important role in the 
compensation mechanism as was emphasized in reference (2). 


From the dynamical equation of motion (5.7) we get for the uncertainty 
of Hy, 


Ap. h 
oT cd, At” 


(7.14) NES 


For any Ax, no matter how small, AH,, can be measured to an arbitrary ac- 
curacy by choosing the total spin 2, sufficiently large. We are interested only 
in the quantum domain, in which 


Lai 


(7.15) AT N 
eci 


with the critical field strength « and the minimum length defined by (4.25) 
and (4.23) respectively. Furthermore, we define a dimensionless factor Ai Dai i 


(7.16) AH ~ Ant , 
which on being substituted into (7.14) leads to the relation 
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The fact that Z,,<1 and Ax< L means that the test body must consist of 


a large number of elementary particles of spin #/2. The number of elementary 


“all 


rei e A! tolé PET Ge ae SA * e PR A 0 . La Le di - L | 
AN CEE F RO Re Le PL Lo RME eo ira AURAS LE ì Goo 
) \ Do ASTE 4 tee « strana ’ 4 
mei NE L Ors | ; ‘ 
LES + | i mes i 
ae . : PRET PA 


MEASURABILITY OF GRAVITATIONAL FIELD STRENGTHS 139 
particles contained in the test body system is given by 


9 3) a Te 
(7.18) fly ele rai 
h A NEILE eT! 


which is very large for 1,<1, L>Ax and L>TL,,,,. 


8. — Measurement of a single field average. 


The measurability of electric-typé gravitational field averages has its basis 
in the fundamental equation of strain-momentum conservation (4.16), which 
relates the strain momentum of a test body to the field strength. This field 
strength actually includes all fields in this region: a given gravitational field in 
the region, a field induced by the self-reaction of the test body covering this 
region, and induced fields resulting from the measurement of fields in other 
regions which may include a part of this region. 

As a preliminary for a more general discussion of the measurement including 
two field averages, we shall restrict our discussion simply to the self-reaction 
of a single field average and then go to that of two field averages. Aside from 
effects influencing the accurate measurement of the electric-type field strength 
at the beginning and end of the measuring time interval (as discussed in 
Section 4) the strain acceleration of the test body produced by the field itself 
may cause the test body not to occupy the region where we want to measure 
the field strength during the measuring time interval. This section affects 
the accuracy of the measurement. Nevertheless, this effect may be ignored 
by choosing the mass of the test body large enough. 

Referring the eq. (6.16) and making the region I and II coincide, we could 
separate a given weighted field average E®, which would be in the region if 
no field measurement were made, and Hf”, the self-reaction of the test body 
resulting from the measurement, in the equation of strain-momentum con- 


servation 


12 


(8.1) nd nd = MATE + E). 


bo 


According to the discussion above and that in Section 4, in order to attain 
an optimum accuracy of the measurement one has to choose the mass suffi- 
ciently large. On the other hand, the heavier the test body is chosen, the 
larger will be the uncontrollable effect of the self-reaction of the test body, 
AE, To reconcile these two opposing requirements, one has to adjust M,, 
L, and the uncertainty of E% to give the minimum uncertainty of AE. 
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Since according to (8.1) 


12 Ar? 


DR alles AES? 
CM LT, TE 


(8.2) AE® ~ 
it follows that, by referring to relation (6.16) and to the uncertainty relation, 


12h e? MIT, Ae® ALD 
CM LT, Ae? LT rit ifs | 1111 | - 


(8.3) AËED = 


From (8.3) the minimum value can be obtained by 


ra Ana Vh AR |; 
which becomes the critical value « in the case Ly > cT, and which may be 
regarded as an absolute limit on the measurability of gravitational field strengths. 
We arrive at the conclusion that the quantum gravitational formalism is 
applicable only in the classical limiting case, but not in the truly quantum 
domain. This conclusion can not be justified unless one arranges things so 
that one can compensate for ies the uncontrollable self-reaction of the test body, 
of which the uncertainty is proportional to the unpredictable strain Ae®; one 
can retain only the field £®%, which may be compensated for by the measure- 
ment except for the inevitable fluctuations. Since we have chosen our test 
body to be an elastic medium, of which the elastic tension per unit area of the 
system in the « direction is equal to — Ye, we can attain these conditions 
by choosing an elastic constant (*) to be 


e? 2 n 
(8.5) Miao & mL) Tr Aïn ; 


as long as the mass density of the test body is large enough and its elastic 
vibration period is longer than 7,: Hence the strain e, is much larger than 
the strain resulting from its oscillation. In the choice of the elastic constant Y 
above, we note that Y is not only dependent on M, and 7,, but also on J,. 
This is a typical result in the measurement of the electric-type field component 
and is caused by taking the weighted average of the electric-type field com- 
ponent. 

Actually the test body regarded as an elastic medium has to be thought 
to consist of many component particles which are connected to component 
particles of a framework body by springs of suitable strength as we have 
discussed in Section 4. Therefore we may regard the oscillation period of the 
framework body to be much longer than the measuring time interval 7,. In 


() The constant Y is actually Young’s modulus in the case of the measurement of En: 
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order to take retardation effects into account, these springs must be suffi- 
ciently short. In the case of a sufficiently heavy test body e” is very small 
and hence these retardation effects can be neglected. 

In the measurement of magnetic-type components, we need also a spring 
mechanism to compensate its reaction which is proportional to the displace- 
ment of the rigid test body as in the measurement of the electromagnetic field 
strengths. 

The spring meachanism for the compensation of self-reactions both in 
magnetic-type components and electric-type components has no direct relation 
to quantum concepts but can, in the limit, be only classically describable. 
We therefore arrive at the conclusion that only the classical description of the 
field action of the test body gives rise to the interpretability of a single field 
measurement. 

Aside from a single field measurement according to equation (4.16), one 
may be able to set up a test body which measures only the weighted average 
of the #,, component. It is in principle impossible to measure the ordinary 
field averages of electric-type components; while the measurement of the 
weighted averages of magnetic-type components has no physical significance, 
for only the ordinary averages of magnetic-type components have physical 
reality according to eq. (5.8). 


9. — Measurability of two averages of the same field component. 


Following the discussion in the preceding section, we shall study in this 
section the more general case in which two test bodies covering the regions I 
and II are so situated that light-signals from at least some points of region I 
can reach region II and vice versa. It is desirable to start from the study of 
the measurability of two field strengths of the same component, which we 
shall do in the rest of this section. 

By considering the weighted field average E,, and by neglecting the limi- 
tation on the classical describability of the field actions of the test bodies 
the equations of strain-momentum conservation of two test bodies are as 
follows: 


| 7 C 72 ga (LD (II,I) 
70 — n= — MITE + HY? + En), 


(9.1) 
TU) vai TR < Malin (Ei EE + per un) Aie Eo ele 


where the uncertainty of #%" is given by (6.16) and Ht”, AY)? by analogous 


expressions. 


“1 
on 
i 
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According to the discussion in the last section, the self-reaction of each 
test body can be cancelled by appropriately adjusting spring constants, which 
are given by expressions analogous to (8.5). By omitting ES” and ES” from 
(9.1) and referring to (6.16) and (4.18), we have the relations for the uncer- 
tainties of the field measurements 


= 12h e? È = 
ARS Sap res + jg MnliluAe | A |, 
IHTUTLI ® I 
(9.2) 
= 12h CA X = i 
AE M I? T Agi 12 IMPERI Ne |A et ls 
È Iaia x Sa 


ORI II 
where we have used the fact that the strains e i 


each other and are known with uncertainties Ae” and Ae. If we simply 
adjust properly M,L?7,Ae® and MI} Ty Ae.” to fit the minimum value of 
the product of AZ® and AH? without further detailed consideration of the 
correlation between subregions of the test bodies, in general, we would arrive 
at an incorrect result (!*) to get the sum of the magnitude of AGI? and AM? 


) 


and e‘ are independent of 


‘in the expression 


(9.3) AE® AB =, flares | Bs PASS |] | 


1111 a is i is 


This result is in contradiction with expression (3.6). Expression (9.3) agrees 


with (3.6) only in the case in which one of the quantities AT? or AMP 


1111 1111 
vanishes and hence one of the two regions is completely unaffected by the 


measuring disturbances of the field in the other. Except for this case, they 
do not agree because (3.6) vanishes even if 445 and A individually differ 
from zero in several important cases such as the case that, when 7, and T, 
coincide, HY and E° can be measured correctly but (9.3) does not vanish. 
The agreement between the quantum gravitational formalism and the meas- 
urability can be attained by the following more complicated compensation 
process which is required because of the fact that the strains of both test bodies 
are not only unknown, but are also entirely independent of each other. 
Since the agreement is attained when either 4°? or 4% vanishes, we 
need only consider the case where two component particles P, and P,, from 
each of the test bodies I and II influence each other. It is desirable to con- 
Sider especially the field components over two different space-time regions I 
and II. Owing to the fact that the retardation of forces admits no direct con- 
nection between P, and P, in order to produce a necessary compensation, 
one can attach to the second test body system another component particle PE 


(1?) L. ROSENFELD: Niels Bohr and Development of Physics (New York, 1955). 
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whose distances from the centers of mass of the first and second test body 
are X, and X, respectively. Pr is located in the immediate neighborhood 
of P, and is joined to P, by a small spring whose spring constant is k,y: 
Even though by the exchange of short light signals between P, and Er 
the total momentum of the second test body system does not change, the 
total strain momentum of the second test body system does change; conse- 
quently we do not know the equation of strain-momentum conservation of the 
second test body system by the exchange of messages (*) between the subsystems. 
This fact, however, suggests that a rather complicated process will be necessary. 
If we denote by e (*) the strain of P,, in the x direction, one can also 
use the same radiation bundle which measures the strain momentum of the 
second test body to measure the strain momentum of P,, through a more 
complex apparatus by adjusting the light path appropriately with mirrors as 
mentioned in Section 4. Indeed, this radiation bundle which measures the 
strain momentum of P,; at the times tir and na has to be arranged to measure 
the strain momentum of P,y at # and #. By this more complex process one 
can attain the uncertainty of e which is exactly the same as that of e". 
also think that P,, is released from the rigid frame at the time ¢, 
and again bound to it at the time # so that during 7, one can measure the 
strain momentum transferred to the second test body system. By choosing 
the force constants of springs joining each test body to its own rigid frame 
to be k,, and k,,, the strain momentum — k,,(Xy 6° Xe0)T,X, will be 


One may 


+ 7 II 1 “a + nc P A 
transferred from P,, to P,, while ky(Xx& 0 — Ae vine will be transferred | 


from P, to P,, during the measuring time interval 7,. Then one finds the 


equation of strain-momentum conservation 


2 
6 pio 


i ni 1) (LT) 
mer Mai Egy Sri 
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(9.4) 


(CLÉ) 
x 
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where the summations Y and Y are extended over all component particles 


{ II 


of each test body. Because one may regard all the springs joining the com- 


(*) One may think alternatively that the exchange of messages is allowable only 
if eqs. (9.5) hold. 
(**) This merely means that by receiving a radiation bundle Py suffers the dis- 


III 
placement e! Xy. 
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ponent particles of each test body to its own framework body to have the 
same spring constants, k,, and k,, respectively, we can take kia and k,, out 
of the summations. Moreover, replacing > af by the average value Lî/12, 
we get 


c? = mal = 
D)’ 2 (1) (L1) (11,1) 
ni — Te = 12 MT (En 3 Ex =F Hi ) NE 
= 


2 
(7 la == a x3) Ti ee = ksa 4; ne Tien 


(9.5) 
e? | = = = 
dt (1I)’ (EEE)! (III)' 2 (IT) (IL II) (LID 
ot eg ee iS METE + ha + En) — 
12 


-(È ese Tees? re bite") — ksa Xi Xn ree . 


We shall now consider the case of partially overlapping regions. There are 
three kinds of adjustable springs; namely, one of them connecting the com- 
‘ ponent particles of the first test body with those of the second test body, and 
the other two joining the component particles of each test body to its own 
framework body. The former kind is confined only to the overlapping spatial 
region and is attached only during the overlapping time interval. Hence in 
this case we obtain 


1)” (D' 2 (1) 71-1) LCL, I) 
Tig ST ee Mi TER + BY? + EMP) — rae? — rage, 
(9.6) 
; ; Ca = = = 
(CRI CUT) De 2 (1) GLI) (LID 
Ty Ty == 12 MELE 4 EX sn En ) Puri MEI TL Aia , 
dl 


where x1a, #24 and %,, are dependent on the geometry of each test body and 
in addition are functions of three independent adjustable spring constants 
kias Koa and ka 


We now turn to eq. (9.5). We choose the spring constants 


4 


. = 
tn = egy M Mobili Tr(Afti + AD), 


ext = — 
lio = 15 MEAT ASR A 5, MMI f (1+ 32) Au + Aah, 
Sad 


and an analogous expression for k,, by interchanging the indices I and II in 


34 
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the expression of k,,. By virtue of (6.16) we get 


MEDI (D 
Ty Ta MAT ES +5 SS ll My Li LT: Tr {€ oi CAE Un — ASSO) 
apr GT)’ Ù , 
(9.8) IO i he gii — Dole Tela E 


es DS M Mn TELAT Ta fe® ( ATI AS 


We therefore obtain two uncertainties 


zx 12h = 
ABI TAO + ggg Mn En Tue] A4 AUD | 
(9.9) 
771) 12h | 2 (1) (LI) 4 (ILD 
A fa Peat Mn Li Ty Ae™ | 53 MI? The; (AS -Aui n 


from which, according to a process similar to that discussed in the previous 
section, we find the minimum value of the product of two uncertainties 
(9.10) ABS ABS? | Asm — AND |, 

(9.10) is in complete agreement with the first expression of (3.6) which has 
resulted from the quantum gravitational formalism. 

In the case of two test bodies having a partially overlapping region, we 
can as well attain the result (9.10). However, in this case we have to connect 
springs joining the component particles of different test bodies to each other 
only during the overlapping time interval, and consequently %,,, %2a and %34 
in eqs. (9.6) are also functions of this overlapping time interval. 


10. — Measurability of two averages of different field components. 


When either of the field components is of the magnetic type, the measur- 
ability is completely the same as that described in reference (*). Hence we shall 
discuss merely three cases: when either of the field components is of the elec- 
tric type, when one component is of the electric type and the rest are of the 
magnetic type. To elucidate this statement we shall study the measurement 
of the weighted average of H,, over the region I and the weighted average 
of E,, or the ordinary average of H,, (or H,,) over the region II. 


Usi 
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Let us start with the first case, namely, when either of the field compo- 
nents is of the electric type. For either of the test bodies covering the space- 
time regions I and II respectively, the equations of strain-momentum con- 


servation are given by 


(I (I)! 
Ne ah — 


Bo 77 (1.1) TI(II,1) 
IiT:(Ei + His + ES ), 


(10.1) 
n ' Ga 
Ta et 75 Mr LA Tu(Es Bos ES »11) 1e Fe 


Ft If we do not take into account the correlation between a component particle 
x of the first test body and one in the second test body, by suitably adjusting 
17 the spring constants for springs joining each test body to its own rigid frame 
we have, by a calculation similar to that used in deriving (9.3) from (9.1), 


n À (10.2) AE® NEO ANA?) se 


1122 


1122 


which is, in general, in contradiction with the results from the quantum gra- 
vitational formalism and agrees with it only in special cases in which either 
A“) or 40 vanishes. The result of the formalism requires that when two 
ay. : time intervals 7, and 7, coincide the different components of the same type 
1. of field (*), say E and E, can be measured correctly, but (10.2) does not 
vanish and hence disagrees with the formalism. 
A | Now we choose from each of the test bodies I and II two component par- 
# à ticles, P, and P,,. Two radiation bundles of photons are arranged to measure 
Ba” the æ component of the strain momentum of the first test body and the y 
> component of the strain momentum of the second test body. We can, as be- 
fore, choose another component particle P,,, attaching it to the second test 
| body system, to be located in the immediate neighborhood of P,: Py, is 
Da Joined to P, through a lever system which consists of two equal and perpen- 
sà dicular arms, one each in the x and y directions, and is rotatable about its 
joint. The joint of the lever is bound to the rigid frame. P, is connected to 
the end of the x arm of the lever with the help of a spring which is itself in 
the y direction, while P, is connected to the end of the y arm of the lever 
through another spring, acting in the x direction and having the same force 
constant ks, and length as that of the first spring employed in connecting P,. 
We use the same radiation bundle which measures the strain momentum 
of the second test body at times #, and #, and also Pi at times t, and # so that 
the uncertainty of e" is exactly the same as that of «0. If we denote by X; 
and Y,, the distances of P, and P,,, from the center of mass of each test body, 


CE D CR 


(‘) The case discussed in (3.18). 
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then CHERE the time interval 7, an amount of strain momentum k (AT eu — 
TO < 
x) will be transferred from P,, to P, and an amount — heey En — 


D 7 
Xe )T; Yi from P; to Py. The equations of strain-momentum conserva- 
tion tee each test body system are then 


9 


nio: an co = a M 137,( Ei BO as Ee EN AI “ONE ae 
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(10.3) 
hae a, ae + i Er TINTE SL oy Mi TETE i ae ROSIE + Bo) ut sd 


Cala sl Vues = Aa ey) IAT RE > ko DE aT ey ee 


By referring to expression (6.16) and using suitably chosen spring constants 
ky, ky» and ky, we have 


2 


Gi ct = 
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from which follow the uncertainties of field averages in each region 


24 0 12% c? > 
AE» ping ta 2: ZIA M LE T. ESE AS +11) — AUD 
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(10.5) 
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and furthermore the minimum value of the product of uncertainties is 


a WO À FD _, Y I) 405,1) 
(10.6) AB AE, hl Aie Arisa | ’ 
which is in harmony with the result of quantum gravitodynamics. 

For the partially overlapping spatial region a discussion can be carried 
out as has been done after expression (9.6) and an agreement can be similarly 


attained. 
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Now let us consider a simpler case: the measurability of E, and Hu. 
The equations of strain-momentum and momentum conservation are given by 


ao” =i! pes N MESI (ED DE E&Y DE Eo) ; 
(10.7) 


AD” 


(D (II) ry (11) ry (11.11) (1.11) 
Pew an Pa = 625 le (Hier eno); 


where AZY, etc., are weighted averages while AH, etc., are ordinary 
averages as defined by (7.12). We get AH! from the corresponding equation 
by making I and II coincide. The self- “Ro of E” can be cancelled as — 
before by choosing k, = c?M,A%?/L, for the force constant of springs joining 
the component particles of the first test body to its own framework body. 
The self-reaction 7° resulting from x, the displacement in the « direction 
of the second test body system (which is considered to be a rigid body), can 
also be compensated for during the interval of measurement 7, by binding 
the component particles of the second test body to the fixed frame through 
springs having the same force constant k,,, which is chosen to be e220" 7,002 /Va 
with OS)? given by an expression analogous to (3.14), so that the force exerted 
on the second test body is equal to — e? LO" 7, 004,1. Omitting each self- 
reaction H°” and Hi from each equation of (10.7) and referring to (6.17) 


and (7.13), according to which E° and HT vanish, we find 


ponte CM BS, 
(10.8) 
Tani DEI ae CAT en A 


from which the minimum value of the product of uncertainties is immediately 
given by 


(10.9) AED AH ~ 0 


in agreement with expression (3.11) derived from the formalism. 
Let us finally consider the measurability of Æ,, and H,,: The equations of 
strain-momentum and momentum conservation are 


RE = = = 
nia = MEAT (ER + BS? + BE), 
(10.11) 
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As previously, if we join only the component particles of each test body to 
its own framework body by springs of suitably chosen force constants ky, 
and k,,, we get 


È ra ele = = 
oe ai Li Gato. M:L31 (ER PU II my Bi) ’ 
(10.11) 


y 


ma G2 = 
(II)” (II)’ (II) (II) 2 (I) D(ALII) 
== 0% == cd, Ta (#9 a 12 Mille, Bap) ; 


where B% is defined by (3.9) and BID is defined by taking the ordinary aver- 
age value over the space-time region II and the weighted average value over 


the space-time region I as 


~ —3K 0? 
CLI __ 4 4m! V2 Drt(x—x'). 
(10.12) ere PT LT. d ofa @' 09 (&)(diaV*) <a D(a—a') 
I I 
From (10.11) it immediately follows that 
(10.13) AHS Aig A [| Baa |+| Bae |], 


which again, in general, disagrees with the result from the formalism. From 
expression (3.7) of the formalism AE,, AH, vanishes when two volumes Vi 
and V,, coincide even if B% and Bam individually do not vanish but (10.13) 
differs from zero. 

The agreement can be attained as follows: as in the foregoing discussion 
we choose from each of the test bodies I and II component particles P, and P, 
respectively. To the second test body is attached another component -par- 
ticle P,, whose momentum is measured at times EE by means of the 
same radiation bundle which measures the momentum of the second test body 
at times #, and t, so that the uncertainty of the position of P,y is exactly the 
same as that of P,,. We join P,, through a spring acting in the y direction 
and having the force constant k,, with the end of the x arm of the lever as used 
in a previous discussion of this section. P, is connected to the end of the y arm 
of the lever with the help of a spring which acts itself in the x direction and 
which has the same length and force constant as that of the first spring. During 
the time interval 7’, the strain momentum — k,,(y" — e? X,)X,7, will be trans- 
ferred from P,, to P,, while the momentum ks (y™ —e®X,)T, will be trans- 
ferred from P, to Py, where X, is the distance of P, from the center of mass 
of the first test body system. Regarding Py as part of the second test body 


10 - Supplemento al Nuovo Cimento. 
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system, we have 
x 
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We now choose the spring constants 
bye = PLM KDE? + Mr Mn Ta TA (AGP + AP), 
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where 
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Equations (10.14) now take the form 
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(10.17) 
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Hence the unknown strain in the x direction in the first equation of (10.14) 
can be completely compensated for by these suitably chosen force constants 
of the spring mechanism. The unknown displacement in the second equation 
of (10.14) is similarly compensated for. The uncertainties of the field averages 
are evidently 
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from which the minimum value of their product is given by 


(10.19) AEDAHSD oF |Bam_ fun 


1122 ? 1122 


4 


which is in complete harmony with the quantum gravitational formalism. 

The measurabilities of #,, and H,, and of Z,, and H,, can also be attained 
similarly as follows. Using a lever having two equal arms which make an 
angle of 45° with each other, springs having the same force constant and 
length are connected perpendicularly to the ends of the arms. The lever is 
again free to rotate about its joint. The other free ends of the springs are joined 
to the component particles P,, belonging to the second test body system and 
to a component particle P, belonging to the first test body system. By this 
spring mechanism one may pursue a similar discussion and attain the results 
in agreement with the formalism. 


11. — Measurability of two field averages with arbitrarily oriented axes. 


Up to this point we have developed the measurement of the gravitational 
field strengths by considering two cubical test bodies of dimensions L whose 
three principal axes are parallel to each other. The principal axes of a cubical 
test body are defined as those which pass through the center of mass of the 
test body perpendicularly to its surfaces. We have limited our considerations 
on the measurements of the test body which is uniformly stretched in the 
x direction to the measurement of the #,, component. We have also limited 
our considerations on the measurement of the test body having the x com- 
ponent of spin density o, to the measurement of the H,, component. 

We shall extend our treatment to the more general case in which it is not 
necessary for the principal axes to be parallel to one another, but in which 
the properties of the test body are kept the same. For this purpose it is con- 
venient to consider merely the #,, and H,, components and to rotate one test 
body keeping the other fixed in the results of the previous treatment. The 
principal axes are assumed to serve the purpose of body co-ordinate axes for 
each test body and will be denoted by x; and x. The direction cosines of 
the j-th axis of each test body may be represented by a,; and b;;. Hence 


ML. 1) D Gal, 5 Di = bi, 0, 
In addition to the orthonormality conditions 


(11.2) dij ig = Oye 5 bi Din = Oix y 
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the cosine and sine of an angle between the j-th co-ordinate axis fixed in the 
first test body and the k-th axis fixed in the second test body may be ex- 
pressed by means of the product of the a’s and b’s as 


(11.3) COS 57 = Gi5bir , 


(11.4) Ay SN 0,x = ep;14;5dw » 


where À, is a unit vector perpendicular to either axis. It is a well known 
fact that cos6;,, can be expressed by the Euler angles since we are concerned 
merely with the orientation of the body axes relative to each other. 

It is desirable to start by considering the commutators between two gravi- 
tational field averages of the magnetic type. We properly rotate the body 
axes of one test body with respect to the space axes in which the second test 
body is fixed. According to (11.1) one finds from (3.3) 


(11.5) [BERT HO = (Sid 


tg 2 pa ijp q iva) 


with 


(11.6) Kisa 413 a a Tae fe bs: (Dip Dia 1 Die Dig — Diz Dy a) D™**(a' iis 


where 
GD = ©? COS Os! ©? 
7. SARE RIO ce oto’ 
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and an analogous expression for 2... 
For the commutators between electric-type components we have 
Wd)’ Eu (I de T OLD 
(11.8) ee doe Ale [LE LE. 
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with 


(11.9 EDI + 4 7 
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where we have considered only weights in the x and x" directions for each 
test body as mentioned at the beginning of this section. 

The commutators between an electric- -type component and a magnetic-type 
component can be obtained by a similar procedure. We rotate the body axes 
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of the first test body to arbitrary directions signified by 7, j and k, (not pa- 
rallel to the space-fixed body axes of the second test body) relative to the 
body axes of the other, which are signified by p’, q’ and 7’. By virtue of (3.2) 
and by the same method as used in deriving the commutators between two 
magnetic-type components, we get 


p20) Pe, Fron) = an Gan, dn) 
with 
3K È ; 
Jalal Gap, A tel G4xl tay: 
( ) jp a Seems | x | d*x w(x) 
si LE 
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Here the expression is obtained by taking the weight in the a direction in 
the body co-ordinates fixed to the first body, for which we measure the electric- 
type field component. 

Following the procedure used in deriving (6.11) we get (*) 


SM Ke? j 
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T 
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Similarly for the magnetic-type components we get 


~ KM Ae? 
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In eq. (11.12), (11.13) we have considered only the field uncertainties in the 
region II resulting from the uncertainty of e” in the region I. 
Putting îi=j=1 in eq. (11.8) and (11.10) and making use of the Un- 


(*) See Appendix II. 
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certainty Principle, we obtain 


DEA i l Rei cr 
HA) APN IR EE fase 42! (0, €')(2D 9 Prat — Dia Ly) DE —%), 
mi 


LT LT ne | 

I 

= 4 3ih 
E OY A FID’ = i FAQ 
(dal) AK AH yg ~ one) al dia 


© © 0 
9G % n : j, NENTI a ret fn! / 
(22, AVRO UE Dia COS On" Ep'atr —— leon D**(4'— x). 
UN 1 n 


Comparing eq. (11.12) and (12:13) with (11.14) and (11.15) respectively, we 
find 


o 12h 
(I° SERRA ae ES 
ene) ME MA 


Hence we have arrived at the conclusion that one is required to measure only 
the weighted average value of the £,, component. 

A similar treatment can’ also be carried out for magnetic-type components 
on the measurement of two cubical test bodies whose principal axes are not 
parallel to one another. Therefore we reach an analogous conclusion, namely, 
that one needs only to measure the H,, component. The measurability of two 
test bodies whose principal axes are not parallel to each other can also be 
attained by an argument similar to Sections 9, 10 and 11 if we replace the 
lever with two perpendicular arms by a lever with two arms which have a 
properly replaced angle between them. 


12. — Conelusions. 


We have shown that the limitations determined by physically possible 
measurements are in complete agreement with those imposed on the measur- 
ability of the quantized gravitational field strength. It has also been shown 
that in order to measure the field strengths one is required to take a field 
average over a space-time region. For magnetic-type components the meas- 
urement procedure and field averages are similar to those of the electric mag- 
netic field strengths while for electric-type components field averages are hyper- 
bolic weighted averages and for the measurement of electric-type field averages 


> 
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it is necessary to measure the strain of an elastic medium which serves as a 
test body. The disturbances of the field averages in one spatial region caused 
by the measurement of some field components in the other can be obtained 
from Einstein’s equation. These disturbances completely agree with the results 
from the quantum gravitational formalism taking into account the dynamical 
equations of motion of the test bodies together with the limitations imposed 
on the measured quantities (the strains and strain momenta for the electric- 
type component and the co-ordinates and momenta for the magnetic-type 
components) by the Uncertainty Principle. The self-reactions of the test bodies 
during a measurement-time interval can be compensated for by introducing 
appropriate spring mechanisms as Bohr and Rosenfeld did in the measurement 
of the electromagnetic field strengths. 
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In the quantum electromagnetic formalism if there is a uniform polariz- 
ation 0Ax per unit volume in a spatial region I the uncertainties of field 
strengths at (4%, Yo, 22, &) are given by 


(A.L1) APyy(&,) = | [(Age) ue — (Aju) y] D™ (ar — a) dt, . 


Without finding four potentials as derived in reference (?), we can easily com- 
pute these uncertainties by means of the method described in Sect. 6 by using 


oP C1 : : 
p ra ì = Res Ai. = Ajs =) , 
(A.I.2) Aj, ae? ATEN n Je f 


where 


j ON eS DE aoe 
P= phv0(x + 5)0(5—#)0(y +5) 9(5—v) ol: AL ou Ott). 
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APPENDIX II 


According to eq. (6.7) we have 


sii [ 7 UR ee Cso 
(A.II.1) AE pa(0') = DI dix (ar. + Bae æ 0%! HD 2AT op om! + 
Osa 0 ret 
+ (ate —22an) 2] (x'— x) . 


By rotating the three principal axes of the second test body through arbitrary 
angles with respect to the space co-ordinate in which the first test body is 
fixed, three direction cosines of the new i-th body axes are given by b,;: 
Multiplying the above equation by b,b,,, summing over p and q and denoting 
ABsabnba by AE, we obtain 


(A.IT.2) AE, (x') = — à fat (2AToo + ATI) = ONT E 7+ 
Oar’ ‘eût 0x, 
ai (2AT;, — AT) Ea Dee 


Substituting ATio=b1Tio, AT1=bab1T;; and relations (6.6) into (A.II.2), 
we get 


(A.IIT.3) AE:(0) = — 4 COS Ou = 


c20x Ot Ox! Ot! 
OH ON OR A, 
mor ar (an dn cer: 


Es 0? GET ANT QUE 


ss 2 COS? Ch 


CF ob 
Referring to the derivation of (6.10) we can write 


KM Ae, 


113 [areorte) D )(2 Dis Dy — Dy Dy) D™*(g'— x), 


ï 


(A.II.4)  AB,(x') = 


from which, by taking the weighted average over the region II, we have 


3KM Ae (I) 
(A.IT.5) AG — a Tu, i dix i dt! co(a, 22 Dr Dy — Da Dry) D™(@'—2). 


I 


Generally speaking, AE” reduces to (11.12). In a similar fashion one can 
easily compute the magnetic-type components induced by an uncertainty in 
an electric-type component in the region II, and the field uncertainties as well. 
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SUPPLEMENTO AL VOLUME XXI, SERIE X N. 2, 1961 
DEL NUOVO CIMENTO 3° Trimestre 


Selected Topics on Analyticity in Potential Scattering. 


A. MARTIN 
CERN - Geneva 


(ricevuto il 22 Luglio 1961) 


Contents. — 1. Introduction. — 2. Analytie properties of pariial waves 
for Yukawa-like potentials. 1. S wave case. 2. Higher angular mo- 
menta. 3. Case where tensor forces are present. — 8. Application of the 
partial wave properties for generalised Yukawa potentials. 1. Numerical 
calculations. 2. The deuteron wave function and the anomalous thresh- 
old problem. 3. The connection between the unphysical cut of a partial 
wave amplitude, the S matrix, and the potential. — 4. Analytie properties 
for fixed energy. 


These lectures were prepared for the summer school of Theoretical Physics 
in Yugoslavia (1961). They are essentially based on work done by the CERN 
group, the Torino group and the Padua group. Though, in the text, reference 
is made to other lecturers, these lectures can be considered as self contained. 
They do not cover all the subject but the interested reader will easily find 
in the references what is missing. 


1. — Introduction. 


In recent years a large amount of work has been done on the analytic 
properties of scattering and production amplitudes in field theory. However, 
at the time where these studies were started, very little had been done about 
the analytic properties in potential scattering except in the case where the 
interaction took place in a strictly finite region (1). The first really important 


(1) E. P. WIGNER: Phys. Rev., 98, 145 (1955). 
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contribution in this field was that of KHURI (?) who showed that the scattering 
amplitude derived from a Yukawa-like potential obeys à fixed transfer dis- 
persion relation. 

More recently, assumptions have been made about the analytic properties 
of the scattering amplitude in field theory, both in energy and momentum 
transfer (*). There does not exist yet any complete proof that these assumptions 
(Mandelstam representation) are valid. However, in this case, the potential 
model is a useful guide which shows that at least for generalized Yukawa 
potential the Mandelstam representation holds (4). 

We shall not speak here about the complete derivation of this result which 
has been gyen by Prof. T. REGGE here. We shall first consider the analytic 
properties of a given partial wave. These analytic properties would be obtained 
by projecting the Mandelstam representation on a given angular momentum 
state. However, a direct method of obtaining the analytic properties of à 
partial wave proves very useful. Indeed, it turns out that the method we 
shall follow for Yukawa-like potentials furnishes many other results. In par- 
ticular it furnishes a very nice integral representation of the wave functions 
(especially of the deuteron wave function) which permits to study certain 
matrix elements such as deuteron electromagnetic form factors and deuteron 
photodisintegration and to show that the so-called anomalous thresholds 
(at least the « normal » anomalous thresholds) are connected with the extension 
of the deuteron wave function. 

In addition, the discontinuities of the 7°-wave scattering amplitude across 
the unphysical cut can be expressed explicitly from the knowledge of the 
potential. It turns out that the reverse is true: given the discontinuity one 
can construct explicitly the potential. Since the discontinuities are at present 
the quantities which the dispersionists calculate first one sees that for a given 
partial wave one can associate a potential and, therefore a wave function. 
In this way one can connect the information about the discontinuities for 
nucleon-nucleon scattering and the external part of the deuteron wave 
function. 

A similar connection between potential and discontinuity can be established 
in the case of fixed energy dispersion relations. We shall have to sketch first 
a method to establish the analytic properties for fixed # which is completely 
different from what you heard from Dr. REGGE. Then we shall be able to 


(2) N. N. KHurI: Phys. Rev., 107, 1148 (1957); see also: A. KLEIN and C. ZEMACH: 
Ann. Phys., 7, 440 (1959). 

(3) S. MANDELSTAM: Phys. Rev., 112, 1344 (1958). 

(4) A. KLEIN: Journ. Math. Phys., 1, 41 (1960); R. BLANCKENBECKLER, M. GoLD- 
BERGER, N. KHuRI and S. Treiman: Ann. Phys., 10, 62 (1960); T. REGGE: Nuovo 
Cimento, 14, 951 (1959); J. Bowcock and A. MARTIN: Nuovo Cimento, 14, 516 (1959). 
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show how one can construct from it the potential. As an application we shall 
show how one can settle the problem of the range of the imaginary potential 
in nucleon-antinucleon scattering. 


2. — Analytic properties of partial waves for Yukawa-like potentials. 


_ 21. S wave case (5). — In the preceding lectures Prof. T. REGGE has shown 
to you the analytic properties of finite range potentials and, also, established 
the existence of an analyticity strip for potentials decreasing faster than an 
exponential. We shall consider here the analytic properties for potentials of 
the type 


(1) Vir) = | O(a) exp[— ar]da, 
le 
which include the Yukawa potential (C(«) = const). 


It proves useful to consider the following two solutions of the Schrédinger 
equation 


(2) a SE eon 7 ees 

uh, 7) = exp kr] f(—k, 2), (CONTE 
(3) with 

Uh, 1) == exp [= kr f(k, 7) , (co) 


whether or not such solutions exist for complex values of % will be established 
by trying to construct f. Notice that f(k, 0) is just the Jost function about 
which you have already heard. 

If the function f(k,r) is known, one can easily construct the S matrix 
(I should say S functions), defined by 


U 


segui — EXP [— thr] — S(k) exp [ékr] , 
(for large r) at least for real k. 

Since 
ut(k, 7)  u-(k, 7) 


Ureguiar = ut(k, 0) u-(k, 0) 


(5) A. Martin: Nuovo Cimento, 14, 403 (1959). 
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we see that 


ay — HE 0) 
(4) DSS OMY. 


For complex values of k, the latter expression will be maintained as a 
definition of S(k). 

Let us now speak about the construction of f. 

In the especially simple case of a pure exponential potential (°) (Co 
=6(a—j)) the solution of the S wave Schrédinger equation is known and 
one can easily recognize that f(k,r) can be expanded as a series of expo- 
nentials of the type 


f(r) =1-+ Ÿ CE) exp [— mur] , 


n=1 


for a potential V = V, exp[— ur]. 
The same can be done for more complicated potentials 


M 
V=)> Vrexp[—uyr]l, 


N=1 
for which no exact solution is known. Here one can write 
EX eee eau 
This strongly suggests that when V(r) is a continuous superposition of expo- 


nentials in which the largest range is 1/u, f(k,r) can be written as 


(5) {(k,9)=14+ | os(a) exp [— oer] der. 


7 


We shall take this formula as a starting point and check later on that our 
assumption was correct. 
Inserting u~ = f(k, r) exp[—ikr] in the Schrédinger equation we get 
(6) f—2ikf — V(r)f =0. 
Replacing V and f by their integral representations (1) and (5) respectively, 


and using the well-known theorem about the product of two Laplace trans- 


(5) R. Jost: Helv. Phys. Acta, 9, 256 (1946). 
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forms (which appear in Vf), we end up with an equation for o,(x): 


(7) a(œ + 2ik)0,(2) = O(a) - | ot B)ox(B) db . 


Now we have to fix the upper and lower limits in the integral. Obviously 
the lower limit is w because O(~)=0 for «<p. From this fact we easily see 
that @,(«) = 0 in the range «= 0, « =, which is consistent with the assumption 
we made. Finally we can write 


au 
(8) a(a + 2;k)or(x) = C(x) + O(a — ZINC — B)ox(B) dp . 


(22 


This integral equation has the great advantage that it can be solved very 
easily by iteration for any finite value of «: 


(Orso ir o=0, 
C (x) 


HN ENT di 
ATL 
pa a _ (Cle = B)C(B) 
dana = (0 +] ta 4 


IZ) 


and, more generally, if 0,(«) is known for «<n one can insert it in the right- 
hand side and obtain o for x<(n+1)u. 

It is clear that if (x+2i%) cannot vanish, i.e. if k lies outside a line on 
the imaginary axis for k=7(u/2) to ico, o,(«), for fixed «, is a well defined function 
of k, is analytic in k. We have now to see what kind of condition should be 
imposed on C(«) to ensure that these properties hold for f(k,7). We shall 
look for sufficient but not necessary conditions. 

Let us assume that there exist a monotonous increasing function B(x) 
such that 


(9) | (a) |< Bla) . 
Let us take & outside the cut on the imaginary A Pe 
axis; more precisely we exclude the domain in- as RARE 


dicated on Fig. 1, where & is an arbitrarily angle 
then, for a>wu 


a + 2ik|> a sine 
| 


la +2ik|>2]|k|sine, Fig. 1. 
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then we get 


(10) | ex(%) |< œ? sin € 


|] 6196]. 


Let us call 


1+ | lox(B)|4B = (À) ; 


then 
y! _ Bla) 
y «sine 
Hence 
“aBB(p 
yx(x) < exp B28 sin e 2 


if 


It is clear now that if fi dBB(B)/B? is bounded 1+ fet «)|da is uniformly 
lt 

bounded with respect to k when k is outside ae neighbourhood of the 

cut. Therefore a sufficient condition for the existence and holomorphy of 


f(k,r)=1+ Jen [— œr]ox(x) dx, for fixed arbitrary r with respect to k is 


(11) | O(a) |< Max, n arbitrarily small: 


another sufficient condition would be 


Ma 
2 | | i —— n 
(22 LEE À + IntM(œ/u) 


This ensures, in particular that f(k, 0) is regular outside the cut i(u/2) — à co. 
Further, making use of the uniform convergence property we see easily that 


since 0;(x)-> 0 as k + co in any direction except the positive imaginary axis 


(13) lim f(k,r)=1. 


1k|—> 0 
Further, in the case when condition (11) is fulfilled a more careful examination 


n 
Le) 
re 
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shows that 


[f( ry — De [ki LOL: meer sl 
In |K| 


Ag ale SEI: for "7 17 
[| \k| OI n >1 


(14) 


The conditions imposed on ©(x) essentially mean that V(r) has to be less 
singular than 1/r? at the origin of coordinates. We know that if this is not 
so difficulties arise, and in particular for sufficiently strong 1/r? attractive 
potentials the Schrodinger equation has no regular solution. The question 
whether one can deal with very singular repulsive potentials using our method 
is not yet settled, but we would like to show that when C(«) is bounded by 
an arbitrary polynomial f(k,r) retains all the desired properties for r +0. 
In this case 


(15) | C(a) exp [— ar], 
is bounded by, say, 6; then 


[22 
” 


1+fox(8) exp [= fr dé : 


Lu 


B 
(16) ox) exp [— ar] |< ——— 
æ sine 


hence one shows easily that 
1 +| lox(a) exp [— ar]|da, 
Lu 


is uniformly convergent and all the above mentioned properties hold. 

Up to now we have established the analytic properties of the outgoing 
and ingoing wave functions for fixed r. In case condition (11) or (12) are 
satisfied it is an almost trivial matter to get the analytic properties of 


k, 0 
S(k es A ? je 
i k, 0) 
It has two cuts from k =ipu/2 to k—ico and from k=—iooto k=— 1/2. 


To show that these are real cuts we should establish explicitly the existence 
of a discontinuity of S(k) across these lines. This will be postponed to Sect. 3. 
In addition to these cuts S(k) has poles which come from the zeros of f(— k, 0). 
You have probably already heard from Regge that the zeros of f(— k, 0) in 
the upper half % plane necessarily lie on the imaginary axis (*) and correspond 


(‘) This is due to the reality of the potential. If f(—%,0)—0, with Imk> 0, 
u+(k, r) is regular at the origin and vanishes at infinity. Hence by combination of 
(d2/dr2--k?)u= Vu and (d2/dr?+k#?)u*= Vu* one gets (k2 k**) futudr— 0. 


Del 
Uci 
~ 
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to bound states whose wave function is just 


(7) exp [— xr]|1 +] ac) exp [— ar]da|, x real > 0. 
Li 


From f(k, 0) —>1 for k — co we get S(k) >1. Using the asymptotic behaviour 
of f(k, 0)—1 we can show that, provided condition (11) is fulfilled, one can 
write wnsubtracted Cauchy integral representations (or dispersion relations) for 


S(k) —1 


= exp [10] sind. 
Finally, taking into account the symmetries of f: 

(18) f(k, 0) = f*(—k*, 0), 

we get the following symmetries for S 

(19) S(k)S(—k) =1, S(k) = S*(— k*), 


and the results are summarized in fig 2. 
As a last comment of this section 
we would like to say that when O(a) 


a 2\ . . . e 
mesa increases like an arbitrary polynomial 
ufo f(k, 0) and f(— k, 0) do not exist. How- 
k plane ever, we believe (without being able 


to prove it at present) that 


SUR Dit) 
= ST cay nen 


might still exist if C(«) is positive (re- 

Fig. 2. pulsive potential at short distances) for 

large a. The analytic properties would 

remain unchanged except perhaps the behaviour at infinity. Our belief is based 
on the remark that for a hard core potential: 


Vit) eo, for Le 
(20) i 
V(r) =e) exp [— ar] da for EE To) 


LU 
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where C(x) is bound by a polynomial the analytic properties still hold, but, 
at infinity, S(k) ~ exp[2ikr]. Then, we believe that letting » go to zero 
will not seriously change the picture, or, in other words, since here we have 
replaced the singular region by a worse singularity in the potential, nothing 
worse than what happened in the hard core case can happen. 


22. Higher angular momenta (7). — Here all the trouble comes from the 
centrifugal terms in the Schrédinger equation 


(21) ne — + k2 — rog, =), 


In this case we have two possibilities: 


— either to keep that same integral representation 


ao 


| 


0] 


(22) u-(k, r) = exp [— ikr] È ~ | 0,(%) exp [— ar] da 


but then, since the range of the centrifugal potential is infinite we expect 
that the lower limit in the integral representation will be zero rather than w 
and the iterative property of the integral equation for o is lost at least in 
appearance. 


— or to try to generalize the S wave case by using functions which satisfy 
the free equations (with centrifugal term included), In the S wave case we had 


u-(k, r) = exp [— tkr] + | 0x(x) exp [— (tk + a)r]da, 
u 
in the 10 case one can take 
(23) u-(k, r) = kr HP) ACC + ia)r]da|, 
| È 


where À is defined as in SCHIFF: Quantum Mechanics, p. 79. 
The latter method which was used by DE ALFARO and ROSSETTI (8) avoids 
the difficulties associated with the centrifugal term itself but, one has to be 


(7) A. Martin: Nuovo Cimento, 15, 99 (1959). 
(8) V. De ALraro and M. RossertI: Nuovo Cimento, 18, 783 (1960). 
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able to write the product w-(k,7)V(r) as a superposition of Hankel functions, 
which is possible but rather tedious and difficult. 

Here we shall prefer to use the first method and show how one can get 
rid of the centrifugal term. Let us first merely include the centrifugal term 
in the potential. Since 

a 
1 


02, 


—|lxexp[-ar]de, 


= 


0 


it is not difficult to see that starting from 


eni + fot exp [— ar]da, 


one gets 


(24) x(a + 2ik) oa) = Cla) + WL De + JC B) + Ul + 1)(x — B)]0x(R) dB. 


What one would like to do is to transform this equation into an iterative type 
integral equation. One can notice that when C(x) —0 the solution of this 
equation is known, because it can be obtained by means of Laplace trans- 
formation from the Hankel functions. We shall prefer to show it directly. Let 
us call o°(x) the solution for C(«) = 0 


(25) ala + 21k) 0} (ae) = Ut + 1) 


by differentiation of both sides of this equation we readily get a differential 
equation for the quantity 


yi(a) == bE | on(B\ dd, 


[oc(x + 2ik)y® (= + 1)y%(z) : 


if we compare this with the Legendre equation 


(27) —[(L= mu |= Ub +1)u , 


6 


715 


SELECTED TOPICS ON ANALYTICITY IN POTENTIAL SCATTERING 167 


we readily see that P,(1+(a/ik)) and Q,(14(«/iX)) are solutions of the equation 
for y;(x). 
Actually, since 0?(«) is regular at «= 0, and since y,(0) =1, we conclude that 


>, (1-4 (a/ik 
(28) EE): 
da 


Let us notice that P(1+-(a/ik)) is regular in the whole complex k plane, 
with a pole at k=0, which means that even for k =ir the equation for oy, has 
a well behaved solution. It will be useful to notice that the other solution 


Q,(1+(a/tk)), for fixed x, has — and it can be shown from 
Stal 
Lei ydt 
Gala) = 5 | a ’ 
fe 


a cut along k=iX where YX is a real number going from 1 to +00. 
Let us now return to the equation for o(x) in the case C(x) #0. It may 
be written 
x 
x fl b) 0:(8) dé = 


0 


(29) a(a% + 21k) o,(~) — U(14 1) 


= O(a) + {Ct — B)ox(B) 4B = g(2) - 
0 
Notice that g(a)=0 for «<y; so for «<p, ox) = 9 (x). We can as before 


use y(œ) = 1 + [o,(8)dB and get 
j Ra dg 


(30) [x(a + 2) el = Ul + yi) 


g will be treated as a known right-hand side according to standard methods. 
Then one can formally solve the equation and obtain 


D 
(31) yale) = P, (i elke: 
à | M DIR le | = | dg a 
+] P.(1 +) a(t t | 21L 1 i) gi {1 ik} | ap dp 


In this solution a choice of the homogeneous solution has been made: it is 


such that for «<p, Yy(x) = yî(a) since for «<p dy/da = 0. 
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At this point we could still transform a little bit the equation, but the 
main point is already exhibited here: 


x— pt 


pa) = O(a) + | O(% — B)yi(B) AB , 


0 


(because C(a)=0 for x< w). 

Therefore, if we know (8) for B<nu we know g(6) for 6<(n+1)u and 
inserting this in the r.h.s. of eq. (31) we obtain y,(6) for B<(n+1)u. In this 
new form the iterative property is re-established. 

One might think that a troublesome underlying assumption of our method 
of transformation of the equation is the derivability of g(a) especially at the 
point «= (which is not valid for a pure Yukawa potential). Actually, by 
partial integration and derivation one can transform the equation into an 
equation for o,(x) where (8) rather than dy/df appears. It is not difficult 
to believe that the lack of derivability of y(f) will not alter the final result. 
_ We shall not, however, enter into these details. From what we know about 
the analytic properties of Q and P it is clear that o,(x), again, is analytic pro- 
vided % is outside the origin and outside the cut k= iu/2 —ico. To prove 
that this is also true for f(k, r) is tedious, because one needs to introduce upper 
bounds of the various functions appearing in the equation, but nevertheless 
possible. All the results which have been established in the S wave case 
can be proved in the present case. To establish the properties of S(k) = 
lim f(+k,r)/f(— k,r), one can notice that S(k) is also Lui 0x(2)/0-x(2). 


23. Case where tensor forces are present (9). — Since the actual nucleon- 
nucleon forces are not central but contain spin orbit and tensor forces, it is 
desirable to extend our method to the case where L=J—1 and L=J +1 
states are coupled through a non-central potential. Then the starting equa- 
tions are 


d? JI) n 3 
| a ki re) Voir) | Walt) = Volr) vs), 
(32) 
d? , (J+1)(J +2) a 
Aviat Kees. DE = Ft) tee tn) Uri 


Let us remind first a few facts about scattering by tensor forces in general | 
because, though these results are well-known by specialists, they are rather 
unfamiliar to the non-specialist. 


(9°) A. MARTIN and R.VinH MAU: Nuovo Cimento, 20, 390 (1961). 
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When in eq. (32) the boundary condition w,_;(0) = %,,, 


the system admits 4 independent solutions. Here we shall select the inde- 
pendent solutions by their asymptotic behaviour 


(r) = 0 is disregarded, 


Ua | exp [ikr] exp [ikr] © ‘exp [— ikr] 
(33) = vi vi a | 
Visa exp [ikr] — exp [ikr], exp [— ikr], 
exp [— thr] 
We 


— exp [— thr], 


(Of course other linearly independent sets of solutions could be used). If we 
impose the boundary conditions at r=0O we still have two linearly inde- 
pendent solutions. 

The scattering matrix can be written as a two by two matrix: 


y J—1 Sg J+1 
(34) 


al SY. 
S y J—1 Sy J+T 


S,_,,-, measures the amplitude of the outgoing J—1 wave when the initial 
J —1 ingoing wave is a pure J—1 wave; S,_,,,, measures the amplitude of the 
J—1 outgoing wave when one starts from a pure J+ 1 ingoing wave, and 
so on. So to construct the scattering matrix elements one has to fulfil at the 
same time certain boundary conditions at the origin and at infinity. For 
instance to determine S,_,,_, and S,,, ,_, we should start from a combination 


(35) avi + av + yy + yy, 


where no L =J-+1 incoming wave is present. Then the L—J—1 asymptotic 
wave is 


(36) | (a, + a») exp[ikr] + 2c exp [— kr], 


hence, by comparison with —$,_,,_, exp[tkr] + exp [— dr] we get 


: Ar + Ay 
Sy = — ve ’ 
and also, clearly 
5 do = ay 
DFE AVE D0 ’ 
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once yi yp,y, are known, a,/¢ and a,fe are determined by the condition 


(37) apt (0) + ay, (0) + e[y, (0) + y, (0)] = 0, 


which furnishes two equations (*). 
The S matrix constructed in this way has to be unitary and symmetric: 


for real energies this means: 


Ù ro 
Pai SERRE 
& MERS Ips 
38 Spa gail yes Dot 1, 
(38) a ates men 
SEI ge |S 1 441 353, 
O ii En yu SR: 
PEG pat pa TOS Jaca 0. 


The latter equation is already satisfied in modulus if the three first equations 
are fulfilled. Therefore the 8 real variables appearing in the four matrix 
elements are connected by 5 independent linear equations and the scattering 
matrix can be expressed in terms of 3 real quantities. 

This is probably the reason why many people prefer to express the infor- 
mation on the scattering in an other way where no superabundant quantities 
appear. In the representation above considered the S-matrix connects the 
ingoing states 

exp [— thr] i 0 
and x 
0 exp [— 2kr] 
with the outgoing states 


exp [tkr] 0 
and 
0 exp [kr] 


One could try to choose another orthogonal set of ingoing and outgoing 


states such that the S matrix becomes diagonal, i.e. make an orthogonal trans- 
formation on S. We take the ingoing states | 


ae | exp [— ikr] Le 1 — n exp [— ikr] 
V1 + n exp | — ikr] ui exp [— ikr] i 
(39) and outgoing states 
Yen ATE exp [tkr] | re po md 
V1 +7? \n exp [ikr], VI+#\ exp [ikr] 


(") The effect of the centrifugal term will be considered more seriously later. 
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then the unitarity condition remains valid after the transformation and if 


we choose 7 such that 8,,—S,,—0 then |S,,|=|S,,|=1. n is determined 
by the equation 

(40) Spa ptt — 9) (8 jag ge — Sa) 

which expresses the fact that 7 can take two values 7, 7; such that 7,4,;=—1, 


which just corresponds to the two ingoing solutions we have considered. One 
can check, using the conditions from unitarity and symmetry that 7 is real. 
Obviously once one has 
Ù ra, = RO ‘ Y AY Dj, *\ 6 
(41) Der EXP 2t0.] and Sop = exp [2705] (©); 
with 7, one can reconstruct the S matrix in its original form if one wishes. 
For complex energies (if the S matrix makes sense), the symmetry property 


of the S matrix still holds, but the unitarity condition becomes S*(k*) S(k) =1. 
It is still possible to diagonalize the S matrix; the equation for 7 still holds, 
with 7,7, ——1 but 7, becomes complex, and | Sao |Sge|#1- 


To be complete we should add that in order to compute a cross-section, 
one has to adjust the amount of ingoing wave to what is present in an incident 
plane wave, and using the S matrix obtain the outgoing scattered wave. 

Let us now return to our specific problem where y,y,y,¥y, are obtained 
from potentials of the type 


fee ct exp [— ar]da, 
mn 
(42) Vu =| Coal) exp [= a7] de. 
ft 
= | C,(«) exp [— ar] da. 
si 


Then we try to represent the wave functions as 


‘fi (ky 1) | f.(k, 1) 
(43) wy = exp [— dr] ; is exp ihr] 4 
Gi(k, 7) ga, 7), 


yi and y, are obtained by changing k to — À. 


(*) dx; 0g are called « eigenphases »; #,, 73 are the « admixtures ». 
a? ) 
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We’ assume (as will be a posteriori checked to be right) 
a 


fit, r) = 1+fon(2) exp [— aride, 
0 
falkyr) = 1 + faute) exp [— ar] da, 
0 
(44) È 


gi(k,r)= 1 + fouls) exp[— ar]de, 


gG.(k,r)= — 1 + Joux) exp [— ar]da, 


0 


thenTwe insert in the coupled Schrédinger equations 


à di — Al Zi £ = 
| Je ee yey. 
(45) c 
} J +1)(J +2 Sal i 
jo" vi | a È e a 


The equations can be treated exactly in the same way as in the 140 case. 
It turns out then that they can be transformed to be of the form 


x 


(46) inl) = 081%) + ple) + | Ks te B)9(R)AB, 


0 


and so on. Where 0°(«) = (d/dx)P,_,(1 + (ajik)). w(x) is known and vanishes 
for «<p, and œ(B), vanishing for P<p contains oy(y) and ox(y) for argu- 
ments y smaller than 6 — u. The iterative procedure works provided one 
solves at the same time the equation for On (oe) and ox(x). Eventually, one finds 
for f,(k, r) if: the same analytic properties with respect to k as in the S wave 
case. 

A slight difficulty occurs when one wants to construct from figifeg, the 
S matrix. As it was sketched in the general considerations about tensor forces 
one has to take suitable combinations of Yi Ys, ¥, and to impose regularity 
conditions at the origin on the combination. Unfortunately, due to the cen- 
trifugal terms, f,, and fi» are singular at the origin. One would naively 
think that 


Dim fe (En) and lame 


r— 0 r> 0 
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are finite, and in the same way 


lim r7** 9,(k; 7) and lim 7777 


7 0 r— 0 


Jo(k, 1) . 


If this were true it would not be difficult to apply the regularity condition 
at the origin. This is just what is done in the 1/40 uncoupled case. However, 
it was shown by H. CORNILLE (1) that this behaviour at the origin is only 
correct in one case, which is when the coupling potential V,(r) vanishes at 
the origin. This means that when V,(0)= 0 there is no difficulty to construct 
Sos go Soi pa and Si; ,4,- It V,(0) 40 a serious trouble arises because 
one can find a set of four solutions, equivalent to the set yy ;p,yp,, two of 
which are regular solutions, the other two being irregular solutions with different 
degrees of singularity at the origin. Any kind of limiting process to deal with 
finite quantities will project out the most singular part of f,f: 919 and there- 
fore the less singular of the irregular solutions will not appear. So the number 
of equations to determine the matrix elements will not be sufficient. A more 
elaborate procedure must be used: one must first take into account the worse 
singularity and then subtract out from f,f.gig. the contribution of the most 
singular solution. This will give just the complementary equation needed. 
Let us notice that certain field theoretical potentials such as the Gartenhaus 
potential satisfy the condition V,(0)=0 so that the method can be applied 
without complication. 

In any case the matrix elements S,_, j-1) Sin gas Ogg 741 ave the same 
analytie properties as the S(k) function in the uncoupled case. They have 
the same cuts and further one can show that the poles in the upper half 
k plane lie on the imaginary axis and correspond to bound states. 

The analytic properties of the diagonalized S matrix and of the admixture 
parameter are not simple. This comes from the fact that the diagonal matrix 
element, as well as the admixture are functions with two determinations. 
However, one can construct directly from VY; viw7y, the eigensolutions: one 


starts from solutions 


(pr + Ays) — S(pi + Ay?) , 


these are automatically eigensolutions. The two regularity conditions deter- 
mine A (related to the admixture) and S. This may be useful in calculations 


for real energies. 


(19) H. CorNILLE: private communication. ‘To be published. 
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3. — Application of the partial wave properties for generalized Yukawa potentials. 


31. Numerical calculations. — The method we have developed in order 
to obtain the analytic properties of partial waves may be also considered as 
a practical method of calculation. Let us consider the S wave case, an approx- 


imate S(k) will be obtained by taking 


(47) rio) > De i 


and for real energies 


(48) dira Ara 


1 + fest) da . 


lu 


The question is to know what is the error committed. 
When C(a)< M*", for real k for instance, 


M 1 
1 +f) 0,(%)|da< exp i =. 
n | 
u 
and 
a lo ] 
lost az Ju SÈ ne 
È ORNE ue n u'| 
(49) Nu Nu 
n MMI 1 
Jose ae a ica A o 
|. n (Nu) yan 
Nu 
r goes to zero as N goes to infinity. 
Therefore 
M 
(50) 0(k) = Arg Î _ fox) da| + sin! = : ; 
LN, | 
Ë 11 +] on(a) da 
lt 


So an explicit bound can be put on the error committed, 


(case 7 < 1), 


As an application let us mention this small result if we take the Yukawa 
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potential: the condition to get a bound state at zero energy is 


replacing this by 
a 
V J a | 4 
Leo) == ie lee lon Ur 
; Ze OE u 3 3 
tt 
we get a well depth parameter 
1 
Bae 505% 


to be compared with the exact result 0.5953. 

For 1-0 or for tensor forces the practical use of the method is not so 
obvious. However, the method we use seems to be very suitable for a computer 
n view of its iterative character. 


32. The deuteron wave function and the anomalous threshold problem (577). — 

We have already said that our method furnishes an integral representation 
of the ingoing and outgoing wave functions. In the special case of a bound 
state the wave function f(— ix, r) exp [— xr] is the regular wave function 
because /(— ix, 0) = 0. Hence the deuteron wave function may be written 
(not the reduced one): 


(ce) 


(51) pn) = Ng = =H a3 | Qu(x) exp [— ar] de 


i 


This is the most general form of the deuteron wave function for Yukawa-like 
potentials. In a sense it is a generalization of the Hulthén wave function 


_exp [— xr] — exp [— Br] 
Vu then = N°- . 


The analytic properties of any non-relativistic matrix element involving 


(11) L. BerroccHi, C. Crorrn and M. Tonin: Nuovo Cimento, 18, 770 (1960). 
(12) R. BLANCKENBECKLER and Y. NamBu: Nuovo Cimento, 18, 525 (1960) 

(13) A. Martin: Nuovo Cimento, 19, 344 (1961). 

(14) A. MartIN and R. Vinn Mau: Nuovo Cimento, 20, 246 (1961). 
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a deuteron in its initial or final state can be investigated by inserting in it. 
the deuteron wave function. 
This has been applied to two cases: 


i) The non-relativistic deuteron electromagnetic form factor 


(52) F(q?) = fra exp 


Typical contributions to this will be obtained by replacing |wda(r) |? by 


| ui TRA a ee | LEE 
r pe 


exp [— ar] exp [— fr]/r? (%, 6 >x +4). 
Let us call 


pat 


Nes 
2 


Spies dr 
ra î 


peo [exp PE web ria de dr, 
da : 2) r 
1 
3 ae const . 
(a + 5)? + (q/2) 
Hence 
j dy 


ge By q°) = const | (y +B)? + (gR2)' 


x 


therefore in the (q?) complex plane g(a, B, (q?)) has a branch point at C= 
= — 4(x + B}?. 

From this we conclude that the analytic properties of F(q?) are as follows: 

it has a branch point at g2?=—16#2 = 

——4MB (where B is the deuteron 

F(q?) binding energy) and then a cut from 

PA} 4 fb) tO ¢* === co ee Fig. 3). 

It can be noticed that the discon- 

-4 (22-11)? “1622 tinuity across the cut Starting from 

q? = —16%? (which is placed for con- 

venience in the same direction as the 

free next cut) can be calculated exactly in 
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the region 


— 16%? + 


4(2x4 u)?, 


from the knowledge of x and N alone, i.e. of the asymptotic form of the deu- 
teron wave function. 

If the N'—p potential is partly known or, more precisely if C(x) is known 
in the range 


peo N IF 


(see next section) one can use our equations to calculate o_,,(%) and then the 
discontinuity of F(g?) across the cut till gq? =—4(2x+(N + 1)u)?. 

It is interesting to compare the result of this 
analysis with the field-theoretical treatment. In this 
treatment one considers the deuteron as an ele- 
mentary particle coupled to two nucleons. In per- 
turbation theory the simplest process of interac- 


tion of a photon with the deuteron is given by 
aio ai 

Let me remind you that the standard procedure to look for singularities 
with respect to q?, the 4-momentum of the photon is to look at the lowest inter- 
mediate state in y->d-+d. The lowest intermediate state lies at @?=— (4m)? 
one would believe that between q?= 0 and q? = — (44?) there is no singularity. 
However, if one calculates the singularities by explicit calculation of the 
Feynman denominator, as was done by KARPLUS, SOMMERFIELD and WICH- 
MANN (15), it turns out that the first branch point is much closer to the physical 
region. It is called an anomalous threshold. With very good accuracy it lies 
at gq? = —16x?, which is just (setting qj=9, which is the case if we consider 
the form factor as a part of e-d scattering amplitude in c.m. system) the 
position at which we found the first branch point in the non-relativistic treat- 
ment. Therefore the appearance of anomalous thresholds is directly related 
to the weak binding or the large extension of the deuteron. If the deuteron 
binding energy were much bigger, the relativistic treatment predicts that the 
anomalous thresholds should disappear. The non-relativistic treatment is 
unable to give this result, and this for the simple reason that in a tightly 
bounded deuteron the non-relativistic quantum mechanics would no longer 
be valid. 


ii) A similar but much more complicated case is the deuteron photo- 
disintegration matrix element. We can write this matrix element, neglecting 


(15) R. KaRPLUS, C. SOMMERFIELD and M. WICHMANN: Phys. Rev., 111, 1187(1958); 
Phys. Rev., 114, 376 (1959). 
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spin complications as 


1 a) a(x) der , 


ne Pi OX 

(93) | vei) exp 9 
where p(x) is the Np wave function (incoming wave+ plane wave) with 
em. momentum p;, and q is the photon momentum. One has to study the 


analytic properties of 


pPi (x) E n aS) exp i a) de 


i 


where «=x, or «>x,+u. This object is very similar to a scattering matrix. 


element 


den 


fe exp [— aa] a 


yet) exp ipixl 


Le 


As far as the analytic properties with respect to the angle p,q for fixed 
\p;| and |g| are concerned, both matrix elements can be treated in the same 
way. One can just copy the results for the analytic properties of an off-energy 
shell (| p,|+|p,|) scattering matrix element to get the analytic properties with 
respect to the angle. 

We shall speak later about the fixed energy analytic properties. So we 
shall just give here the result for the deuteron photodisintegration matrix 
element: 


it has a pole at (ps — q/2)} =— #, 
and a cut starting at (p,— q/2)? = —(x + yes 


The fixed transfer dispersion relations for this matrix element are not 
very easily established because of the complications due to kinematics. This. 
however, can be done (1) and by combination with the fixed energy dispersion 
relations for nearly real energies one gets a Mandelstam representation for 
this matrix element. 

A much simpler thing to do is to look at the analytic properties of the 
partial waves for this matrix element (8). Then one can also replace the final 
wave function by its integral representation derived in Section 2. Needless 
to say that the analytic properties obtained in this way coincide with those 
obtained by projecting the Mandelstam representation. 

The field theoretical analogue of this matrix element is described by the 


co 
x 


ti 
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following graphs Fig. 5. The second order graphs correspond to the . poles 
appearing in the energy or in the momentum transfer variables. The fourth 
order graph gives rise to anomalous 


thresholds in agreement with the non- N 

relativistic treatment (cut at (p,—q/2)}— o à 0 

= —(x+u)). d d N 
It is interesting to notice that in DL 

spite of the presence of anomalous 

thresholds the Mandelstam representa- 

tion holds, at least in the non-relati- 

vistic approach. Y 5 
Other applications can be made d oP bn 

using this technique, for instance to N 


study the electromagnetic form factor of Etes, 

the «-particle (1%), which can be consi- 

dered either as a d-d system or a .N°-*He system. Then one gets cuts due to the 
extension of the x-particle but one can show that the closest cut starts relatively 
far from the physical region (q? = 114?) 
and that it is possible to observe the 


x sy 
Z\ of à : : : : 
105 effect of light intermediate particles, for 
PE 
a x x a instance a T=0, J=1 particle of 
mass 2.2u,, which gives rise to closer 
Fig. 6. singularities (Fig. 6). 


33. The connection between the unphysical cut of a partial wave amplitude, 
the S matrix, and the potential (17). — In recent years many theoreticians 
have attempted to relate the nucleon-nucleon scattering matrix to pion-pion 
and pion-nucleon scattering, using Mandelstam representation and unitarity. 
In many approaches, the primary object which is calculated or estimated is 
the discontinuity of the S matrix of a given wave on the unphysical cut in 
the energy complex plane. Once this discontinuity is known one has to construct 
a scattering amplitude such that it has the correct discontinuity on the left 
and fulfils the unitarity condition on the physical cut (*). 

In the simple case of potential scattering (for Yukawa-like potentials) 
one would like to be able to construct the S scattering amplitude (or any 
other partial wave amplitude) from the information on the discontinuity of 
the left-hand cut in the k? (i.e. kinetic energy) plane (Fig. 7). This might prove 


(6) L. Berroccut, 8. Fusivi and M. Tony: private communication, 
(17) A. MARTIN: Nuovo Cimento, 49, 1257 (1961); V. DE Auraro and T. REGGE: 


Nuovo Cimento, 20, 956 (1961). 
(18) H. P. Noyes and D. Y. Woxa: Phys. Rev. Lett., 3, 191 (1959). 
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useful in the more general case, because, if we succeed in doing so, we may 
forget the intermediate assumption of the existence of a potential, and obtain 
an acceptable solution, as long as relativistic 
kinematical effects are unimportant. 
E complex plane ; e 
physical cut In our study of analytic properties of 


+ |} partial waves we have used k rather than k? 


R | Ara Le NE threshold as a variable: the upper half % plane corre- 

sponds to the first k? Riemann sheet (with 
va biso: respect to the physical cut). The cut 
i k =iu/2 — too is precisely the unphysical cut. 
“4 The lower cut of the Æ plane is hidden in the second Riemann sheet of 


the energy complex plane. 
So, what is supposed to be known is the discontinuity S(ix +e) — S(in— e) 
across the upper cut in the k plane: 


à) TARE (54) Six + ©) — Six — e) = D(x)'2ri, 


wee... Where D is a real function. 
: We shall write 


È n as was previously done and remember that for potentials such that O(a) <o1-” 
ie elites Omer 
100 lim DALE — < const. 
M, > |k|—> [x |" 
i This enables us to write an unsubtracted integral representation for f: 
È LE (R(x) al. 
à (55) Hb 0) = il DI 
à « J k—ix 
vi u/2 
or R(x) is a real function because f*(— k*, 0) = f(k, 0). We get 
ae — i/ R( x) dx] (k — ix) 
SE 56 rr es wile 
; (56) SE) = —— 2. 


1 +i [Re ) dæ/(k + ix) 


nl 
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it is a trivial matter to evaluate the discontinuity of S across the upper cut. 


(87) E AIA Fe 
1 sii dæ'}(æ' + x) 
ua 
or 
(58) R(x) = — D(x) ! + free) dx'/(2'+ 2)|. 


/ 
H/2 


This equation can be considered as a Fredholm (or almost Fredholm) type 
integral equation determining (x), and hence S(k) from the information on 
the discontinuity. This equation was obtained by NovEs and Wong (18) for 
the nucleon-nucleon case, at least for the non-relativistic kinematics, in a rather 
involved way, as the result of the combination of two coupled integral equations, 
one expressing the fact that the left-hand cut has the correct discontinuity, 
the other expressing that the physical cut discontinuity (in the energy variable) 
satisfies unitarity. 

Our method can be actually generalized to the general case where no po- 
tential is assumed. This has been done by Omnes (!°). One still writes 

Se aioe ; 
Î(— k) 


7k), as in the potential scattering case has no singularity in the lower half 
k plane, and goes to unity as |k|— oo. A natural definition for i (k) will be 


foe} 


1] [dar 
AT pay? 


0 


(for k in the lower half plane), where v is the kinetic energy. When the Levinson 
theorem is fulfilled and when no bound state is present (this imposes 
6(0) — 6(c0) = 0) F(k) coincides with our f(k, 0). f(k) can be prolongated in 
the upper half plane by 


(*) In ref. (17) the sign in incorrect. 
(1) R. OMNES: Nuovo Cimento 21, 524 (1961). 


12 - Supplemento al Nuovo Cimento. 


771 


182 A. MARTIN 


Then one shows that f(k) satisfies the same kind of integral representation as 
the one we proposed. The equation one ends up with is similar to ours but 
a little bit more complicated because: 


i) of relativistic kinematics; 


ii) the bound states no longer appear as zeros of ae k) for k in the 
upper half plane but as poles of f(k); 


iii) if 6(0) —d(co) #0, but rather na, f (k) has a multiple pole at k=0. 


Let us come back to our original equation. Its solutions have been studied 
by REGGE and DE ALFARO (17) who have shown that if 


the solution is unique; also, if D(x) is positive and A finite the solution exists. 
Further they show that the number of resonances -antibound states is less 
than two plus the number of zeros of D(x), if the above condition is fulfilled. 

The next question we shall consider now is the relation between the po- 
tential and the discontinuity D(x). One might question about the usefulness 
of constructing the potential since we have seen that from D(x) one can 
directly construct the S matrix. However, three points are worth noticing: 


i) when the discontinuity is not completely known S cannot be con- 
structed, but it will turn out that one can construct explicitly the external 
part of the potential, 


ii) information on the external part of the potential permits to deter- 
mine the external (however, inside the interaction region) part of the deuteron 
wave function which as it appeared in b) is so important for the deuteron 
form factor; 


ili) nuclear physicists have to deal with potentials in the many-body 
problem. 


We shall first calculate the discontinuity from the potential. In order to do 
this we have to look more carefully at what happens for k = ix +e (x real> u/2). 
f(k, 7) has no longer a unique meaning but we have some hope to define 


(59) ft (ia, r) = lim f(ix + e, r) , 


and, similarly f (ix, r). Correspondingly we define 0;.(a) and 0,(a). Now we 
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notice that from the equation 


(60) Ha 22 + ie) (a) = C(x) + [Ce — Hot) aB 


i) o(&) —0:(x) — 0 for a<2x; 


li) near x=2%; 


dl 


O(a) +fO(a — B)0(8) ap 


x 


(61) OX) — out) = — 2rid(a — 22) 


There is no ambiguity in the right-hand side because of i). 
ili) in the range 2a+e<a<2x+u—e, 07,(a) — p;,(x) = 0. 


Now there are several possible ways of calculating the discontinuity of S(k) 


) — Seite) — cad 
1 a a) da 


(62) S+(ia) — S-(ix) = 2niD (x 


The most straightforward may consist in noticing that ft(ix, 7) exp [ar], 
f-(ix, r) exp [ar], f(— tx, r)exp[— ar] are not independent solutions of the 
Schrôdinger equation. The only possible relation between these quantities, if 
we take into account the fact that 


fr(ix, co) =f (ix, co) = f(— ta, cc) = 1, 


ft (ia, r) — f-(ta, r) = 2in A(x) exp [— 2ar](— ta, r) . 


If we look at this relation for r —0 we see that A(x) is identical with D(a). 
So we can evaluate D(x) by taking the most suitable value of r; in particular 


(63) 2riD(x) = Jim exp [2er][ft(îx, r) — f-(a, r)] 
2æ-+e 
Op 


+ lim exp [— 2xr] | [Oie(%) — 0;2(2)] exp [— ar] da. 


20+u—e 
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The last term does not contribute to D(x) because it contains a factor 


exp [— (u— e)r]. The contribution of the first term is trivially evaluated. 


It gives just 
2a 
C(20) + [Oa — A)0;r(*) da 


ù 


(64) Da de 
In this formula 0;,(«) is a well defined quantity because the range of inte- 
gration is such that the condition x< 2x is always fulfilled. At this point we 
have really proved that the lines 


7 , ai : 
k=il >i00 and. be i 


a 4 


are cuts of the function exp [2id(x)]. We have an explicit expression of the 
discontinuity. One can notice that if the potential is written 


œo 


Vir) = Jew exp [— ar]da, 


H 


the discontinuity in the range 4=—n(u/2), x —(n+1)(u/2) is given by an n-th order 
polynomial in 4. Therefore it is just S(ix +e) — S(ix —e) calculated till the 
n' Born approximation. However, the expression we have derived is exact, 


even when the Born series fails to converge. 
It is interesting to notice that an attractive potential with C(«)<0 gives 
rise to a purely positive D(x). Indeed looking at 


x 


L(x = 22) ix (2) As C(x) + fota ai B) Cin (B) )dB, 


Lu 
we easily see that C<0 makes 0,,(&«) > 0 for a<2x. Inserting this in eq. (64) 
we see that D(x) has to be positive. Returning to eq. (58) we see that if there 


is no bound state with binding energy larger than B= (u?/4)/M (which is the 
case of the Np system) 


“er L. ; 


is positive for «>/2, and hence R(x), if D(x) is positive, is purely negative, 


7 
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Then 


R(x’) da’ 


can only vanish once in the range æ—0, x=y/2 and there is at most one 
bound state. Hence, if a system interacting by a purely attractive potential 
has two bound states one of these bound states must have a binding energy 
larger than (u?/4)/M. Let us also notice that if R(x) is strictly negative there 
is no zero of f(k, 0) outside the imaginary axis because then f(k, 0) is a Herglotz 
function (Im f has the sign of Re k). This means that there is no resonance 
since there is no pole close to the real axis in the lower complex % half plane 
There may be, however, a pole on the imaginary axis between k=0 and 
k = —i(u/2) which would correspond to a «virtual state » (as in p-p system), 
but then the true bound state cannot exist. 

Let us now exploit formula (64) in order to construct the potential from 
the information on the discontinuity. 

If we know D(x) in the range «= p/2 > «= 2(u/2) we get 


C(2x) = — 2a D(a). 


Let us go now to the range =p, æ = 3u/2 we know from the preceding 
result that for «<2u 


O(2x) = — 2aD(x) + fre — 8/2) D(B/2) dB. 


u 


So, here again we are able to construct O(a) from D. 
The procedure can go on: if we know O(«) in the range «=p, a = nu 


and D(x) in a larger range, say 
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We know 0;.(«) for a<nu and we can calculate the integral appearing in 
the r.h.s. of (64); the 1.h.s. is known by hypothesis; hence C(2x) can be deter- 
mined for nu<2x<(n+1)w. 
In general, when the discontinuity is known for 2x< Nw the contribution 
to the potential 
Nu 


Jew exp [— ar]da, 


is known. 

We see that the further the discontinuity is known, the better the inner 
part of the potential is known. If the discontinuity is completely known and 
if the conditions imposed by Regge and De Alfaro are fulfilled, it seems that 
the complete potential 


foe} 


few exp [— ar]da, 


u 


is acceptable because we know anyway from previous work by GEL'FAND, 
LEVITAN and others that one can fit the phase shift by a potential. 

The condition that the maximum range of the potential is 1/u is a very 
important one. It permits, as long as there are no bound states with «> u/2 
to distinguish between the zeros of f(— %, 0) and the singularities of f(k, 0). 
This is the reason why in our treatment we do not get various phase equiv- 
alent potentials. 

In case a hard core is present, the treatment as given here breaks down 
because the scattering amplitude has an essential singularity at infinity. 
However, if the range of the hard core is known the procedure can be adapted 
to this case. 

One starts from 


V(r)= + co for Pot 
V(r) =|C(a)exp[- a(r— »)]da for rire, 
u 
one defines 
S'(k) = S(k) exp [2ikr;], 
and 
D'(x)= D(a) exp [— 2xr,], 


from D'(x) one can construct O(«) in the ordinary way. 
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The connection between the discontinuity and the potential can be extended 
to an arbitrary partial wave. Formula (64) can be written as 


(65) D(a)=— lim 0;;(a)(a— 22). 


>2x 


Then the expression of D(x) can be obtained in terms of C(«) and g(x) which 
itself can be calculated from C(x). When one tries to go the reverse way, 
one finds that at each step of the calculation one must solve a Volterra equation 
to construct 0 from D. In this latter case it is clear that D must satisfy 
stringent conditions to ensure that the phase shift has a correct low energy 
behaviour. These conditions are not apparent in our treatment because they 
can always be satisfied by a convenient choice of D(x) for large +. 

The treatment can as well be extended to the construction of potentials 
in presence of tensor forces. Starting from the discontinuities Doo Die Dez of 
Soo(k) Soo(k) So(k) one can construct the inverse Laplace transform of Vo, 
Van iV. 


tensor * 


4. — Analytic properties for fixed energy (*). 


The complete scattering amplitude is a function of two variables f(H, t), 
where t = —2k?(1— cos 0), k is the c.m. momentum and 0 the scattering angle, 
and H=k? is the energy. 

The most straightforward approach to the analytic properties for fixed 
physical Æ with respect to ¢ is to make use of the Born expansion. However, 
it will turn out eventually that the analytic properties of the exact amplitude 
can be obtained from those derived by using the Born expansion in a very 
easy way. The only question this approach cannot solve is the behaviour of 
j(E,t) for too. Among all the approaches to this problem only the one 
by Regge furnishes the solution. However, this does not mean that the present 
approach is useless, because it gives another insight of the structure of the 
discontinuities. 

For convenience we assume for the potential 


V(r) = | Ta) HSE da, 


be 


i.e. we take it as a superposition of Yukawa potentials rather than expo- 


nential potentials. 


ti 
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The n-th Born term can be written, dropping numerical factors: 


(67) <il VORRAI with n V’s, 
and 
Gor = i 2 a 7 "I , in coordinate space, 
pe 
(05) kia — ka — 8 : 
= ; in momentum space. 
il 


Inserting the representation for V(r) we get a typical contribution to 
the n-th Born term: 


1 1 il iL 
a? + (k'— ky)? ki — k® — te 03 + (hy — ka)? ah + (kn ke)?” 


(69) far ae … 


with %,%,..-,%,>p this formula should be integrated again on @%,..., &n, 
with the weight function J'(x). 

To study the analytic properties with respect to the angle (k', k) which is 
equivalent to the study of the analyticity in { we just need to make the 
angular integrals. Let us consider first the integration over ky. 


foo 1 1 = [22 1 Livre ins 
og d-(k'— fey)? Of + (hey ka)? J4bkkiu,— (lE) Ua (ky bea) 


with 


HER DL ot ta 
sa 


U, = 
i 2k'k, à Qik,” 


and the ~ mean unit vectors. 
This integral by means of Feynman parametrization can be transformed into 


sige 


et U, | 9 (kK, ka) 


2 2 ; 


dl 


enr | dy 
4k'k,k? AE 2) OPEN Ne 


or 


cia 


dd 


27 | dy 
ar ele) È a; “| e ((1 + y)/2) + U, ((1— y)/2)P — (14+ y?)/2 
2 22 
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Let us rewrite the bracket as U,,— kb, and minimize it first with respect 
to |k,| (since k, goes from 0 to co) 


al œ + k'2 Oke (LL Ey) (= y) Lg 
1 | 2 2 1 È È 
12 È 
(1 —y2)/2 


Now since (4/X)+BX=>2VAB 


E (gl 7 
E Ex y) ke ane 
1 — y? 


1 2 a1 2 2 2 / = 
= G Ly) a; +k X + =F ry n de di i (1 + 92) 


Letting t= (1 + y)/(1—y) this becomes 


1 fe +, | % + ha) [1 a) 17 
k' AN RETIRE" 


which, again can be minimized with respect to t. Eventually 


(x, + a)? + hk? + k3 


di 2kk, 
Therefore 
[iene ier tee) a + ke 
122 2kk; , 


with A4,,>0; A, is a function of y and k,. U,, has the same form as U, and U,. 
We can now perform the integration on d2,. We have to integrate 


1 1 
USS a R'È, US = 7 DIA 


; 


obviously the same reasoning can be applied. Eventually the n-th Born term 
will be written as 


? 


[ET variables) d(lengths of k’s) x Green’s functions x 1 
products of ((1— y2)/2)X[Utss..n — È] 
with 


(0, + Oa vs + an)? + Bi? + Re ney? + he + 
Dkk’ 2kk" i 


UT, n > 


i 
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The remaining integrations give no trouble because they do not give any trouble 
in the physical range (*). 


The structure of the denominator exhibits a cut in the cos 6 plane from 


nur + k'2 + k? 


(70) Cos O = Eh! to +00, 
or in ¢=— (k'—k)?, from t= nu? to +00. 
Notice that we do not assume |k|=|k'|. Notice also that the result 


holds if the k2 in the free Green’s function is different from % and k’. Let us 
now consider the full scattering amplitude 


(71) DUR E RE LA, 
where 
(72) I, =r — TI, k) GTR, k, k) deh" , 


in coordinate space 
(7 — T,) (k', k’, k) - [exp [— k'x] V(x) G(xx') V(x') exp [ik"a']d30d?x', 


where G, is the exact (not the free) Green’s function. G,(a, x’) is regular and 
bounded by const/(|#’—«#|) in the upper complex half plane in k. This still 
holds for real k. If | V(x)|< Ctexp[— ux] which is the case here, the con- 
vergence of this integral is dependent on the factor 

exp [— th'x'— u|x'| + ik’x— |x|]. 
. Let us take 
k"= U cos 20 + V sin 20, 
k’= U cos (1 — 2)0 + V sin (A—1)0, 


with UV orthogonal unit vectors. Then the real part of the argument of the 
exponential is at most 


| |[k" V (Im cos 20)? + (Im sin 40): — u] + 


+ |e' |[k' V[m cos (1 — 2)6}? + [Im sin (1 — 2)0]}? — u], 


(*) The singularities at y—+1 may be shown to be harmless. 
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cos 70 = cos (0, +10.) = cos 70, cosh 70, — à sin 70, sinh 20,; to make the first 
bracket negative we get 


k" sinh 20, —u<0, 
and the second 


k' sinh (1 — 2)0, — u <0. 


Therefore the integral exists in the smaller of the ellipses in the cos 0 plane 
with foci —1, +1 and semi major axis 


IE 
cosh 6, = cosh 3 sinh- Li? 


1 : IZ 
108h 0, = cosh ——— sinh —. 
iin A PIE) k' 


The maximum domain is obtained when these ellipses coincide, by convenient 
choice of 2. 
Then 
x IZ 3 u 
Se =1 1 4 —1 ! 
0, = sinh pt sinh pr? 


and the semi major axis is 


| a 2 Ue 
a= |/1 na IR i a kik"? 


This means that T(k'k",k) can be expanded in Legendre polynomials 


Tk, k’, k) = > APAE, &"), 
with 


1 
lim su TA re 
eed à Ca+Va—1 


Now 7,(k",k,k) can also be expanded in Legendre polynomials because from 
its analytic structure this expansion is certainly convergent in an ellipse with 
foci —1, +1 and semi major axis 


np + e + ke 
RE) 
Therefore 

dhe k, k) a > BP i(k"; k) ’ 
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with 


Il 
lim sup (B,)!" = —- , 
nr e 


One can insert these expressions in the expression of k,. The angular inte- 


grations can be performed because 


È ES SA 3 AT 
[ra dQ = dw TTI 


Hence 


k? dk’ AG Ay By SIA 
= —— P,(k'k). 
RA fa) 21 21 Le 


The expansion under the integral is convergent inside an ellipse with 
foci —1, +1 whose semi major axis c is just given by 


1 Il al 
——— = lim sup (A,B,)"* = — — 
AI) <& DV a The Vip: ia 


or 


C= Abe V ay bete 


Now 


e = cosh |sinh=! = + sinh - + cosh 
y Ù 


NE ee 
2kk" i 


We have to minimize it with respect to k”. Since cosh and sinh are mono- 
tonous we can minimize first 


sinh |sinh-! 12 a Coni np? + Je2 sa ki 
k" Ohh! a 
If 
2u? 2 112 
e REA 
2kk 


This quantity is XY + V¥2+1V ye—1 but this is just the minimum with 


respect to A>0 of 
il Yi 1 
sl + a 
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We minimize first this with respect to hk”: it gives 


1 k2 + mu? 1 A il 1|2u | nu] (n +1) 
io Sea le ee 


x 2 a2 k 
Hence 
Vira 2 / TA 21 
(73) Croinimnm = Vi fal Vi (+? (n+ le 
minimum = | k'2 | ! k2 | kk’ | 


The region of regularity of R, is an ellipse with foci —1, +1 and semi- 
major axis c increasing with n. One can make this ellipse arbitrarily large. 
Inside this ellipse the only singularities are those given by the Born expansion, 
i.e. a series of cuts starting at 


k? + k'? + u? , Na? 
2kk' a 2kk' 


cos 0 = 
(ORI 


We have now proven that outside the cuts there is no singularity. The 
only question we cannot answer is the problem of the behaviour at infinity. 
As we increase the dimensions of the ellipse we must take into account more 
and more Born terms and it looks impossible to obtain in this way an idea 
on the asymptotic behaviour. 

At a given point of the cut the discontinuity can be in principle calculated 
since it consists of a finite number of contributions. As was done in the case 
of the analytic properties of the partial wave, one can here also try to con- 
struct I'(«), the weight function appearing in the potential from the discon- 
tinuity in the t plane (*). The reasoning is very similar to the one we have 
already done: in t=? t—4u? the only contribution to the discontinuity 


comes from the first Born term 


The discontinuity is easily seen to be proportional to 1(v t)//t. Therefore 

I'(«) can be constructed in the range a= y, x= 2u because 7, and the remaining 

terms have singularities starting at least at t=4u?. Let us now generalize this. 
Let us write the potential as 


V= Vee + ere a= Viana ’ 


(2) A. Martin and Gy. TARGONSKI: to be published in Nuovo Cimento. 
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where 
nu (n+1) pi [ I 
i cp [— exp [— ar 
Wa = |P) SMI da ’ Per = [rc Si ue da, 
i nu 


Vase | ra) 


(n-1) 


exp [27] 4, 
Yr 


Assume now that we have been able to construct V,,. We want to get Vin 
from the discontinuity of 7 in the range t= n?u? t=(n+1)?u?. In this range 
the discontinuity is given by 


Now the terms in which V,:,, appear do not contribute to the discon- 
tinuity because they give rise to a cut starting at least at t= (n+1)?u?. The 
contributions 7, +... + 7, in which V,,4, appears not to contribute to the 
discontinuity either, because they are combined at least with one potential 
of range 1/u and the singularities start at least at t= (u+nu)?= (wo 1) 202) 
Therefore the discontinuity AT is given by 


AT = AT HAT, ++ L)(Vaoa) 


? 


where 7,... T, are calculated from V,, alone which is supposed to be known. 
This equation determines AT, in the range t= nu, t= (n+1)?u? and, because 
of the special form of 7, it determines /«) in the range «= nu, «= (n+1)u. 
One can then iterate. Eventually one can construct I'(«) from «=p to «= Nu 
if AT is known in the range t=u? ->t=N?u?. Here again the farther we 
know the discontinuity the better we know the potential. 

This procedure furnishes an unambiguous definition of an energy-dependent 
potential fitting at a given energy a scattering amplitude having the prop- 
erties implied by Mandelstam representation. In the low energy region, for 
the nucleon-nucleon case, FUBINI, CHARAP and TAUSNER (21) have shown that 
at least the external part of the potential (one and two-pion exchange) is energy- 
independent. 

A possible application of this is the study of the range of the imaginary 
potential in nucleon-antinucleon scattering (°°). For this one makes the pre- 


(21) J. CHARAP and S. FUBINI: Nuovo Cimento, 14, 540 (1949); 15, 73 (1960); 
J. CHARAP and J. TAUSNER: Nuovo Cimento, 18, 316 (1960). 


(22) A. MARTIN: CERN preprint. To be published in Phys. Rev. 
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liminary remark that a purely elastic scattering amplitude (no other channel 
open) will be automatically associated with a purely real potential. In this 
case the unitarity condition reads 


Ju) Im [V(x)lyvz(x) dba = p(t) = 0, 


for |k|=|k’|. q(t) can be easily seen to have the same analytic properties 
as the scattering amplitude in t; since it has to vanish in the physical range 
it vanishes everywhere and one con- 
cludes from this that Im V(a)=0 = 


(assuming for Im V(#) the same re- dI N A di 
presentation as for V(2)). a 

We can obviously apply our for- ae PT 
malism to NN-NN, but then à 
the potential we get is clearly com- /, K Ke ae 
plex. Two kinds of graphs contribute Ta TM 
to NN NN. T,is purely elastic Fig. 8. 


below the meson production thresh- 

old, i.e. if we use an energy such that the reaction NN > NN +7 is impos- 
sible. 7, has the following analytic structure in the cos? plane (see Fig. 8 
and 9). 


two pion cut 
ne 


(nue 
cos@ : ar 
4 ni + + per rr nee 
to =I +1 | 
NN cut one pion pole 
LL 
4M = 
-(1+ ) ake 
2k? a 
così 
+ + uv 
=(14 44) vd #i +(1+ 4047) 
2k? Ts 2k 
mig Ob 


The singularities of 7, are symmetrical. 
With such analytic properties we have to use a potential 


i (0) a da + P| Te) SIR si da, 


Ù 


where P is a Majorana exchange operator. There is no difficulty to extend 
our treatment to this case: to construct I°, and J’, one has to progress sym- 
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metrically with respect to cos 6. The important point is that obviously /.(a) =0 
for «<2M, and /,(a) is given by the discontinuity of T, alone for a<2M. 
But since 7, is a purely elastic amplitude, from the above argument /'(x) 
has to be purely real in the range «=m, «=2M. Hence the range (defined 
by the behaviour exp [— 2Mr]) of the imaginary potential cannot exceed 1/2M. 

One could provide many other applications of these fixed energy properties, 
in particular the study of the deuteron photodisintegration matrix element 
which has already been mentioned. Some people are trying to make an explicit 
calculation of the discontinuity for a Yukawa potential to see from what point 
the oscillations predicted by Regge start. 
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1. — Introduction. 


Present-day quantum electrodynamics attempts to unify the dynamical 
properties of a charged particle system with quantum mechanics. The former 
is assumed to follow, in strict accord with Einstein’s principle of relativity, 
the Maxwell-Lorentz field theory. The latter depends on the association of 
a probability amplitude with the trajectory of each of the constituent charged 
particles. The probabilistic interpretation is inherent in the measuring process, 


(*) A preliminary account of this paper was given at the Berkeley and New York 
meetings of the American Physical Society, Bull. Amer. Phys. Soc., 5, 505 (1960); 
6, 60 (1961). 

(*) Part of this work was carried out while the authors were at the Department 
of Physical Science, San Jose State College, San Jose, California. 
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as expounded by Bohr’s complementarity principle. The theory is therefore 
non-deterministic. 

The apparent incompatibility between Hinstein’s deterministic, completely 
field theoretical approach and the non-deterministic approach of quantum me- 
chanics has been the subject of much debate ever since the initial successes 
of quantum mechanics (1). Although this incompatibility is unimportant as 
far as the mathematical description of non-relativistic phenomena is concerned, 
it is indeed all-important in the exact (relativistic) description of quantum 
electrodynamics. 

The deterministic approach of EINSTEIN has not been taken seriously by 
the majority of physicists for the two important reasons that non-relativistic 
quantum mechanics has bad a great amount of success in its predictions and 
also because of the lack of an alternative deterministic theory which would 
duplicate the predictions of the properties of atomic and other elementary 


particle systems. 


In order to explain high energy phenomena, it becomes necessary to extend 
the non-relativistic theory of quantum mechanics into the relativistic domain 
if relativistic covariance is to be taken as a principle underlying the description 
of all physical phenomena. The development of present-day quantum electro- 
dynamics represents the construction of a field theory which attempts to in- 
corporate the non-deterministic approach of quantum mechanics with Ein- 
stein’s principle of relativity. A compilation of many of the important papers, 
representing the development of present-day quantum electrodynamics is given 
by SCHWINGER (?). Also, FoLDY (*) discusses the compatibility of relativistic 
covariance with the standard quantum-mechanical theory. 

To display some of the successes of the present-day theory, one need but 
mention the calculations of the Lamb shift and the anomalous magnetic mo- 
ment of the electron (*). On the other hand, the theory does suffer from some 
undesirable features. First, there are divergences which appear and must be 
removed by additional physical interpretation and further mathematical tech- 
niques in order to predict the measured properties of the system. Second, 
the conventional formalism does not give results in closed form. All of its 
predictions are inherently contained in expansions which represent the per- 
turbing influence of the interactions within the system of particles on free- 
particle fields. Coupled with these objections is the question of the consistency 
of the present formulation of quantum electrodynamics. 


(1) See, for example, the article by N. Bonr and the reply by A. EINSTEIN in: 
Albert Einstein Philosopher-Scientist (Evanston, Il., 1949), pp. 201, 666. 

(©) J. SCHWINGER: Quantum Electrodynamics (New York, 1958). 

(*) di. Lis Fotpy: Phys. Rev., 422,-275 (1961). 

(4) For a survey of these calculations see A. PETERMANN: Fortschr. Phys 20 
505 (1958). | 
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These difficulties in the present-day theory have their origin either in the 
mathematical formulation or in the conceptual framework. A great deal of 
work has been done on the exploration of the mathematical formulation in 
attempts to obtain a more consistent theory (5). It is our purpose in this work 
to present the results of a study that took as its starting point a reinvestigation 
of the conceptual framework of quantum electrodynamics. The approach that 
we take, similar to that of EINSTEIN, attempts to describe natural phenomena 
in terms of a completely covariant, non-singular and deterministic field theory. 

An implication of the deterministic nature of the theory is that the solutions 
of the field equations provide a complete description of observable phenomena. 
This requires that the theory must necessarily be self-consistent. By self- 
consistency we mean that the system is closed (i.e., the measuring apparatus 
as well as that which is observed is included in the description). 

In the section which follows, the theoretical approach to our theory will 
be discussed, with emphasis given to our reinterpretation of the Dirac and 
Maxwell field variables. Section 3 then presents the Maxwell field equations 
in the two-component form developed by the authors (9). In Section 4 we 
develop a self-consistent fleld theory from a Lagrangian formalism. The Lag- 
rangian density is expressed in terms of the bispinor Dirac field variables 
and the spinor variables of the Maxwell field. From the variational procedure 
we recover the Maxwell equations (in their spinor form) and a set of coupled 
non-linear Dirac equations. 

The remaining sections in the paper apply the field equations to special 
cases. In Section 5 a result is derived for equivalent particles that is inter- 
pretable in terms of the Pauli exclusion principle. Section 6 is then devoted 
to the determination of the properties associated with a particular state of 
positronium. An exact solution is obtained for the coupled non-linear Dirac 
equations in the special case when the state of motion of the electron and 
positron are the same. It is found that this bound state of the particle-anti- 
particle pair exhibits all of the conserved properties that are normally asso- 
ciated with the annihilation of a pair and their simultaneous creation of a pair 

(5) One such approach of current interest attempts to provide a description of 
interacting particles without the use of «field » at all. Instead, physical events are 
described in terms of the scattering matrix that connects the state funetion of the 
particles in the remote past to the state function in the distant future, both of these 
times corresponding to the limit of no interaction. This phenomenological approach 
was originally introduced by HEISENBERG (Zeits. f. Phys., 120, 513, 673 (1943)) and 
most recently developed by LEHMANN, SYMANZIK and ZIMMERMANN (Nuovo Cimento, 
1, 205 (1955); 6, 319 (1957)) and others. For a recent compilation of papers on this 
approach see L. KLEIN: Dispersion Relations and the Abstract Approach to Field Theory 


(New York, 1961). ee | 
(9) S. L. Scuweper and M. Sacus: Bull. Amer. Phys. Soc., 5, 505 (1960). Also, 


manuseript in preparation. 
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of oppositely circularly polarized photons. Thus, in our theory the pair is 
not truly annihilated, but rather has a ground state with associated properties 
which make it appear so. It is this state of positronium (or any particle- 
antiparticle pair) which replaces the vacuum state in the conventional theory. 
It is also shown how the electromagnetic coupling to a third charged particle 
can give rise to the phenomenon interpreted as pair production. 

In Section 7 we show how the hydrogen spectrum (including the Lamb 
shift) emerges from our field equations. The Lamb shift A(m’S,—n’P,) in 
hydrogen is calculated for n'= 2, 3,4 and shows very good agreement with 
experiment and the conventional calculations, although no divergences appear 
at any stage in our calculation. 

Finally, Section 8 examines the limitations of the present work and dis- 
cusses its extension to other types of interactions. 


2. — Theoretical approach. 


The theoretical approach to the development of the covariant self-consistent 
field theory to be presented here is based on the following fundamental pos- 
stulate : 


The laws of nature must be described in terms of field variables which may 
be associated only with elementary interactions. 


It follows that the most elementary physical system is a two-body system 
because an interaction takes place between at least two physical particles. 
Thus we replace the absolute concept of the elementary particle by the operational 
concept of the elementary interaction. 

The dynamical behavior of the system is described in terms of the solu- 
tions of the quantum mechanical equations. The solutions of these equations 
are then related to the interaction of the constituent particles of the system 
at the space-time point x. Thus, x does not locate any individual particle; 
rather, if defines a point at which all of the particles interact. As a result of 
our requirement that the interaction rather than the free particle be considered 
elementary, the Dirac field variables, for the simplest system must be deter- 
mined from a set of two coupled non-linear differential equations rather than 
the single (linear) free particle equation. As the number of particles in the 
system increases, the number of coupled Dirae field equations must corres- 
pondingly increase. Since however, the solutions of each of the n equations 
are not independent of each other, we postulate that a single interaction field 
function ¥Y(y,..., y") exists, whose absolute Square is a measure of the 
weighting of the elementary interaction at the space-time point x. We thus 


/ 
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require a conservation law (*) 


D) 7 US == 

(2.1) TO 

to hold and interpret it as conservation of interaction. The normalization 
of Y, giving 


(2.2) [Fa ARE 


© 


is thereby interpreted to mean that the integrand weights the interaction 
between the constituent particles of the system throughout space. Such a 
result would then imply that as time progresses, particles could not be created 
from a vacuum, nor be annihilated with the subsequent creation of non- 
interacting particles. Such cases alter the value of (2.2). The first increases 
the total interaction while the second decreases it. Thus, our theory must 
predict the experimental results which are interpreted as pair creation and 
annihilation without actually creating or annihilating physical particles. This 
result is indeed obtained when we consider a particular state of the bound 
particle-antiparticle pair (Section 6). 

In essence, our interpretation of the solutions of the Dirac field equations 
replaces the description of an n-particle system, in terms of space-time co- 
ordinates tied to each particle, by a four-dimensional space-time with an asso- 
ciated interaction field. 

Once y(x) and y?(x) of the two particle system (or p(x), ..., y(x)) of 
the n-particle system) are determined, the conserved properties of the system 
may be obtained from the invariance properties of the Lagrangian density, 
in accordance with Noether’s theorem. Thus, to determine the constants of 
the system it is only necessary to have the solutions py, y? and does not 
require the form of the Ÿ field variable as a function of y? and y?. On the 
other hand, a knowledge of the functional form of the Y field gives us infor- 
mation pertaining to the actual weighting of the interaction throughout space- 


(*) The notation used is 


id 0 : 
Vente vy 5 Vp IA 
x (703 Vics Vo 0 1 


p= vr, vis = Vii k=1,2,5, 


and o, are the usual Pauli matrices. 
Also, we take h—c—1 throughout. 
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time. The derivation.of this field function, for the case of equivalent par- 
ticles in the same state of motion (Section 5), leads to the Pauli principle. 
It is also used (Section 6) to describe the dynamics of a particular state of 
a particle-antiparticle pair and to identify it with the standard interpreta- 
tion of the annihilation of a pair with the simultaneous creation of a pair of 
oppositely circularly polarized photons. 

Not only the interaction field 7, but also the individual Dirac functions py” 
satisfy the continuity condition (2.1). This is, in general, a consequence of 
the invariance of our Lagrangian (Section 4) with respect to gauge trans- 
formations of the first kind (7.e., a phase change in w®). We interpret 


2, y rp! di 2.0 
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within the framework of our theory, to mean that each of the sources of the 
total interaction is conserved. This, once again, supports our earlier statement 
that within our theory particles cannot truly be created or annihilated. In 
view of (2.3), we then interpret the individual y°° field variables to describe 
the kinematic behavior of the constituent interacting particles through the 
velocity field poy, yp. 

Thus far, the discussion has developed our interpretation of the Dirac field 
as a consequence of the basic postulate of the elementary interaction. As a 
next step in this discussion, we must reinterpret the Maxwell field variables 
in a way which differs from that of present-day quantum electrodynamics. 

In view of our assertion that the interaction rather than the free particle 
is elementary, we can accept the statement of the existence of a particle only 
in so far as it can influence other particles in an observable way. In particular, 
the Maxwell field variables are another way of representing a particle in an 
elementary interaction. Then, Maxwell’s equations give us a covariant means 
of determining this representation of the particle in terms of these field vari- 
ables. Based on this interpretation, we accept only the particular solutions 
of Maxwell’s equations as physically meaningful. For if there are no sources 
then there can be no field variables and vice-versa. It then follows that the 
concepts of a source-free radiation field and of a photon, as an elementary 
particle, must be abandoned. These are replaced, in the present theory, by 
the concept of the propagation of interactions between particles. Such a con- 
cept was discussed by TETRODE (7) and by WHEELER and FEYNMAN (ey 

According to our interpretation of the Maxwell field equations, the field 
variables are just another way to representing the particle’s behavior. Thus 
in an elementary interaction, the influence of one particle on another can be 


() H. Terrope: Zetis. f. Phys., 10, 317 (1922). 
(8) J. A. WHEELER and R. P. FEYNMAN: Rev. Mod. Phys, 17, 157 (1945). 
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described either in terms of the variables associated with its motion or in 
terms of the related Maxwell field variables, but not both. Consequently, such 
a concept as that of a particle acting on itself has no place in the present 
theory. (Nevertheless, when describing an elementary interaction we obtain 
formulas which are generalizations of those found in conventional electro- 
magnetic theory (°). If we relax the distinction between the particle and its 
representation by field variables, then the usual expressions are obtained.) 

The present theory is deterministic. Let us consider the nature of an ele- 
mentary interaction. No preference is given to either of the constituent 
particles; either may be the « observer » and the other the « observed ». It is 
impossible to separate the observer from the observed, and therefore a meas- 
urement on the observed (which is the function of an observer) is an integral 
part of the elementary interaction and completes its description. This requires 
a formulation of the theory which is symmetric with respect to an interchange 
of the field variables representing the observer with those of the observed. 
It should be emphasized that the completeness of the description results from 
the union of quantum mechanics with Einstein’s theory of relativity. 

Within the framework of our theory, quantum mechanics is interpreted as 
a statistical theory of elementary interactions. This is because the measuring 
process contemplated by the quantum theory is made up of a distribution of 
elementary interactions, and not of elementary particle trajectories. It is for 
this reason that we have not found it necessary to describe our field theory 
in terms of second quantization. 

In essence, our theory returns to an Einstein-Maxwell concept of a field 
theory by replacing the particle with the interaction as an elementary concept. 


3. - Maxwell’s equations in spinor form. 


It is shown by the authors (°) that the four conventional Maxwell equations 


may be expressed in terms of a pair of uncoupled two-component spinor equa- . 


tions. If we make the following identification with the conventional field 


variables . 
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then the four Maxwell equations can be expressed in the form of two uncoupled 


spinor equations 


x 


(3.2) GP x == te (x 1: 2), 
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where o,==i1 (I is the unit 2x2 matrix) and o, are the usual Pauli matrices, 

The covariance of this form of the Maxwell field equations is discussed in 
the Appendix of reference (°). 

We now specify that there is a set of Maxwell field eq. (3.2) for each con- 
stituent particle of a system. Thus, the p-th particle of an n-particle system 
is described in four-space by the source terms Y(a) or the associated solu- 
tions g(x) of the field eq. (3.2). It is emphasized that x is not a function 
of the p-th particle. The four co-ordinates (4, #1, 4, 23) are to be regarded 
as one does in a purely field theoretical description, as a set of four parameters 
which define a metric and a corresponding space which maps the field variables 
Y (x) and @,(x). It should also be emphasized that the conventional field 
variables (E, H,j,) are introduced in (3.1) only for the purposes of identi- 
cation with our field variables in one particular Lorentz-frame. The variables 
of the field theory which is developed here are w,, g, and the Dirac variable w 
(identified below with Y,). The transformation properties of the Maxwell field 
variables which leave the field eq. (3.2) Lorentz-covariant are therefore those 
of two-component (first-rank) spinors (see the Appendix in ref. (5)). Eq. (3.2) 
may then be regarded as a factorization of the vector (spin-one) field equations. 
No inconsistency is introduced by this factorization since we have abandoned 
the concept of the photon as an elementary particle. 

Finally, we must identify the source terms Y, in (3.2) with the Dirac field 
variables. To do this we note the identification of Y, with the components 
of j, in (3.1). The usual identification (for the p-th particle) is then made 
as follows: 


(p) 


(3.3) eg ess ieMy®@y, yp : 


(p) 


where e” may be regarded here as a proportionality factor. With (3.3), the 


Maxwell field eq. (3.2) take the form (for the p-th particle), 


(3.4) 5,0% PL (x) CE ey (x) I yp? (x) 7 
where 
(D) — 1 (D) mf 1 ©) 
(3.54) YO yp = Agi fe a aay | x (È 1 (1)y 
POs Li por (2)p® 


(3.5) yoL yw? = Ari 


{ po (— 191 va Va) pl? He yy 
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VI, (2)p® 

Finally, à comment should be made regarding the solutions of (3.4). In 
view of our interpretation of these equations (discussed in the preceding sec- 
tion) we accept only the particular solutions of (3.4) as physically meaningful. 


Thus, if the right-hand side of (3.4) vanishes, we must take the solution q 
a 
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to be identically zero. Similarly, if the solutions are zero, the source on the 
right-hand side of (3.4) is taken to be identically zero. 

In view of the formulation of the Maxwell field equations presented here, 
we are now in a position to construct a covariant self-consistent field theory 
in terms of the coupling of spinor fields. 


4. — A self-consistent field theory. 


It is our purpose, in this section, to construct a Lagrangian density function 
which leads, in accordance with the principle of least action, to the relativistically 
covariant self-consistent field equations for a system of two charged spin one- 
half particles. We require at the ouset that such a Lagrangian satisfy the 
requirements of Lorentz and gauge invariance, of hermiticity and be sym- 
metric with respect to the interchange of the field variables associated with 
each particle. From its variation we require a recovery of the electromagnetic 
field eq. (3.4) and a generalized Dirac equation. 


41. The principle of least action. — Since we are dealing with a e-number 
field theory, the standard application of the variational principle to classical 
fields is valid. We now specify that the two particle Lagrangian density 
(&.e., the most elementary form of the Lagrangian) depends explicitly on the 
Maxwell field spinors g®, 9%, g®, g©*, on the Dirac field bispinors y, y®, 
y®, p® and on their respective first derivatives. Thus we shall consider as 
fundamental the spinors gP, gÈ and yy and shall investigate the equations 
of motion of the system in terms of these functions. 

In the limiting case in which the electromagnetic interactions become small 
in comparison with the rest energy of the system, the terms contributing to 
the Dirac equation are obtained from the following portion of the Lagrangian 
density (*): 
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The portion of the Lagrangian density which gives rise to the Maxwell 
field eq. is expressed analogously to (4.4) and (3.4) has the following general 


form : 


x . 


2 2 
n x È ) fy ‘ € (a) | 
(4.5) AT D aes ap (a 0,02 — 2e ww y) +h. c., 
pa x=1 


q=1 


(*) The arbitrary constant that multiplies the righ-hand side of (4.4) must be, (to 
agree with experiment) fc, This is where h enters the theory. 
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f 
si where g, (a real constant with the dimensions of length) represents the one 
+ LA extra fundamental constant of our theory and where 4 and 4, are dimension- 
i ne G à qins 
Jess numbers, to be determined below. It should be noted that £,, contains 


‘an interaction term which couples the Maxwell field variables y, to the Dirac 
I" ‘ > ragtricti mn £ 
field variables through the source term w/'’,y. Note that our restriction p+ q 
follows from our discussion in Section 2 which prohibits a particle to act on 
itself. 
The remaining part of the interaction is then obtained from the part of 
the Lagrangian which couples the velocity fields of the particles. It has the 


form: 


(4.6) P,=— 2 S e MEDP@y, yp (x) [RYO x') dita’ + 
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is the Green’s function for D’Alembert’s equation, %.e., 


(4.3) LIS(x — x) =— Ô(r —r') dt — 1) (= V?— 07). 


The symbol e,,,, is the antisymmetric unit four-tensor in (uvZ0) and e, acts 
as an odd permutation operator on the particle indices (1) and (2) of the 
yp functions. The subscript I on the last integral in eq. (4.6) indicates an 
indefinite integral and the symbol d* indicates an integration over all four-space. 


The first part of £, is equivalent to the scalar interaction term 
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which appears in the conventional theory, except for the particle aspects of 
this theory, indicated by the superscripts p and q. Similarly, the second part 
of £, is equivalent (except for the p,q notation) to the pseudoscalar inter- 
St, action term 
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which appears in the generalized Dirac equation (*. The replacement of 
S(e—x') in eq. (4.6) by the right-hand side of eq. (4.7) describes the coupling 
of velocity fields in terms of the advanced and retarded solutions of the field 
equations. 

Although there exist many invariant interactions expressible in terms of 
the field variables of our theory, we have chosen the Lagrangian density La 
because it has a classical counterpart in the generalized theory. 

Summarizing, the total action function whose vanishing variation gives 
the behavior of a system of two charged particles is 


t 


x 


(4.9) I= [ (ie, + Lu+ L,) 0? x dt = | | Zarar, 
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where Y,, ,, and &, are given by eqs. (4:4), (4.5) and (4.6), respectively. 
The extension to many particles is then straightforward. 


42. The field equations. — The derivation of the field equations follows from 
the principle of least action. The vanishing of the variational derivative of 
the Lagrangian density with respect to the spinor fields (g,, qi) and (y, 7) 
give respectively the Maxwell eq. (3.4) and the Dirac equations. 

Taking the variational derivative of Z with respect to the independent 
Dirac field variable y™ (with p=1 and 2), we obtain the following coupled 
field equations: 


(4.10a) OL, 2)p? = 0, 
(4.100) O(2, ly? = 0, 


where -the operator 0(1, 2) has the form 
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and where 
(4.12) gL, =o, (1) LL) + pi(2) (2), 
(/"(1), /',(2)) are the 4 «4 matrices defined in eqs. (3.54 and b) and (g,(1), 94(2)) 
are the components of the Maxwell field spinor @,. 
(®) M. Sacus: Ann. of Phys., 6, 244 (1959). 
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The coupled set of Dirac field eq. (4.10), together with the Maxwell field 
eq. (3.4) and the interpretation of the field variables provided by the present 
theory, constitutes a completely self-consistent and relativistically covariant 
field theory of the elementary electrodynamical interaction. 

It should be noted that the non-linearity of the Dirac field equations is 
a consequence of the requirement that its solutions represent a complete de- 
scription of the elementary interaction. In contrast, the Maxwell field equa- 
tions are linear because they are descriptions of individual sources. 


43. The continuity condition. — Because of the hermitian property of the 
Hamiltonian operators associated with eq. (4.10), it follows in the usual manner 
that the following conservation law is valid (*): 


(4.13) 2, (pay wp”) = (p=1, 2). 


However, in view of our operational description of elementary interactions, 
È a > a) ne Cie is D (2) € 4 Q i a a 

the four-vectors yWy,y and yy,y” may not by themselves be interpreted 

as functions related to the observable measurements. As was pointed out 


earlier, the function which provides a complete description of the system must 
depend on both y™ and y™ and also satisfy the continuity equation 


(4.14) OP (8 yy AE less Wm) ieee ie 


It is the time component Pr (x) that weights the interaction at the space- 
time point x. The condition (4.13) follows from the requirement that the source 
of the particle force field be conserved. We now wish to consider some of the 
properties of the field variable Y as a function of y™ and y® for a special case. 


5. — The Pauli exclusion principle. 


Consider, once again, the form of the generalized Dirac equations given 
by (4.10). Multiplying eq. (4.10a) on the left by w® and the hermitian con- 
jugate of eq. (4.10b) on the right by y”, subtracting these two equations and 
then repeating the operation with (1) and (2) interchanged gives us the fol- 
lowing result 
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() This, of course, also follows from the invariance of Z with respect to gauge 
transformations of the first kind. 
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According to the Maxwell eq. (3.4), we find that 
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Therefore, if the state of motion of each of the particles (i.e., p/w) is the same, 
the right-hand side of eq. (5.4) vanishes. This implies that 


yy w(x) = poy wp (x) 


for all », so that both P, and @, vanish. 

Eq. (5.4) is just of the form of the Maxwell eq. (3.4). The reason that Q, 
vanishes is a consequence of our interpretation of these equations. It is only 
when sources are present that the field variables are non-zero. 

If we further assume that the two particles have the same inertial mass, 
eq. (5.1) reduces to 


(1) 


(GS) 0,(pey,y? + way y?) = 0. 

Combining egs. (5.5) and (4.13) and also taking into account the symmetry 
with respect to an interchange between (1) and (2) imposed by the theory, 
we obtain the following result: 

(2) 


+ (pOy,y Eye, go)]= 0, (#7, = 0, 


(5.6) 0, (poy yp Je yay wy 
where 
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The state function W(x) provides a complete description of the elementary 
interaction between two charged particles which have equal masses and are 
in the same state of motion. 

In the conventional linear theory, it would suffice to note that when the 
masses of each particle and their interactions with a background force field 
are the same, y and y'® would be solutions of the same differential equation, 


so that the linear combination (5.7) would also be a solution, The present 
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non-linear theory differs, therefore, from the conventional theory in the neces- 
sity of specifying that the solution (5.7) is valid for a system of equal mass 
particles only if the state of motion of each of the particle is the same. 
Summarizing thus far, we see that the function (y +w?)(x) satisfies the 
continuity eq. (5.6). Asserting that the interaction field variable has a unique form 
(which follows from the specification of a complete description) it is clear that 
the solution must either be (y®+wy®) or (y?—w?). The result obtained 


‘above applies to two possible cases. These are the system of equivalent par- 


ticles, each in the same state of motion, and that of the particle and anti- 
particle, each in the same state of motion. "The respective interactions are 
characterized by +e? and —e?. To resolve which of the solutions applies, 
we consider the case of two equivalent particles. The determination of Y from 


this case must necessarily apply to the second case. 
If the two particles are equivalent then 


(5.8) w (x) = (2) 


and thus, either 


(5.9a) Ys — nye ae yp == 2y 
or 
(5.90) PA NES 


The solution y (5.8) is determined from the Dirac equation 


O(p)p = 0 


which has all of the mathematical consequences of a particle acting on itself. 
The acceptance of the non-zero result (5.9a) is then incompatible with the 
basic premise of our theory which rejects the notion of a particle acting on 
itself. On the other hand, the solution (5.90) indicates that if we approach 
a situation in which the field equations describe a particle acting on itself, 
the interaction field automatically vanishes at all points in space-time. Thus, 
as a logical consequence of the fundamental premise of our theory and the 
assertion of the uniqueness of Ÿ we choose the solution TE EN A ne 
result then shows that the interaction between two identical particles, each 
in the same state of motion, is not an observable. 

The foregoing result is, therefore, in accordance (mathematically) with the 
statement of the Pauli exclusion principle. It should be emphasized, however, 
that the physical statement of the Pauli principle and our result are quite 
different because of the different interpretation of the wave function. In the 
conventional theory, y'(x,) and y®(x,) are interpreted separately as proba- 
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bility amplitudes which locate the physical particles, while in this theory, the 
state function (y'? —y™)(x) provides a complete description of the elementary 
interaction. It is the replacement of the description in terms of locations of 
particles by a description in terms of an interaction field which renders the 
present theory a field theory, in accordance with the fundamental postulates 
of Einstein’s theory of relativity. 

Since the mathematical result obtained here is equivalent to that of PAULI, 
the generalization to the many-electron system follows in the same way, and 
all consequences of the exclusion principle (such as atomic structure) are in 
complete agreement. 

The special case of the particle-antiparticle system, in which the state of 
motion of both particles is the same, is shown in the succeeding section to be 
characterized by 

pa) = — p(x) . 


This situation corresponds to the symmetric function which occurs in the 


conventional theory. 


6. — Positrorium and the pair annihilation-creation processes. 


The particle-antiparticle system is the simplest one to treat, owing to the 
fact that the Dirac function for one particle can be obtained directly from that 
of the other particle by applying the Wigner operator (!°) to the latter. Thus, 
if one of the coupled Dirac equations (describing the electron) has the form 
(4.10a), 


. (6.14) He e-I (et) + m}yp(e”) = 0 


the other equation (describing the positron) has the form 


(6.15) Mi, et F(e-) + m}p(et) = 0, 
where 
(6.2) pet) = Wyle), et=—e-=e>0, 
(6.3) Se y= ety. fr [rv e)y,p(e) S(a — x') dta’— 
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(19) E. P. WIGNER: Group Theory and its Application to Quantum Mechanics of 
Atomic Spectra (New York, 1959), chap. 26. 
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and 
FRED 0 0, —1 
a = 0 0 1 0 2 
Fr a 0 1 0 Bee 
— 1 0 0 0 


is the Wigner operator. The operator È, denotes the complex conjugation 
operation. 

Equation (6.15) (for the antiparticle) is obtained from eq. (6.1a) (for the par- 
ticle) by applying the Wigner operator to the latter. This is precisely the 
companion eq. (4.100) required by our theory for the particle-antiparticle 
system. The complete description of positronium is therefore provided by the 
self-consistent eq. (6.1) together with the Maxwell field eq. (3.4). 

We now wish to solve this set of equations for the special case in which 
the state of motion of the particles is the same. The Maxwell field eq. (3.4) 
for the electron and positron are 


(6.4a) G,0,0,,(€7) = epaler)l'uler})s 
(6,40) 0,0,Palet) = CWn (CL pn(er) , 


where the quantum numbers » and n', associated with the eigensolutions of 
the coupled Dirac eq. (6.1), specify different states of motion of the respective 
particles. If we now substitute Wy,,(e-) for w,.(e*) in eq. (6.46), we have 


(6.4e) Oyun SEO AI 
The right-hand side of eq. (6.4c) follows from eq. (6.2), since 
(6.44) Palet)y np (01) = (Wa) Poule) = 
= — (Pal) Yn Po Pel) = Papale) 
and /, depends on linear combinations of the Dirac matrices. 


It then follows that if the quantum numbers n and n’ are the same, the 
sum of eqs. (6.4a) and (6.40) is 


(6.5) COPI 0 
where 
Ha, Pa = Pal) + pale) 
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is the total electromagnetic spinor force field variable. Thus, we see that in 
this particular case the total electromagnetic field appears as a source-free 
radiation field. 

We now observe that when n= »/, 


(6.7) ET P(es) == e-F (er). 


This result follows from using eq. (6.4d) in eq. (6.3). It follows that if ple) 
is a solution of the non-linear Dirac field eq. (6.1a), then (et) =Wy(e-) is the 
solution of the coupled equation, and the desired interaction amplitude is, 
from considerations of the preceding section, 


(6.8) Y= (1—W)y(e-). 
We find below a solution for which Wy=— y so that eq. (6.8) has the form 
(6.9) Y= 2y 


which describes the ground state of the electron-positron pair. This will be 
referred to as the annihilation state, since this state of positronium will be 
shown to exhibit all of the properties usually associated with the annihilation 
of an electron-positron pair. 


61. A solution of the field equations. — A solution of the non-linear field 
eq. (6.1a) is 


exp [— imt] 


(0) 
(6.10) p= | 
0 
exp [èmt| 
if a ——a, in the interaction term (eq. (6.3)). Thus, taking a =—a =+1, 


the unique form of this interaction term is determined for all other cases. 
That (6.10) is a solution is readily verified by substituting eq. (6.10) into 
eq. (6.1a). As we show in the Appendix, with the solution (6.10), the inter- 
action term operating on yp is 


I(e)w(e) = 0 


while the interaction operator 4% itself does not vanish. On the basis of the 
present theory, this result was anticipated for the annihilation state of the 
particle-antiparticle pair. 


14 - Supplemento al Nuovo Cimento. 


803 


214 M. SACHS. and. S. L. SCHWEBEL 


Let us now consider the components of the four-vector 


(Py P (Y= (8 V) ‘y =. Den) 


where (8)? is the normalization coefficient for Ÿ in a volume V. 

According to eq. (4.14), we interpret the time component as the «interac- 
tion weighting function » for the particles to interact. From eq. (6.9) and (6.10) 
we find that 


(6.114) Py C= iV 

and that the space components are 
w 1 5 
CE a PAL 2mt , 
: LUE 

(6.116) LIT ES 7 sin Imi, 


VP PRE ON 


$ Using eq. (6.10), the source terms (3.5) of the Maxwell field equations become 


it 1 1 
epr,y = — 8nie; | 
ae e exp [2imt] 
i — exp [— 2imt] 
egI,y = saie( ; | é 


The interaction velocity field, given by eq. (6.115), indicates that the particles 


interact in a plane transverse to the 3-axis. The source terms (6.12) are de- 


scribed by two independent (right- and left-handed) circularly polarized par- 


ticle currents. Thus, the dynamical properties of the particle-antiparticle pair 
(in the annihilation state) may be compared with the conventional description 
of the annihilation process in terms of the creation of two photons. 

We now proceed to calculate the energy-momentum of this state of po- 


| sitronium. 
62. The energy-momentum of positronium in the annihilation state. — The 


invariance properties of the Lagrangian density imply conservation laws (11). 


In particular, invariance under translations gives the conservation of the 
energy-momentum vector of the system. 


(11) See, for example, N. N. BocoLiuBov and D. V. SHIRKOV: Introduction to the 
Theory of Quantized Fields (New York, 1959), Sect. 2. 
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bo 
pari 
or 


The variations of the space-time coordinates and the field variables are given, 
respectively, by the equations. 


(6.13a) DTA EUR 00 0 (x, >, =x, + de), 


and 


= N 3 / A! / > 
(6.135) SA, = 11, do, ; (A(x) > A,(x’) = A,(x) +3A4,(2)), 
where the » summation extends over the number of independent parameters ©,. 
The vanishing of the variation of the action function under the varia- 


tions (6.13) yields the conservation law 


(6.14) ,9,=0, 
where 

x 

gv 
(6.15) OE AT 1. Se 


EN E 
CPE 8) 
The independent field variables A, are the Maxwell variables q, g®* and 
the Dirac field variables y” and ym, (p,«=1, 2). 
When the co-ordinates are subjected to an infinitesimal translation 
! N Po A 16 È. RS 
o, >%,—=%,+0w,, from eq. (6.13a) we have 


1) 


SITI 


Il,,= 


(6.16) 
Substituting eq. (6.16) into eq. (6.15), we obtain the energy-momentum tensor 
(6.17) By se dA — P by, - 


From eq. (6.14), the energy-momentum vector which is conserved is 


(6.18) (pz to, dr. 


It is shown in the Appendix that, for the source terms (6.12), the Maxwell 
field eqs. (3.4) have as their solutions 


op) I 
(6.194) Dire ne | 
2190 i 4 ; 
m \exp[2imt] 
n di È 43) 
à ( Dis = . 
(6.190) Pò “a 
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With egs. (6.10) and (6.19), P, is readily evaluated from (6.15) and is found 
to vanish, (i.e., P, is a null-vector). Thus, we have shown that the solution 
(6.10) represents the ground state of positronium characterized by its energy- 
momentum vector being zero. 

The binding of the two particles in this state exhibits properties which 
conventional theory would interpret as annihilation. However, as egs. (6.19) 
indicate, positronium in the annihilation state can couple electromagnetically 
to another charged particle. With sufficient energy it is possible for the third 
particle to excite the system into a state which is interpreted conventionally 
as pair production. The result that the total energy associated with this 
(ground) state of positronium is zero means that the ground state of this atom 
is (2m) below the energy associated with the «limiting » state of the two free 
particles. It is this (ground) state of the particle-antiparticle pair which in 
our theory replaces the vacuum of the conventional theory. This is because 
it is from this state that a pair is created and into this state that a pair is 
annihilated in the conventional theory. 


7. — The hydrogen spectrum. 


We now consider the field eqs. (4.10) applied to the electron-proton system 
and show how the hydrogen spectrum evolves from their solutions. Whereas 
the application of the theory to the cases considered thus far necessitated the 
utilization of the full non-linearity of the coupled field equations, the com- 
plexity of the application to hydrogen has led us to consider a linearization 
procedure as a first approximation. This is done below by allowing the proton 
mass to become infinitely large. It is our purpose in this calculation to exhibit 
that the hydrogen spectrum (including the Lamb shift) is contained in our 
field equations. We will then calculate the Lamb shift A(n’S,—n'P,) for 
n = 2, 3,4 and determine the magnitude of the fundamental length g,,. 


71. The linearization of the field equations. — Using the same notation as 
that of the preceding section, the coupled field equations for hydrogen are 
expressed as follows: 


(7.1a) ud, CHE in a= 02 


(7.1b) {ru 0, + eI + M} y® = 0, 


where m and M are respectively the inertial masses of the electron and the 
proton and where the operator .Z is defined according to eq. (6.3). (The 
superscripts (e) and (p) refer to functions of the electron and proton solu- 
tions y° and y). 
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=] 


The approximation used by Dirac in going from his covariant wave equa- 
tion to the wave equation for hydrogen was to assume that the proton had 
infinite mass. This in turn implied that the space parts of the four-potential 
shall vanish, leaving the time component A,=ie/r. Adopting this same ap- 
proximation in our field theory, it is expected, on physical grounds, that as 
M becomes infinitely large, the space parts of the proton velocity field yoy, y™ 
Shall vanish, while the remaining time component yoyoy” should lead to the 
Coulomb potential in the corresponding part of the Hamiltonian. 

Guided by this physical argument, we consider a solution of (7.1b) to have 
the limitine form 


ILA 
0 


Md) y® = exp[—1Mit]f(r) 
Dividing (7.15) by M, we have 


1 M 1 a e dr | 
ia 7% —+ RAI È we ! il yi — (i. 


Clearly, as M approaches infinity, (7.3) takes the form 


I 
a lex Mt pe leso 
(7.4) i pe [— +Mt] f(r) Yd peer 
0) 


and (7.2) correspondingly approaches an exact solution. 
Next, it is observed that the space part f(r) of y may be chosen arbi- 
trarily, still maintaining the validity of (7.4). Once again, appealing to physical 


arguments, it is observed that taking 


f(r) 33 V(r) ; 


— 
= 
Qt 
— 


leads to the appearance of the Coulomb potential in the electron eq. (7.14). 
The square root of the 6-function means that a bilinear form of the proton 
wave function is proportional to the Dirac Ô-function. This is, in fact, the 
normalization coefficient for the proton wave function for the infinite mass 
case, and is another reason for choosing the functional form (7.5). 

With (7.2) and (7.5), we have 


(7.6) pop" — () ; poyyy = O(r) : (k =13, 0) 
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As indicated in the Appendix, we now break up % into three parts 
I= PT PT Sa. 


, : WP Ier 
Inserting (7.6) into (A.8), the interaction operator J,” takes the form 


aoe SEO 
(ot) I = er. 
This is the conventional Coulomb potential term which appears in the Dirac 
equation for hydrogen. 

In a similar fashion, the substitution of (7.6) into (A.13) gives the con- 
tribution 


(7.8) FREE [ 


This is just the added contribution coming from the pseudovector part of the 
electromagnetic four-potential in the generalized Dirac theory which was de- 
rived in an earlier publication (*). The solution of the Dirac equation for 
hydrogen with both ee and BA present in the Hamiltonian, was carried out 
exactly by the authors (12). It was found that the contribution of .4{” to the 
measured Lamb shift is indeed small and contributes to this value only be- 
cause of the presence of an external magnetic field in the experiment. It was 
assumed in this preceding work that the major portion of the Lamb shift was 
due to the radiative effects which arise in the conventional theory of quantum 
electrodynamics. The present theory, however, does not involve any such 
radiative effects and therefore the explanation for the Lamb shift must reside 
primarily in /y” (eq. (A.18)) which depends on the constant length g,,- 
Before 4,” can be determined, the two two-component spinors 7, must be 
calculated. 
Inserting (7.6) into (3.4), we have 


Dy ed ; 
0,0,p = — dried(r) fs ’ 


0 
A (p . 
oops = Axieb(r)|. } . 


The solution of these equations is carried out by the Fourier transform 


(2) M. Sacus and S. L. SCHWEBEL: Ann. of Phys., 8, 475 (1959). 
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method (18); giving 


ie ( V3 
is , 
T° \æ, + ix, 
(7.10) 
ie (È a A 
Pa = WE | . 
pe ur 
Finally, substituting (7.10) into (A.17) gives (using a, ——4a,— +1) (see Sec- 
tion (6'1)) 
Pr > Si € 
(411) ID = 167i9% 1 CO 


Combining (7.7), (7.8) and (7.11), the linearized Dirac eq. (7.14) for hydrogen 
can be written in the form 
n + ? Dr 3 be X 
(7.12) ab in BÈ +O + (Long yer) © a 


ip ) ip 


pin oe Bly 04 


where the notation used is 


NY AK = Bio-(L—er x B) +h, y=y® 
and 
7.13) B=—é&e geo p ih Cae ae 
ula TREAT Le 


i 


Following the usual procedure, we make the transformation 


UE 
where 


= se Vine ; Yi =m SE E, TUT E. 


The wave eq. (7.12) then takes the form 


(7.14) Hy= (ot Vy= iv; 
where 

4 Gg ‘di 1 SE | pm, 
(7.19) int o =i ls, | a P — À puerta puo 


(13) The details of these calculations are contained in an unpublished report by 
the authors, prepared at San Jose State College, San Jose, California. 


S 
77 


+f. +1 
pa 
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is the unperturbed (generalized) Dirac Hamiltonian for hydrogen (1?). The 


symbol y will be used to denote the fine structure constant, (e?/he). Making 
N . D 
use of our previous results (1%), the eigenvalues of 27, are 


E m jeer ni 
(7.16) ee Say ia) ee | ; 
n n (s+ n)? 
with 
Ra K?— y? 


and 7 is the half-integer angular momentum quantum number. (It is recalled 
that as £>0, Td, K +k of the ordinary Dirac. theory.) 

The eigenvalues of #, were found (12) to maintain the accidental dege- 
neracy of the excited states of hydrogen, even with +0. However, the per- 
turbing term (in dimensionless form) 


lifts this degeneracy. The strength of this term is measured by the constant 


IA = 167 9 ; 
(7.18) M LOR Vi Nos n + yi (02) yy 


where (inserting /, €) 


= hime 


is the electron Compton wave length. 

It is observed that V diverges at the origin as 07? while the ordinary Cou- 
lomb term diverges only as o-*. In order to ensure that the solutions of the 
unperturbed equation have the proper behavior at the origin to justify a 
rapidly convergent perturbation expansion of operators which depend on o ” 
(n>2), we redefine the zero order and the perturbing energy operator by 
adding and substracting the term 


10g 
0? 

to the Hamiltonian. Thus, we take 

(7.19) ik = i [%,— =) 
\ Ra 
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as the unperturbed Hamiltonian operator. The remaining part of i, 


(e <A) 


0 


a i x È 
(7.20) VARESE a ic 4 


0° 


is now considered as the perturbing term. Clearly 


The wave equation 


“ ) 
Ù al 
Hy = — - Hy . 
1 
‘an be solved exactly giving 
(7.21) y= Yn exp[—x/0] 


where w, are the eigensolutions of the original Dirac Hamiltonian %, 


(eq. (7.15). y 
The perturbation operator 7V’ (7.20) has the form 


“ 4 Weil ; Oe <d5\) 
(7.22) iV > {sin 0 exp [ig] 6 0) + è così i 4 ) ; 
where 
LU 0 
7.2: = 24 
(7.23) A i, al 
and 


is the usual Pauli matrix. 


72. The Lamb shift. — The double (accidental) degeneracy of the first 
excited state of hydrogen is described by the two orthogonal wave vectors 
|+|K|> and |—||>, corresponding to the same energy state (K is defined 


in eq. (7.16)). Since the alteration of the wave function due to E#0 is very 
slight (its influence on the calculated Lamb shift for the isolated atom is beyond 
the experimental accuracy of the measured value) (!°) we will take these eigen- 


; a a BAe ay TO ee or eer 
vectors proportional to the ordinary Dirac solutions (4) for 7=4, 3 —%%; 


(14) See, for example, H. A. BETHE and E. E. Sarperer: Quantum Mechanics of 


One- and Two- Electron Atoms (New York, 1957). 
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n=1,2,3. Using (7.21), they are, for js=-+3(°), 
\ Pm) 
0 
(7.24) \—|K n= (ile) exp [— x/o], 
VA — i cos0G_(n) 


— à sin0 exp [ip|G_(n), 
cos OF, (n) 


Ween (n/0) sin 0 exp [i@] PF, (n) 
VAri — 1G,(n) 


fa 
Ÿ 
DI 
+ 
= 
| 


exp [= x/0| 9 
0 
where |K|=1. 7F,/0 and 7G,[o are the large and small components of the 


radial part of the Dirac solution, with the appropriate + A values. 
It is readily verified that after the integration over angular variables has 


been carried out, 


(7.26) (+|E||iP|F|E|=0 
and 

ee 2 4 PACS 
(0.27) 4+|K||iV'|+|K|=+ sai “A exp [— 2x/o]do , 

0 
where, inserting h, € 
167% (or, 

1.28 n= mex = SR E) ume) 
(7.28) T yhex foo | i y?(me*) 


Thus, the energy splitting between the two states, conventionally labelled 
(n-+-1)S, and (n+1)P,, is 


(7.29 Ake SIE (GE 
where 
; P(n)G.(n 
(7.30) PETE ees Da ) exp [— 2x/o] do . 


Evaluating J, (%) for n=1, 2, 3 (corresponding to the principal quantum 


() The same results obtained below are obtained as expected for the corresponding 
states with 7, = — 1: 
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bo 
wy 


« 


apa ‘) ra ., À . . . è 
numbers 2, 3, 4) we find (to the accuracy desired in this calculation) 


(7.31a) In 
(7.315) RESTA MEZ: 07 
(7.31c) [REA RAS 


Combining (7.28), (7.29) and (7.31) we find that to order ve 


(7.324) A™ = (162/9)y4(g,,/4-)(me?) , 
(7.320) (ARES (= 0.296 ...) , 
(7.32c) AA = LS (= 0.125) . 


Experimentally, (1164) 


(7099) ZIE LOTTE 2985 


exp 
The significance of this ratio is that it is independent of the magnitude of the 
one extra fundamental constant g, that appears in our theory. The theore- 
tical value for this ratio according to the self-consistent field theory presented 
in this paper is therefore within .7°, of the experimental value. 

If we now substitute the experimental value for the Lamb shift (1°) 


ALOT ie 1079) 10 here’. 
n) 


into (7.32) we find that the fundamental length g, has the value 


(7.34) Ge 08 (a O01) «Oss cr 


Thus far we have seen that the theory predicts the two Lamb shift meas- 
urements for the hydrogen energy levels (with principal quantum numbers 2 
and 3) to a quite high degree of precision, and that these results are compatible 
with the magnitude (7.34) of the one extra constant g, inherent in our for- 
malism. 

The Lamb shift 


A (48, — 4P,) 


corresponding to that state of hydrogen with n—3 has not yet been meas- 
ured. We have calculated this quantity in order to compare it with the result 
of the conventional calculation in terms of radiative effects. 


(5) S. TREIBWASSER, E. S. Daynorr and W. E. Lamp jr.: Phys. Rev., 89, 98 (1953). 
(16) W. E. Lamp jr.-and T. M. SANDERS jr.: Phys. Rev., 103, 313 (1996). 
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With (7.32c) and (7.34) we have 
AP i= 18220 Mean 
The result obtained with the conventional theory is (*) 
A® jh = (133.10 + .02) MHz. 


Thus, while both theories seem to give the same value for the Lamb shift in 
the states with n=1 and n=2, they apparently differ by the order of 1 out 
of 133 MHz in the state with n=3. This would then indicate that if our 
linearization procedure is good te .01, an investigation of the Lamb shift 
in the higher principle quantum number states would indeed be significant. 

It is also quite interesting to note that our theory and the conventional 
one, while being so different with respect to both formalism and interpretation 
should give results for the Lamb shift which are so close. At this stage of our 
calculation, we can only conjecture that this is not pure coincidence. 

In view of our self-consistent theory, the Lamb shift for D and He* should 
be determined respectively from sets of three and five coupled non-linear equa- 
tions because of the number of particles in their respective nuclei. To treat 
these atoms as we did hydrogen would be to assume that their nuclei were 
a single particle with increased mass and charge (for Het). It is not clear at 
this stage that such an assumption is an accurate one. If, however, such an 
assumption were made, then in the approximation made thus far in this 
section, D would have the same Lamb shift as H (since À in eq. (7.29) does not 
depend on the nuclear mass). Similarly Het would have Lamb shift 16 times 
greater than hydrogen (for the n=2 state) since A « y* and in Het, y* must 
be replaced by (27) The experimental results are (4) 


Z au AGR 
ZIO ARI and qu = 14. 


8. — Discussion. 


An extension of the present work would be a study of the other states of 
positronium and of systems composed of particles of different masses. Such 
a Study would indicate whether additional terms have to be added to the inter- 
action operator /(e), since the only property of this operator needed to deter- 
mine the annihilation state of positronium is 


(8.1) SF (é) pe) = 0 . 
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For example, if in £, (eq. (4.6)), Vy, is replaced by wy,(a4 by;)y, then 
eq. (3.1) is satisfied for the same solution (6.10). Although fF, was originally 
constructed to have a counterpart in classical electromagnetic theory, such 
additional terms, without a classical counterpart, may be needed to provide 
a complete description of the electromagnetic elementary interaction. 

As mentioned in Section 7, a study of the fine structure of Het and D 
should be a very significant one, especially since the energy level structure 
of these atoms, including the Lamb shift, is known to high precision. Such 
a study would begin, in view of our theory, with an investigation of the non- 
linear equations, coupling the neutron-proton and electron to form D and of 
the four-nucleon «-particle coupled to the electron to form He*. The respective 
field equations for these systems would be three and five coupled non-linear 
equations of the form (4.10), extended to 


TRO for D, 
1,2,3,4,5 for Her. 


~ 


(8.2) Cg any eee 


Also, other interactions besides those considered in the explicit expression 
(4.11) would have to be considered in (8.2) to account for the nuclear coupling 
of the nucleons of these many particle systems. 

With regard to nuclear forces, the present theory should be extended to 
include interactions which correspond to the existence of the strong, short- 
range interaction observed in high energy nuclear experiments. It should be 
noted, however, that the various elementary particles, with their associated, 
Masses, spins, modes and rates of decay would be interpreted in our theory 
not as properties of each of the constituent particles of a system but rather 
in terms of various modes of interaction between a system of particles. 

One of the implications of the theory discussed in this paper which could 
possibly be subject to experimental corroboration by more precise measure- 
ments is the denial of the spontaneous decay process. The theory predicts 
that if a system is in a particular state, it will remain in this state unless another 
system enters, causing a transfer of energy to or from the latter. Thus, for 
example, the present theory implies that the decay of an excited nucleus occurs 
only because of its interaction with its environment. It would then appear 
that the «half-life » of a nucleus in an excited state is not truly a constant 
associated only with that nucleus. It is clear, however, that measurements 
of deviations from the constancy of the half-life of the excited nucleus should 
be extremely small due to the relatively small reaction which the macroscopie 
observer has when observing the microscopie nucleus in its environment. These 
qualitative considerations have not yet been applied quantitatively. 

The interpretation of the Maxwell field equations as a covariant represen- 
tation of the dynamics of charged particles in terms of field variables is a fun- 


Usi 


ot 
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damental tenet of our theory. The role played by the electric charge e of the 
classical theory enters here as the coupling constant ¢?, measuring the strength 
of the electromagnetic interaction. In view of this approach, the appearance 
of the inertial mass m in £,, (eq. (4.4)) is an anomaly, for it would be expected 
that m as well as e should play a role in the elementary interaction of the 
unified field theory (17). Hence, the present theory would require that the 
Einstein field equations be interpreted as a covariant representation of gravi- 
tational matter in terms of field variables. Thus, in the unified field theory, 
the inertial mass would not appear in #%,, but rather in a coupling constant 
as a measure of the strength of the elementary interaction. The latter would 
be expressed in terms of the Einstein, Maxwell, and Dirac field variables. The 
authors are currently investigating this aspect of the theory. 


APPENDIX 
The annihilation state. 


The Maxwell field eq. (3.4), with the source terms given by eq. (6.12), are 


/ 


1 
(A.la) Ou On G1 = — Errie aay 
— exp [2imt| 
_ f— exp [— 2imi]| 
(A.1b) Cao, On = SOULE | = ?,. 
1 


We proceed to solve these equations by using the Fourier transform method. Let 


(A.2) x) lot) exp [tk,æ,]dtk . 
Carrying out the indicated operations, we find that 


id, kg e E Ma f i À 
(Qn) ke? exp [— 1h,&,|(— 821) Jato 1 


— exp Zoni 


(A.3a) ilk) = — 


and 


. iogkze . (— exp [— 2imt] 
(A.3b) p(k) = — ats [exp [ 1k, | (820) | 3 1 d'a 


(7) For the solution of the two-body problem of gravitational interaction, according 
to the general theory of relativity, see A. EINSTEIN and L. INFELD: Can. Journ. Math. 
1, 209 (1949); A. EINSTEIN, L. INFELD and B. HOFFMANN: Ann. Math. 39. 66 (1938). 


A. EINSTEIN and L. InreLp: Ann. Math., 441, 797 (1940); L. INFELD: Rev. Mod., 
Phys., 29, 398 (1957). 
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where 


= yo e 
ol 0}. 


The requirement of our theory that the field variable must vanish when there 
are no sources leads to a complication which is resolved below. 
From eq. (A.3) we obtain the following solutions: 


Brie [ ra (iogk L 
(A.4a) qa) = — ——_ [exp [th p(@,— æ,)] ( 291 6) dix dig. 
(27): ki im’ 
eNO x exp [21m | 
Sie ria (— 40 gkg) (— exp [— 2umt 
(A.4b) (x) = sa [exp [ak (a, — a) | say 1 CT EH 
29 De Il 


Carrying out these integrations, we find that 


: dare 0 
(A.5a) pia) = — | Dar? 
m \exp[2imt], 
Are {exp [— 2imi] 
(A.5b) g(x) = 2 |: 
m 0 


The solutions (A.5) differ from (A.4) in that the latter satisfy the differential 
eq. (A.1) but the former do not. The apparent difficulty is associated with 
the treatment of the constant part of the source term in solutions (A.4). To 
obtain a consistent treatment of the singularities, we appeal to the invariance 
properties ot the spinors involved. 

In general, the Maxwell field eq. (3.2) can be written as follows (5): 


x = gill) gl); dea 
(A.6) Ou Oy ) mir 
P1(2) Pa(2), 7, (2) 7, (2), 


The determinant of the matrix on the right-hand side of this equation is the 
invariant (according to the spinor algebra) (f) 


Yr(i0,) Y, . 


For the particular source terms given in eq. (A.1), this invariant is zero. 
This imposes the same value on the matrix of the field variables, which is 
the invariant 

gi (ios) Ge - 


A consequence of these considerations is an alteration of the solutions (A.5) to 


dove | i 
(A.Ta) Pia) | | ; 


mM \exp[2imt] 


doe (exp [— 2imt] 
(A.71) (x) = | | , 


m | 
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which are the solutions (6.19). Thus, we define the solutions in such a way 
that the integral form (A.4) is to be used when derivatives are to be cal- 
culated; otherwise, the functions are given by eq. (A.7). 

A similar situation occurs with the velocity field coupling terms in #,: 
The scalar term gives the following contribution to I: 


n 


(A.8) I, = ty, [eww S(x — x')d4x’, 


J 


where S(a—a') (eq. (4.7)) can be written in the factored form 


1 È li 1 CALE 
(A.9) Six — x) = fare dtk (a) 


Da, a ee el UND 
“exp [ ik, (Lp — Ly) dix dik (— a exp [thy (x, — Lu) = 


A comparison of this equation with the kernel of eq. (A.4) shows S(x — x‘) 
to be equivalent to a double application of that kernel. Thus, in the event 
that any of the functions yy,y in eq. (A.8) are constant, the procedure outlined 
above for g, is doubly followed here. Carrying out this procedure for our 
specifie solution (6.10), we obtain for (A.8), 


€ 


(A.10) I, = —— {yy — ty, COS 2mt — ty, sin 2m} . 


2m 
It is readily verified that for w given by eq. (6.10), 
(A.11) (ty, Cos 2mt + vy, sin 2mt)p = yoy - 
Thus, from eq. (A.10) and (A.11), 

(A.12) Lane 0). 


The pseudoscalar velocity coupling term in £, contributes to the Dirac 
equation 


(A.13) SM i cc cers] Bray dix! fa, S(x x') a, : 


I 


To evaluate this term, we consider the indefinite integral 


ROSE 1 - exp [ik.(e,— x, 
(A.14) lè, S(x =X ) dx, = Cr) [ ar, Co vw P Le Di £ 1] dk == 
> a ve “B 
s 1 k, exp lik, (a2— om, 
SA a à Cp [ ae )] dik ; 
Cr) J ky fe 
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For the solution y given by eq. (6.10), Pyay is either a constant or a function 
of x. Therefore, the integration over the space variables can be carried out. 
We find that 


1Ee EN dix 1 
(A.15) VE = Cm E uvdo aM | Wyaw dx, | 


“ ~ 


kg (2x) - exp [thio (% — a) exp [ik ; r |: 


“S(hy) Sky) d(ks) - 


Since- 0 cannot be equal to w, it is clear that the function (k,/k,)3(k)3(k:)0(k3) 
eq. (A.15) is an odd function in one of the k, components (x=1, 2, 3). Therefore, 


(A.16) Ixy =0. 


Finally, we evaluate the interaction term 


(A.17) “ps ue — UW) y vi agi I, ty) Ti vg x) i 
x=1 


70 


Substituting the solutions (A.7) into (A.17) and using (A.11), we find that 


2 
(A.18) Fay = — 327 (“| Gud, + 2) Y3V - 
Choosing a, = — 4, we have 

SANTE AE 
Thus, 
(A.19) (4, +4, +Is)y=Iy=0. 
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